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ABSTRACT 

We identify 73 z ~ 7 and 59 z '--^ 8 candidate galaxies in the reionization epoch, and use this large 
26-29.4 AB mag sample of galaxies to derive very deep luminosity functions to < —18 AB mag and 
the star formation rate density at z ~ 7 and z ^ 8 (just 800 Myr and 650 Myr after recombination, 
respectively). The galaxy sample is derived using a sophisticated Lyman-Break technique on the 
full two-year WFC3/IR and ACS data available over the HUDF09 (--29.4 AB mag, 5cr), two nearby 
HUDF09 fields (-29 AB mag, 5a, 14 arcmin^) and the wider area ERS (-27.5 AB mag, 5a, ~40 
arcmin^). The application of strict optical non-detection criteria ensures the contamination fraction is 
kept low (just ^7% in the HUDF). This very low value includes a full assessment of the contamination 
from lower redshift sources, photometric scatter, AGN, spurious sources, low mass stars, and transients 
(e.g., SNe). From careful modeling of the selection volumes for each of our search fields we derive 
luminosity functions for galaxies ai z ^ 7 and z — 8 to < —18 AB mag. The faint-end slopes a at 
z ~ 7 and z ~ 8 are uncertain but very steep at a = —2.01 ± 0.21 and a = —1.91 ± 0.32, respectively. 
Such steep slopes contrast to the local a ~ —1.4 and may even be steeper than that at z ~ 4 where 
a = —1.73 ± 0.05. With such steep slopes [a < —1.7) lower luminosity galaxies dominate the galaxy 
luminosity density during the epoch of reionization. The star formation rate densities derived from 
these new z ~ 7 and z 8 luminosity functions are consistent with the trends found at later times 
(lower redshifts). We find reasonable consistency, with the SFR densities implied from reported stellar 
mass densities, being only ^^40% higher at z < 7. This suggests that (1) the stellar mass densities 
inferred from the Spitzer IRAC photometry are reasonably accurate and (2) that the IMF at very 
high redshift may not be very different from that at later times. 
Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

One of the primary frontiers in extragalactic astron- 
omy is observing the build-up of the galaxies in the first 
1-2 billion years of the universe. With the recent instal- 
lation of WFC3/IR on HST, we are in position to study 
galaxy evolution out to z ^ 8 (e.g., Bouwens et al. 2010b; 
McLure et al. 2010; Bunker et al. 2010; Finkelstein et al. 
2010; Trenti ct al. 2011) and even plausibly to z - 10 
(e.g., Bouwens et al. 2011a). Through the detection and 
characterization of galaxies at such high redshifts, we can 
place important constraints on the accretion rate of gas 
onto galaxies at early times, feedback effects, and per- 
haps even the dark matter power spectrum. These stud- 
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ies also provide essential inputs to reionization discus- 
sions (e.g., Bolton et al. 2007; Pawlik et al. 2009), galaxy 
modeling (e.g., Finlator et al. 2011; Lacey et al. 2011; 
Salvaterra et al. 2010; Choi & Nagamine 2010; Trenti 
et al. 2010) as well as our understanding of how quickly 
metals pollute the IGM (e.g., Ryan- Weber et al. 2009; 
Becker et al. 2009; Simcoe 2006). 

One very useful way of characterizing the galaxy popu- 
lation over an extended period of cosmic time is through 
a careful quantification of the UV luminosity function 
(LF) of galaxies with redshift. The UV LF is a partic- 
ularly useful tracer of the build-up of galaxies at early 
times, given the observed correlation of UV light with 
stellar mass (Stark et al. 2009; Gonzalez et al. 2010; 
Labbe et al. 2010a, b) and halo mass (e.g., Lee et al. 2006, 
2009; Ouchi et al. 2005; Overzier et al. 2006; McLure et 
al. 2009). Comprehensive studies of both the UV LF 
and the differential evolution of the UV LF with cosmic 
time are available at z^4-5 (Yoshida et al. 2006), z^4- 
6 (Bouwens et al. 2007: see also Beckwith et al. 2006), 
z-5-6 (McLure et al. 2009), z-4-7 (Bouwens et al. 2008), 
and z-2-4 (Reddy et al. 2008; Reddy & Steidel 2009: see 
also Sawicki & Thompson 2006). The general result of 
these studies at z < 6 has been that the volume density 
of galaxies evolves much more rapidly at the bright end 
of the UV LF than at the faint end. 

Recently, as a result of the installation of the WFC3 /IR 
camera on HST, it has been possible to extend these 
studies to z > 6. The first ultra-deep WFC3/IR obser- 
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vations over the Hubble Ultra Deep Field (Beckwith et 
al. 2006) as part of the HUDF09 program (GO 11563: 
PI Illingworth) made it possible to identify samples of 
>20 galaxies at z~7-8 (Oesch et al. 2010a; Bouwens et 
al. 2010b; McLure et al. 2010; Bunker et al. 2010; Yan et 
al. 2010; Finkelstein et al. 2010) and even out to z^8.5. 
These early WFC3/IR samples have already been used 
to derive astrophysically-interesting results on the sizes 
and structure (Oesch et al. 2010a), [/F-continuum slopes 
(Bouwens et al. 2010a; Finkelstein et al. 2010; Bunker et 
al. 2010), stellar masses (Labbe et al. 2010a, b; Finkel- 
stein et al. 2010), and specific star formation rates (sS- 
FRs: Labbe et al. 2010b) of z > 7 galaxies. 

Here we significantly improve upon these early stud- 
ies by taking advantage of the full two year observations 
over the HUDF and two other ultra-deep WFC3/IR fields 
taken as part of the HUDF09 program, as well as the deep 
wide-area (~40 arcmin^) WFC3/IR Early Release Sci- 
ence (ERS) observations (Windhorst et al. 2011). These 
observations allow us to significantly expand the sample 
of z 7 and 2; ~ 8 galaxies and substantially improve LF 
determinations. The ultra-deep HUDF09 observations 
are valuable for extending the previous results to even 
fainter levels and adding substantial statistics - allowing 
for significantly improved determinations of the LF shape 
and its faint-end slope a. Meanwhile, the wide-area ERS 
observations provide important constraints on the rarer, 
luminous population of 2^7-8 galaxies. The ERS fields 
reach '^1 mag deeper than other recent wide- area surveys 
in the near-IR (e.g., Castellano et al. 2010a, b; Hickey et 
al. 2010; Bouwens et al. 2010c; Ouchi et al. 2009). Anal- 
yses of shallower versions of this data set are provided by 
Wilkins et al. (2011), McLure et al. (2011), and Loren- 
zoni et al. (2011). 

We begin our manuscript by describing the data sets 
we use for our z ^ 7 and 2 ~ 8 selections (§2). We then 
move on to describing our procedure for constructing cat- 
alogs, selecting our Lyman Break samples, and estimat- 
ing contamination (§3). In §4, we compare our samples 
with those available in the literature and discuss the ro- 
bustness of current Lyman Break selections. In §5, we 
use the present Lyman Break samples to quantify the LF 
at z 7 and z ^ 8 and then compare our LFs with those 
in the literature (§6). In §7, we discuss the implications 
of our results for the evolution of the LF at high redshift, 
the SFR density and stellar mass density, and the reion- 
ization of the universe. Finally, in §8 we summarize our 
results. The appendices provide a detailed description of 
how we utilize the optical data to ensure our selections 
are as clean as possible, as well as compiling a number of 
other useful but secondary results and tests. Through- 
out this work, we find it convenient to quote results in 
terms of the luminosity i*^3 Steidel et al. (1999) derived 
at 2 ~ 3, i.e., M1700 ab = —21.07. We will refer to the 
HST F435W, F606W, F775W, F850LP, F098M, F105W, 
F125W, and F160W bands as B435, Veoe, 1775, Z850, ^098, 
^105, "^125, and HiQo, respectively, for simplicity. We as- 
sume ^0 = 0.3, r^A = 0.7, Ho = 70km/s/Mpc. We 
express all magnitudes in the AB system (Oke & Gunn 
1983). 

2. OBSERVATIONAL DATA 

The primary data set for this analysis is the two-year 
ultra-deep WFC3/IR observations that were taken as 



TABLE 1 

A SUMMARY OF THE OBSERVATIONAL DATA USED TO SEARCH FOR 
Z ~ 7 AND ^ ~ 8 GALAXIES. 



'*assband 


Detection Limits^''' 
(5<x) 


PSF FWHM 

(arcsec) 


Area! Coverage 
(arcmin^) 




HUDF09 (WFC3/IR HUDF) 




B435 


29.7 


0.09 


4.7 


Veo6 


30.1 


0.09 


4.7 


*775 


29.9 


0.09 


A 7 


^850 


29.4 


0.10 


4.7 


no5 


29.6 


0.15 


4.7 




29.9 


0.16 


4.7 




29.9 


0.17 


4.7 




HUDF09-1 (WFC3/IR P12) 




Veoe 


29.0 


0.09 


4.7 


i775 


29.0 


0.09 


4.7 


2^850 


29.0 


0.10 


4.7 




29.0 


0.15 


A 7 


J125 


29.3 


0.16 


4.7 


H16O 


29.1 


0.17 


4.7 




HUDF09-2 (WFC3/IR P34) 






28.8 


0.09 


3.3 




29.9 


0.09 


4.7 


«775'' 


29.3 


0.09 


4.7 




29.0 


0.09 


3.3 


^850'' 


29.2 


0.10 


4.7 


no5 


29.2 


0.15 


4.7 


J125 


29.5 


0.16 


4.7 


^^160 


29.3 


0.17 


4.7 




WFC3/IR ERS Fields (CDF-S GOODS) 


-B435 


28.2 


0.09 


39.2 


Veoe 


28.5 


0.09 


39.2 


*775 


28.0 


0.09 


39.2 


^850 


28.0 


0.10 


39.2 


5^098 


27.9 


0.15 


39.2 


J125 


28.4 


0.16 


39.2 


Him 


28.1 


0.17 


39.2 



^ 0.35"-diameter apertures 

^ No correction is made for the light outside of the 0.35"-diameter 
aperture. This is in contrast to several previous WFC3/IR studies 
by our team (e.g., Bouwens et al. 2010a, b) where our quoted depths 
were corrected for the missing light (which can result in a ~0.3 mag 
and ~0.5 mag correction to the quoted depths for the ACS and 
WFC3/IR, data, respectively, but depend upon the profile assumed). 

Our reductions of the ACS data over the HUDF09-2 field include 
both those observations taken as part of the HUDF05 (82 orbits: 
see dark blue region covering the HUDF09-2 field in Figure [TJ and 
HUDF09 (111 orbits: see cyan "PI" region covering the HUDF09- 
2 field) programs. The latter observations add ~0. 15-0.4 mag to 
the total optical depths and are important for ensuring our Lyman 
Break selections are largely free of contamination to the limits of our 
selection (~29 AB mag). 

part of the HUDF09 program (GO 11563: PI Illing- 
worth). These IR data were obtained over the HUDF 
and the two ACS HUDF05 fields. We also include results 
from the new wide area WFC3/IR Early Release Science 
(ERS) Observations over the CDF-South GOODS field 
(GO 11359: PI O'Connell). A summary of the observa- 
tional data we use for the present analysis is provided 
in Tabled The layout of these deep WFC3/IR obser- 
vations over and around the CDF-South GOODS field is 
illustrated in Figure [TJ 

2.1. HUDF09 observations 

The full two-year WFC3/IR observations from the 
HUDF09 program consist of 192 orbits of ultra-deep 
WFC3/IR data over the HUDF09 (111 orbits), HUDF09- 
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Fig. 1.— Deep WFC3/IR data over the extended CDF-South 
GOODS field used to search for 2 > 7 Lyman Break galaxies. The 
red rectangles over the HUDF, HUDF05-1 (P12), and HUDF05-2 
(P34) fields (labelled HUDF09, HUDF09-1, and HUDF09-2) show 
the position of the three ~4.7 arcmin^ ultra-deep WFC3/IR ob- 
servations (~29-29.5 AB mag at 5a) that make up the HUDF09 
program (§2.1). The upper ten rectangles show the WFC3/IR 
pointings that make up the Early Release Science observations 
(GO 11359: PI O'Connell) and cover ~40 arcmin^ (§2.2). The 
dark blue regions show the position of the ultra-deep 100-1— orbit 
ACS/WFC data over the HUDF (Beckwith et al. 2006) and 
HUDF05-1/HUDF05-2 (P12, P34: Oesch et al. 2009) fields, while 
the blue regions show the position of the ACS observations over 
the GOODS fields (Giavalisco et al. 2004). The cyan region la- 
beled PI show the position of the deep Ill-orbit ACS observations 
we obtained over the HUDF09-2 WFC3/IR field while observing 
the HUDF09 field with the WFC3/IR camera. A convenient sum- 
mary of the properties of these search fields is given in Table [T] 

1 (33 orbits), and HUDF09-2 (48 orbits) fields. The ob- 
servations from tliis program are now complete. 

All 111 orbits of imaging data over the HUDF were 
obtained on a single ~4.7 arcmin^ WFC3/IR pointing 
and were distributed over three bands Y105 (24 orbits), 
•^125 (34 orbits), and i?i6o (53 orbits). The reductions of 
those data are similar those already described in Oesch 
et al. (2010a) and Bouwens et al. (2010b), but now in- 
clude the full two years of observations of the HUDF. 
Great care was exercised to ensure the image registra- 
tion to the optical HUDF was as accurate as possible. 
Not only were latest distortion solutions (10/11/2010) 
utilized, but each source in our WFC3/IR reductions 
was cross-correlated with the corresponding source in 
the vl.O HUDF ACS zsso-band observations (Beckwith 
et al. 2006). Small corrections to the distortion solu- 
tion were required to obtain excellent registration (i.e., 
RMS differences of <0.01") for the Yios-band data (no 
WFC3/IR distortion solution were explicitly derived for 
this filter by STScI). Pixels affected by source persisi- 
tence were explicitly masked out. This masking was per- 
formed by remapping our initial reductions of the data 
back to the original frames, subtracting them from the 
original frames, coadding these subtracted frames for all 
exposures within a visit, smoothing, and then flagging 



all pixels above a 3ct threshold. The final reduced frames 
were drizzled onto the vl.O HUDF ACS reductions re- 
binned on a 0.06"-pixel frame. 

Our HUDF09 WFC3/IR reductions reach to 29.6, 29.9, 
and 29.9 AB mag (5cr: 0.35"-diameter apertures) in the 
5^105, J125, and Hieo bands, respectively. The optical 
ACS imaging over the HUDF reach to 29.7, 30.1, 29.9, 
and 29.4 AB mag (Str: 0.35"-diameter apertures) in the 
-B435, Vgoe, *775, ^850 bands, respectively. These depths 
do not include any correction for the light outside of the 
0.35"-diameter aperture and on the wings of the PSF, so 
that the depths can be readily reproduced. 

The availability of the WFC3/IR observations over the 
HUDF09-1 (33 orbits) and HUDF09-2 (48 orbits) fields 
allows us to substantially extend our 2; ~ 7 and z ~ 8 
samples. In both fields, the WFC3/IR observations were 
concentrated in single ^4.7 arcmin^ WFC3/IR pointings 
and distributed across the three bands Y105, ^125, and 
HiGo- In the HUDF09-1 field, 8 orbits in the Yios-band, 
12 orbits in the Ji25-band, and 13 orbits in the ifigo-band 
were acquired. In the HUDF09-2 field, 11 orbits in the 
Yi05-band, 18 orbits in the Ji25-band, and 19 orbits in the 
i?i6o-band were acquired. These new observations were 
reduced using the same procedures as for the HUDF. 
The WFC3/IR observations were registered and drizzled 
onto the same frame as our reductions of the ACS data 
on the two fields rebinned on a 0.06"-pixel frame. As 
with our UDF reductions, care was taken to ensure that 
the registration to the optical data was as accurate as 
possible (< 0.01") and that even issues such as velocity 
aberration were properly treated. Our reductions of the 
WFC3/IR over the HUDF09-1 field reach to 29.0, 29.3, 
and 29.1 AB mag in the Y105, J125, and Hieo bands, 
respectively. These depths arc 29.1, 29.4, and 29.3 AB 
mag, respectively, over the HUDF09-2 field (see Table[T]). 

For the ACS observations over the HUDF09-1 field, 
we could only take advantage of those from the HUDF05 
(GO10632: PI Stiavelh) program, and so we used the 
publicly available vl.O reductions of Oesch et al. (2007). 
For the HUDF09-2 field, on the other hand, deep ACS 
observations are available as a result of two programs: 
the HUDF05 and HUDF09 programs. The ACS obser- 
vations over this field from the HUDF05 program (102 
orbits in total: 9 orbits Vgoe band, 23 orbits ^775 band, 
70 orbits zgso band) were obtained over an ^^18 month 
period in 2005 and 2006 at a number of different orien- 
tations, as well as a 1' coordinate shift for 20 orbits from 
the primary pointing. Ill additional orbits of ACS data 
(10 orbits -B435 band, 23 orbits Veoe band, 23 orbits 2775 
band, 16 orbits /814 band, 39 orbits z^bq band) were ac- 
quired over this field in 2009 and 2010 from our HUDF09 
program. 

Reductions of the recent ACS observations were con- 
ducted using the ACS GTO apsis pipeline (Blakeslee et 
al. 2003). Special care was required to cope with the sig- 
nificantly reduced charge transfer efficiency in the new 
post-SM4 ACS observations and to correct for row-by- 
row banding artifacts. Our recent ACS observations 
(from HUDF09) cover -70% of the HUDF09-2 foot- 
print and add ~0. 15-0.5 mag to the depth of the earlier 
HUDF05 observations (after coaddition). Reductions of 
the HUDF05 ACS observations are described in Bouwens 
et al. (2007). The depths of the WFC3/IR and ACS ob- 
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servations over these two fields are given in Table [T] and 
reach to >29 AB mag (5cr). 

2.2. WFC3/IR ERS Observations 

The WFC3/IR ERS observations cover ~40 arcmin^ in 
the upper region of the CDF South GOODS field (Wind- 
horst et al. 2011). These observations include 10 sepa- 
rate ^4.7 arcmin^ WFC3/IR pointings (see Figure [1]). 
2 orbits of near-IR data are obtained in the F098M, 
F125W, and F160W bands, for a total of 6 orbits per 
field (60 orbits in total). These near-IR observations are 
reduced in a very similar way to the procedure used for 
the WFC3/IR observations from our HUDF09 program 
(Oesch et al. 2010a; Bouwens et al. 2010b). To keep 
the size of the drizzled WFC3/IR frames manageable, 
we split the output mosaic into 8 discrete pieces - corre- 
sponding to different segments in the ACS GOODS mo- 
saic (i.e., S14, S24, S25, S34, S35, S43, S44, 845)0 For 
each segment, we first aligned the observations against 
the ACS data binned on a 0.06"-pixel scale and then 
drizzled the data onto that frame. For the ACS data, we 
made use of our own reductions (Bouwens et al. 2006; 
Bouwens et al. 2007) of the deep GOODS ACS/WFC 
data over the GOODS fields (Giavahsco et al. 2004). 
These reductions are similar to the GOODS v2.0 reduc- 
tions, but reach ^0.1-0.3 mag deeper in the zgso-band 
due to our inclusion of the SNe follow-up data (e.g., 
Riess et al. 2007). Our reduced WFC3/IR data reach to 
27.9, 28.4, and 28.0 in the Fogs, J125 and iJieo bands, re- 
spectively (5cr: 0.35" apertures). The ACS observations 
reach to 28.2, 28.5, 28.0, and 28.0 in the B435, Vene, «775, 
Z850, respectively (5ct: 0.35" apertures). The FWHMs of 
the PSFs are ~0.16" for the WFC3/IR data and -0.10" 
for the ACS/WFC data. 

3. LYMAN BREAK SELECTION 

In this section, we describe our procedure for select- 
ing star-forming galaxies at z — 7 and z ~ 8 using the 
observational data (§2). We begin by detailing our tech- 
nique for constructing source catalogs (§3.1). We then 
describe our procedure for selecting galaxies at z ^ 7 
and z — 8 using Lyman-Break Galaxy selection criteria 
in the HUDF09 (§3.2) and ERS (§3.4) fields. In §3.5, we 
discuss possible sources of contamination and attempt 
to quantify their importance, and in §3.6, we provide a 
brief summary of our final samples of z ~ 7 and z ~ 8 
galaxies. 

3.1. Catalog Construction 

We generate separate catalogs for each of our Lyman 
Break selections, to obtain more optimal photometry for 
each selection. This is done through the use of square 
root of detection images (Szalay et al. 1999: similar 
to a coadded inverse noise-weighted image) constructed 
from only those bands that are expected to have flux for 
a given Lyman Break selection. Specifically, this image is 
constructed from the Yios, ^125 and -ffieo band images for 
our z ~ 7 HUDF09 zgso-dropout selection, the J125 and 
HiQo band images for our z — 8 IIUDF09 Yio5-dropout 

^ The configuration of these segments within 
the CDF-South GOODS mosaic is illustrat ed in 
|http: / /ar chive. stsci.edu/pub /hisp / goo ds /vl /h_cdfs_vl .Osects-plt.jpg| 




(^105 '^125)aB 



Fig. 2. — (left) zgso — Y105 vs. Y105 — J125 two-color diagram 
we use to identify 2; ~ 7 zgso-dropouts over the three ultra-deep 
WFC3/IR HUDF09 fields (14 arcmin^). The z - Y /Y - J colors 
required for our zgso-dropout selection is indicated by the gray re- 
gion. 2850-dropout candidates identified in the HUDF09, HUDF09- 
1, and HUDF09-2 fields are shown with the blue solid circles, blue 
open squares, and blue solid squares, respectively. The error bars 
and lower limits are la. The blue lines show the expected colors of 
star-forming galaxies with t/y-continuum slopes /9 of —3 and —2 
while the red lines show the expected colors of low-redshift con- 
taminants. The colors of low-mass L,T dwarf stars (e.g., Knapp 
et al. 2004) are indicated by the green squares. In addition to 
the presented two-color Lyman-Break selection, we also enforce a 
very strict optical non-detection criterion involving a Xopt quantity 
(§3.3; Appendix D). Through extensive simulations, we have found 
that a two-color LBG selection combined with strong constraints 
on the optical flux allow for the robust selection of star-forming 
galaxies at 2: ~ 7. 

selection, and from the J125 and iJieo band images for 
our ERS zgso and logs-dropout selections; the square 
root of image therefore includes all deep WFC3/IR 
observations redward of the break. 

Object detection and photometry is performed using 
the SExtractor (Bertin & Arnouts 1996) software run in 
dual image mode. Object detection is performed off the 
square root of image. Colors are measured in small 
scalable apertures (MAG_AUTO) defined using a Kron 
(1980) factor of 1.2 (where the Kron radius is established 
from the square root of image). All of our imag- 
ing data (optical/ACS and near-IR/WFC3) are PSF- 
matched to the WFC3/IR F160W imaging data before 
making these color measurements. Fluxes in these small 
scalable aperture are then corrected to total magnitudes 
in two steps. First a correction is made for the additional 
flux in a larger scalable aperture (with Kron factor of 
2.5). Then a correction is made for the light outside this 
larger scalable aperture (a 0.7"-diameter aperture is typ- 
ical) and on the wings of the PSF. The latter correction 
(typically 0.2 mag) is made based upon the encircled en- 
ergy measured outside this aperture (for stars) O The 
typical size of the total correction (including both steps) 

See 'http: // www.stsci . edu /hst / wfc3 / documents /handbooks/] 
currentlHB /cUY.irOY .html. 
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Fig. 3. — (left) Yio5 — J125 vs. J125 — Hmo two-color di- 
agram used to identify z ~ 8 Vios-dropout galaxies over the 
three ultra-deep WFC3/IR HUDF09 fields (14 arcmin^). The 
^105 — 1^125/^125 — i^l60 colors required for our yio5-dropout se- 
lection are indicated by the gray region. The other symbols and 
lines are as in Figure [3] The overlap of our color selection region 
with that occupied by L and T dwarfs is not a concern - given that 
an exceeding small fraction of the sources (<2%) in our selection 
appear to be unresolved (see §3.5). In addition to the presented 
two-color Lyman-Break selection, we also enforce a very strict op- 
tical non-detection criterion involving a Xopt quantity (§3.3; Ap- 
pendix D). As noted in Figure |2] our extensive simulations show 
that these two color criteria - combined with strict optical non- 
detection requirements - allow for the robust selection of z ~ 8 
galaxies. 

is ^0.4-0.8 mag. 

3.2. Selection Procedure (HUDF09 Fields) 

Our primary 2;~7-8 Lyman-Break samples are based 
upon the ultra-deep WFC3/IR observations over the 
three HUDF09 fields. The availability of ultra deep 
three-band (^105, >/i25, and TJieo) observations over these 
fields allow us to use traditional two-color Lyman-Break 
Galaxy selection criteria to identify star-forming galax- 
ies dX z > 7 . Spectroscopic follow-up has shown that 
these criteria are generally quite robust in identifying 
star- forming galaxies at 2;~3-7 (e.g., Steidel et al. 1996, 
2003; Popesso et al. 2009; Vanzella et al. 2009; Stark 
et al. 2010) though admittedly spectroscopic follow-up 
of the very faintest sources at z > 4 has been ambigu- 
ous since only those sources with Lya emission can be 
confirmed. 

We favor the use of Lyman-Break Galaxy (LBG) 
"dropout" selection criteria over the use of photometric 
redshifts for the identification of star-forming galaxies at 
z 7 and z 8 due to (1) the simplicity and trans- 
parency of LBG criteria, (2) the ease with which optical 
non-detection criteria can be engineered to optimally re- 
move contaminants (e.g., see Appendix D), (3) economy 
with which these criteria capture the available redshift 
information (while minimizing redshift aliasing effects), 
and (4) accuracy with which it is possible to model the 
effective selection volumes. 
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Fig. 4. — Redshift distributions predicted for our 2 ~ 7 2850- 
dropout (solid magenta lines) and 2 ~ 8 Yio5-dropout (solid red 
lines) selections over our ultra-deep HUDF09 fields (see §3.2). Also 
shown are the redshift distributions for our 2 ~ 7 2850-dropout and 
2 ~ 8 ybgg-dropout selection over the wide-area ERS fields (§3.4: 
dotted magenta and red lines, respectively). Different WFC3/IR 
y-band filters are employed for the HUDF09 fields (Vios) and for 
the ERS search area (lo98); the J125 and Higo filters are the same. 
The selection criteria are carefully chosen to allow for a better 
match between the HUDF09 and ERS redshift selection windows. 
The result is that the redshift distributions are similar at 2 ~ 7 
and 2 ~ 8. The mean redshift for our 2g5o-dropout selections is 
6.8 and 6.7 for our HUDF09 and ERS selections, respectively. The 
mean redshift for our Yio5/yo98-dropout selections is 8.0 and 7.8 
for our HUDF09 and ERS selections, respectively. 

The color criteria we utilize for our selection are sim- 
ilar to those already used by Oesch et al. (2010a) and 
Bouwens et al. (2010b). In detail, these criteria are 

(Z850 - ^105 > 0.7) A (Fios - J125 < 0.45)A 

(^850 - ^105 > 1 .4(^105 - J125) + 0.42) 

for our 2^7 zgso-dropout sample and 

(1105 - J125 > 0.45) A (J125 - Hieo < 0.5) 

for our z ~ 8 lio5-dropout selection, where A represents 
the logical AND symbol. These two-color criteria are 
illustrated in Figures [2] and |3l respectively. The criteria 
were chosen to allow for a fairly complete selection of 
star-forming galaxies at z ~ 7 and z 8 and so that our 
HUDF09 samples have a very similar redshift distribu- 
tion to those galaxies selected with the ERS yo98>^i25-f^i60 
bands (see Figure [4] and §3.4). They are also carefully 
crafted so that all star-forming galaxies at z > 6.5 are 
included and do not fall between our z 7 and z ~ 8 
selections. To ensure that sources are real, we require 
that sources be detected at >3.5i7 in the J125 band and 
3(7 in the Y105 or ffieo bands. This is equivalent to a 
detection significance of 4.5(t (when combining the J125 
and YiQz/HiQQ detections). 

For both of our z 7 zg^Q-dropout and z ^ 8 Yios- 
dropout samples, we enforce very stringent optical non- 
detection criteria. Not only do we reject sources detected 
at 2(7 in a single passband or 1.5cr in more than one band, 
but we also compute a collective xlpt value from all of 
the optical data together and eliminate sources above 
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Fig. 5. — Surface densities found in our 2~7 zgSO-dropout {upper left) and z ~ 8 ybos-dropout {lower left) selections over the wide-area 
ERS observations {black histograms, with la errors) and in our 2^7 zgso-dropout {upper right) and 2^8 Yios-dropout {lower right) 
selections over our HUDF09 observations {histograms). For comparison, we show the predicted surface densities based upon the Bouwens 
et al. (2007, 2008) determinations of the UV LF at 2 ~ 6 {cyan lines) and 2 ^ 7 {magenta lines). The solid red lines are the surface 
densities predicted extrapolating the Bouwens et al. 2008 LF results to 2 ~ 8 (based upon the fitting formula in that paper). The two 
^^^160 AB ~ 26 yio5-dropouts seen in the lower right panel - while exceeding the predicted surface densities for 2 ~ 8 star-forming galaxies 
— may be part of an overdensity (see §5.3). 



specific thresholds. We outhne this procedure in detail 
in §3.3. 

To ensure that our selections do not suffer from signif- 
icant contamination from SNe or low mass stars (both of 
which have point-like profiles), we examined each of our 
candidate z^7-8 galaxies with the SExtractor stellarity 
parameter to identify those sources consistent with being 
pointlike. The only source that we identified in our fields 
that was pointlike and satisfied our dropout criteria was 
the probable SNe (03:32:34.53, -27:47:36.0) previously 
identified over the HUDF (Oesch et al. 2010a: see also 
McLure et al. 2010; Bunker et al. 2010; Yan et al. 2010). 
No other point-like sources were found (or removed). 

We use the above criteria to search for candidate z >7 
galaxies over all areas of our HUDF09 fields where the 
WFC3/IR observations are at least half of their maxi- 
mum depth in each field (or ~4.7 arcmin^ per field). In 
total, 29 z ^ 7 zgso-dropout sources and 24 z ~ 8 Y105- 
dropouts were identified over our ultra-deep WFC3/IR 
HUDF pointing. 17 z ~ 7 zgso-dropout sources and 
14 z ^ 8 Yio5-dropouts were found over our ultra-deep 
WFC3/IR HUDF09-1 pointing, while 14 z - 7 zgso- 
dropout sources and 15 z ^ 8 Yio5-dropouts were found 
over our uhra-deep WFC3/IR HUDF09-2 pointing. As 
in the Oesch et al. (2010a) and Bouwens et al. (2010b) se- 
lections (see also McLure et al. 2010, Bunker et al. 2010, 
Yan et al. 2010, Finkelstein et al. 2010), the sources we 
identified have ffieo-band magnitudes ranging from ~26 
mag to ~29 mag. A catalog of the z ~ 7 zgso-dropouts 
in the HUDF09 fields is provided in Tables [Ml2] of Ap- 
pendix E. A similar catalog of Yio5-dropouts is provided 
in Tables \lM6\ of Appendix E. Figures [20ll2T] of Ap- 
pendix E show image cutouts of all of these candidates. 
Figure [5] shows the approximate surface density of these 



candidates as a function of magnitude in our HUDF09 
search fields. 

As a result of the very limited depths of the ACS data 
over our fields, it was necessary to exercise considerable 
caution in selecting z > 7 galaxies over our fields. Our 
conservative selection procedure resulted in our eliminat- 
ing a modest number of sources from our still very large 
z ^ 7 and z 8 samples. For comparison with other 
studies, we have included a list of these sources (which 
have a reasonable probability of being contaminants) in 
Tables [MH] of Appendix E. 

3.3. Controlling for contamination using the measured 
values of xlpt 

Identifying lower-redshift contaminants in our 
IIUDF09 selections is very challenging due to the 
limited depths of the available optical data. This is 
especially true at the faintest magnitudes. We have 
found that the only truly effective way to identify these 
contaminants is to make full use of the optical data. 

Accordingly, we have developed a rather sophisticated 
multi-step procedure to obtain tight controls on contami- 
nation. The first step is rather standard for high redshift 
selections: remove all the sources detected at 2a in a sin- 
gle optical band or at 1.5cr in more than one optical band. 
The second step is more complex, but very effective. In 
this second step we compute a collective x^pj from all of 
the optical data, as 



S^SGN(/0(/./a, 



(1) 



We reject all sources with Xopt greater than a specific 
value ^hat we establish through simulations (Ap- 

pendix D). At brighter 0.5(Ji25,ab + -ffiecAs) < 28.5 
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Fig. 6. — (left) ^850 ~ J125 vs. J125 — Hiao two-color diagram we use to identify 2: ~ 7 zgso-dropouts over the ^40 arcmin'^ ERS 
observations (first part of our two part zgso-dropout criterion). Symbols and lines are the same as Figure [2] {right) 2350 — J125 vs. 
^098 — J125 two-color diagram used to identify 2 ~ 7 2850-dropout galaxies (second part of our two part 2850-dropout criterion). Lines and 
symbols are as in the left panel. The colors required for our 2850-dropout selection are indicated in gray. We base our 2 ~ 7 2850-dropout 
selection on a 2850 — ^125 dropout criterion rather than a 2850 — 5^098 dropout criterion (e.g., as used by Wilkins et al. 2010) to extend our 
selection over a wider redshift range (see §3.4). In addition to the presented two-color Lyman- Break selection, we also enforce a very strict 
optical non-detection criterion involving a Xopt quantity (§3.3; Appendix D). 



magnitudes, we take Xum equal to 5, 3, and 5 for our 
HUDF09, HUDF09-1, and HUDF09-2 selections, respec- 
tively, while at fainter 0.5(Ji25,ab -I- ^^leo.As) > 29.2 
magnitudes, Xum taken to be half that. At magnitudes 
between those limits, the limiting Xum interpolation 
between the two extremes. More stringent limiting val- 
ues on xfim adopted for the faintest sources or where 
the optical data are shallower; more permissive limiting 
values on Xum adopted for brighter sources or where 
the optical data are deeper. Note that fi is the flux in 
band i in our smaller scalable apertures, Ui is the uncer- 
tainty in this flux, and SGN(/i) is equal to 1 if /i > and 
— f if /i < 0. The filters included in the Xopt summation 
are B435, Vqog, and 1775 for the zsso-dropout selection and 
S435, V(306, 1775, and ^850 for the Fio5-dropout selection. 

These particular limits on Xopt were chosen to minimize 
contamination for our Lyman-Break selections while 
maximizing the completeness of the z ~ 7 and 2; 8 
sources in our samples (see Appendix D for the relevant 
simulations). We determine xlpt i'^ three different sets 
of apertures (0.35"-diameter apertures, 0.18"-diameter 
apertures, and scalable Kron apertures with typical and 
maximum radii of 0.2" and 0.4") to maximize the in- 
formation we have on possible optical flux in our can- 
didates. Larger aperture measurements are useful for 
flagging more extended low-redshift galaxies in the se- 
lection. Meanwhile, the smaller aperture measurements 
provide higher S/N measurements from the higher spa- 
tial resolution ACS data. Making measurements in such 
small apertures is meaningful, given that the alignment 
we achieve between sources in the WFC3/IR and ACS 
images is better than 0.01". 

Overall, we have found that a Xopt criterion is ex- 



tremely effective at reducing the contamination rate from 
low-redshift galaxies that scatter into our selection. The 
reductions in the contamination rate are very substan- 
tial, i.e., factors of >2-3 over what it would be exclud- 
ing only those sources detected at 2a in 1 optical band 
or >1.5(T in >2 optical bands. Finally, we emphasize 
that use of the Xopt statistic and an examination of the 
Xopt distribution allows us to verify that our samples are 
not subject to significant contamination (since contami- 
nating sources would show up as a tail to positive x^ 
values: §4.2 and Appendix D.4). 



opt 



3.4. Selection Procedure (ERS Fields) 

To complement our ultra-deep searches for z~7-8 
galaxies in the HUDF09 observations, we also take ad- 
vantage of the wide area (^40 arcmin^) ERS observations 
over the upper portion of the CDF South GOODS field 
(Figure [J). 

Ideally we would use the same selection criteria for 
identifying z~7-8 galaxies over the ERS observations as 
we use over the three ultra-deep HUDF09 fields. This is 
not possible, however since the ~40 arcmin^ WFC3/IR 
ERS observations use a different F-band filter (logs) 
than in the WFC3/IR HUDF09 observations (Yios). The 
J125 and HiQQ filters remain the same. We can never- 
theless do quite well in selecting galaxies over a similar 
redshift range and with similar properties by judiciously 
choosing our selection criteria (as can be seen in Fig- 
ure E]). 

The multi-color selection criteria we use to select galax- 
ies at z ^ 7 and z ~ 8 are illustrated in Figures [B] and [71 
respectively, and nominally correspond to 2350 dropout 
and lo98 dropout selections. We base our z ~ 7 zgso- 
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dropout selection on a zgso — J125 dropout criterion rather 
than a zsso — ^098 dropout criterion to extend our selec- 
tion over a wider redshift range than is possible using a 
-2^850 ~ 5^098 criterion. Dropout selections based upon the 
latter color criterion (e.g., as used by Wilkins et al. 2010) 
result in galaxies starting to drop out of the logs band 
at z > 6.7. This makes it more difficult to identify the 
Lyman Break at 2; > 7, resulting in a redshift selection 
window that is artificially narrow. 

In detail, we require that galaxies show a strong 
•^850 — J125 > 0.9 break and satisfy zgso — J125 > 0.8 + 
l-l(Ji25 --ffieo) and zsso- J125 > 0.4+ 1.1(1098 - -^125)- 
The latter two criteria exclude sources with redshifts 
z < 6.5. To bound our redshift window selection on 
the high end to redshifts z < 7.4, we require that their 
^098 "^^125 colors be bluer than 1.25. In addition, sources 
are required to have J125 — i?i6o colors (redward of the 
break) bluer than 0.5 to exclude intrinsically red sources 
(which otherwise might contaminate our selection). Can- 
didates are required to be detected at >4.5cr in the J125 
band and at >3.5cr in the -ffieo band (to ensure that 
they correspond to real sources). A higher detection 
significance is required for our ERS selections than our 
IIUDF09 selections (§3.2) to compensate for the fewer 
exposures that go into each ERS field (and hence less 
Gaussian noise). 

For our higher 2^8 logs-dropout selections, we re- 
quire that galaxies show a strong Yogs — J125 > 1-25 
break and, similar to the previous 2; 7 zgso-dropout 
selection, have J125 — i?i60 colors bluer than 0.5. As we 
will later see, our use of a Fogs — ^125 > 1-25 criterion 
selects for star-forming galaxies at 2 > 7.4 and is conve- 
nient for combining the ERS dropout samples with the 
HUDF09 samples (selected using a slightly different F105 
filter). Candidates are again required to be detected at 
>4.5tT in the Ji25-band and Scr in the iJieo band. 

To minimize the contamination from low-redshift 
galaxies that scatter into our dropout color windows, 
we utilize very stringent criteria to exclude sources that 
show any evidence for detection in the optical. This in- 
cludes being detected at 2a in a single optical band, 1.5cr 
in more than one optical band, or having an optical Xopt 
value >2.5 (see §3.3 and Appendix D for our definition 
of xlpt)- 

To control for potential contamination from low-mass 
stars or SNe - both of which have point-like profiles - 
we examined each of our candidate 2 7 zgso or 2 8 
logs dropout galaxies with the SExtractor stellarity pa- 
rameter to identify those sources consistent with being 
pointlike. One of the sources (03:32:27.91, -27:41:04.2) 
in our selection had measured stellarity parameters con- 
sistent with being pointlike (with stellarities > 0.8 in the 
J125 and HiQo bands). This source had a J — 
color - which is much bluer than typical L* galaxies at 
2 ~ 7 (Bouwens et al. 2010a), but consistent with the 
colors of a T dwarf (e.g., Knapp et al. 2004). We there- 
fore identified this source as a probable T dwarf and re- 
moved it from our 2850-dropout selection. For all other 
sources in our selection, the measured stellarity parame- 
ters were substantially less, with typical values ^0.05-0.3 
indicative of extended sources. From a quick inspection 
of Figure [20] and [211 it is clear that essentially all the 
sources in our samples are extended. 
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Fig. 7. — ^098 — J125 vs. J125 — Hiao two-color diagram used to 
identify 2 ~ 8 Vogs-dropout galaxies over the ~ 40 arcmin^ ERS 
search area. Regions and lines as per Figure [3] In addition to 
the presented two-color Lyman-Break selection, we also enforce a 
very strict optical non-detection criterion involving a Xopt quantity 
(§3.3; Appendix D). 

The above selection criteria are only applied to those 
regions of the WFC3 /IR ERS observations over the CDF- 
South GOODS field where the optical and near-IR obser- 
vations reach within 0.2 mag of the typical ACS GOODS 
field depth and ERS WFC3/IR depth, i.e., >3 orbits 
(ACS S435 band), >2 orbits (ACS V^oe-band), >3 or- 
bits (ACS i775-band), >7 orbits (ACS zsso-band), >2 
orbits (lios), ?J2 orbits (J125), and >2 orbits (Hiqo). 
The total effective area is 39.2 arcmin^ (~3 arcmin^ of 
the WFC3/IR ERS fields extend outside of the deep ACS 
GOODS data: see Figure [1]). 

With the above criteria, we selected 13 2 ~ 7 zs5o- 
dropout galaxies and 62^82^8 logs-dropout galax- 
ies. The 2 ~ 7 Z850-dropout candidates have magnitudes 
ranging from 25.8 and 27.6 AB mag, while the 2^8 
lo98-dropout candidates have magnitudes ranging from 
26.0 to 27.2 AB mag. The properties of the sources are 
given in Tables [T3| and [ITj The surface densities of our 
2 ~ 7 2850 and 2 8 logs-dropouts are shown in Fig- 
ure [5] 

As with our HUDF09 dropout selections, we have com- 
piled a list of other possible 2 > 7 candidates that nar- 
rowly missed our ERS selection. These sources are given 
in Table [20) along with our reasons for excluding them 
from our samples. 

3.5. Contamination 

We have carefully crafted our selection criteria to op- 
timally identify star-forming galaxies at 2 ~ 7 and 2 ~ 8 
while minimizing contamination. However, some amount 
of contamination is almost inevitable, and so it is impor- 
tant to try to quantify the contamination rate. The five 
most important sources of contamination for our 2^7 
and z 8 selections are (1) low-mass stars, (2) SNe 
and other transient sources, (3) AGN, (4) lower red- 
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shift sources and photometric scatter, and (5) spurious 
sources. We consider each source of contamination in the 
paragraphs that follow. 

3.5.1. Low-Mass Stars 

Low-mass stars have 2:850 — 5^105) -2850 — "-^125 and 
-^125 — HiQo colors very similar to that of 2;~7-8 galax- 
ies, and therefore potentially contaminate our z ^ 7 
and z ^ 8 selections. Fortunately, low mass stars can 
be identified by taking advantage of the high-resolution 
WFC3/IR data, which permit us to determine which 
sources are extended and those which are not. Since 
essentially all of the sources seem to be extended ex- 
cept for one SNe candidate in the HUDF09 (§3.2) and 
one probable T dwarf in the ERS observations (§3.4), it 
seems unlikely that our sample suffers from significant 
contamination from low-mass stars (i.e., ^2%). While it 
is difficult to determine whether our faintest z > 7 can- 
didates are resolved (due to their low S/N level), the fact 
that only 2 out of our brightest ^^50 z > 7 candidates are 
consistent with being pointlike suggests that contamina- 
tion from low-mass stars is not particularly important. 
Our conclusions here are similar to those of Finkelstein 
et al. (2010: see their Figure 5). Simple arguments based 
on the surface density of stars from our galaxy at high 
latitudes (e.g., Ryan et al. 2005) also suggest that dwarf 
star contamination must be small. 

3.5.2. Transient Sources 

Since the optical ACS HUDF and GOODS data 
were taken some 5 years earlier than the new infrared 
WFC3/IR data, our samples may suffer from contami- 
nation as a result of transient sources like SNe that ap- 
pear in the new IR data but are absent in the old optical 
data. That is, they mimic high-redshift dropouts! Since 
essentially all known transient sources likely to show up 
in our searches as dropouts are unresolved, we would flag 
any such contaminants at the same time as we identify 
possible T dwarfs in the new observations. Only one such 
source is found in our search fields, and that is an appar- 
ent SNe over the ultra-deep HUDF09 pointing (Oesch et 
al. 2010a; McLure et al. 2010; Bunker et al. 2010; Yan et 
al. 2010), and it was indeed unresolved. The only other 
such source found in our selection had colors consistent 
with being a T dwarf (see §3.4). It therefore seems un- 
likely that our samples are subject to significant contam- 
ination from transient sources. 

3.5.3. Active Galactic Nuclei (AGN) 

Another possible source of contamination for our selec- 
tions could come from z > 7 AGN (e.g., Mciksin 2006). 
z > 7 AGN would show similarly strong Lyman break in 
their optical/near-IR colors and also very blue near-IR 
colors. However, we would also expect AGN to be unre- 
solved, and as we have already remarked in discussing 
possible contamination from low-mass stars and SNe, 
we only find one such source in our selection that has 
a profile consistent with being unresolved and which is 
not transient. It is interesting that obvious AGN (those 
that are comparable to, or brighter than, the underly- 
ing galaxy) are not evident (i.e., <2%) in these z ~ 7-8 
samples. 



3.5.4. Lower Redshift Sources 

Are there lower redshift sources that mimic the colors 
of z ~ 7 and z ^ 8 galaxies? They would need to be 
small but extended, faint objects that show a distinct 
and deep spectral break with rather blue colors at wave- 
lengths longward of the break. It is not at all obvious 
what such objects would be, but one possible contami- 
nating source might be a relatively-compact, low surface 
brightness, very faint (and hence low mass) z~1.5-2.5 
galaxy with a Balmer or 4000A break. Yet, whatever 
they might be, they are unlikely to dominate the num- 
bers in high redshift samples given (i) the size of the spec- 
tral breaks routinely seen in higher S/N samples and in 
stacked samples (e.g., Labbe et al. 2010a, b; Gonzalez et 
al. 2010) and (ii) the lack of large numbers of such sources 
in spectroscopic samples of z^A-Q galaxies (e.g., Vanzella 
et al. 2009; Stark et al. 2010). Instead, they likely con- 
tribute to contamination at a low level and most often 
as a result of photometric scatter. 

3.5.5. Lower Redshift Sources and Photometric Scatter 

Perhaps the largest source of contamination for our 
2 ~ 7 and z ~ 8 samples is from lower-redshift galaxies 
that satisfy our selection criteria due to the effect of noise 
on their observed colors. Clearly, this is a much greater 
problem for sources within the faintest 1-2 magnitudes 
of the sample where the effect of noise is substantial, and 
so the contamination rate is expected to be the highest 
in the magnitude interval just brightward of the selection 
limit. 

The procedure we use to estimate the importance of 
this effect relies upon actual sources from our observed 
datasets and so automatically takes into account the 
properties of "real" high-redshift galaxies. We estimate 
the number of contaminants in each of our galaxy sam- 
ples by the following multi-step process. For each faint 
source in our search fields, we (i) randomly choose a 
brighter source from our three ultra-deep HUDF09 fields, 
and (ii) then make that bright source match the faint 
source by scaling the flux of the brighter source and 
adding noise. Having done this for all of the faint sources 
in our real catalogs, we then have a "simulated" catalog 
with the same number of galaxies, magnitude distribu- 
tion, and noise characteristics. We can use this catalog 
to test the reliability of our selections by applying the 
same z ~ 7 zgso-dropout and 2; ~ 8 yios/yoos-dropout 
selection criteria as we applied on the real data. 

The brighter sources are taken from only slightly 
brighter magnitude ranges so they have higher S/N but 
are otherwise expected to be as similar as practical to 
the sources being tested for contamination. Thus, for 
our dropout samples over our three HUDF09 fields, our 
estimated contamination rates are made assiiming that 
the faint photometric samples (-ffiecAS > 27.5) have the 
same color distributions as that found in the higher S/N 
HUDF09 26.5 < Hi60,ab < 28.0 subsample. For our 
ERS selections, our contamination estimates are based 
upon the higher S/N HUDF09 26.0 < Fieo.As < 27.0 
subsample. Noise is then added to the photometry, and 
our Z85o-dropout or Y105/ loos-dropout selection criteria 
are applied. Any source that is selected by our dropout 
criteria (after adding noise), but detected in the optical 
in the original high S/N observations is counted towards 
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the total contamination level. We repeat this simulation 
200 X for each data set. 

We find 5.0 zsso-dropouts and 3.3 Yio5-dropout con- 
taminants per simulation from our three HUDF09 fields 
and 2.9 zgso-dropouts and 2.3 Yogs-dropout contaminants 
per simulation for our ERS observations. We present 
these contamination rates as a function of magnitude for 
each of our search fields in Table [9] of Appendix B. These 
simulations argue for an overall contamination rate of 
8% and 6% for our zgso-dropout and Fogs-dropout selec- 
tions, respectively, over our HUDF09 observations and 
22% and 38% for our zgso-dropout and Ybgs-dropout se- 
lections, respectively, over the ERS observations. 

Our extensive testing has given us some useful insights. 
For example, it is striking how much lower the estimated 
contamination levels are in our HUDF09 selections than 
in our HUDF09-1/HUDF09-2 selections. The differences 
are greater than a factor of ^^3 typically and are the 
direct result of the much greater depth of the optical 
data over the HUDF09 (by -1 mag) than the HUDF09- 
1/HUDF09-2 selections. The present tests also show the 
importance of the ACS parallel observations obtained as 
part of the HUDF09 WFC3/IR over the HUDF09-2 field 
(the cyan "PI" region shown on Figure [T]). Those regions 
with the deeper ACS data have an expected contamina- 
tion rate that is 2.5 x lower than those regions without 
these data (see Table [9]). 

3.5.6. Spurious Sources 

Spurious sources are unlikely to be a concern for our 
samples, given that our z > 7 candidate galaxies show 
detections in at least two bands. One of these bands is al- 
ways the Ji25-band where a >3.5(t detection is required. 
The second band is either the Y105 or -ffigo band and that 
detection must be > 3(t. Obtaining such significant de- 
tections in two independent images is extremely improb- 
able. One simple technique for estimating the probable 
contamination from such sources is to look for similar 
> 3.5(7 Ji25-band and > 3cr i?i6o-band detections on 
the negative images. No such sources were found on the 
"negative" images. 

3.5.7. Summary 

Based upon the above discussion, tests and simula- 
tions, it seems clear that the only meaningful source of 
contamination for the current selections are low-redshift 
sources that enter through photometric scatter. Over- 
all, the estimated contamination rate is 5.0 zgso-dropouts 
and 3.3 Yio5-dropouts in our HUDF09 selections and 2.9 
2;850-dropouts and 2.2 loQS-dropouts in our ERS selec- 
tions. This works out to an estimated contamination rate 
of ^11% for our z ^ 1 zgso-dropout selection and ~9% 
for our z ^ 8 Yios/Fogs-dropout selection. We there- 
fore quote a contamination rate of ^11% overall. As 
expected, most of this contamination is for the fainter 
sources. The estimated number of contaminants is tab- 
ulated as a function of magnitude in Appendix B in Ta- 
ble [S 

3.6. Summary of our z ^ 7 and z ^ 8 HUDF09+ERS 

Samples 

In total, 73 2; ~ 7 zgso-dropout candidates are identi- 
fied over all of our search fields and 59 z ~ 8 Y105 /Yogg- 
dropout candidates are found. A complete catalog of the 



TABLE 2 

Summary of z ~ 7 and z ~ 8 samples.^ 





Area 




2 ~ 7 




z ~ 8 


Sample 


(arcmin^) 


# 


Limits^ 


# 


Limits'' 


HUDF09 


4.7 


29 


J < 29.4 


24 


H < 29.4 


HUDF09-1 


4.7 


17 


J < 29.1 


14 


H < 29.0 


HUDF09-2 


4.7 


14 


J < 29.2 


15 


H < 29.0 


ERS 


39.2 


13 


J < 28.0 


6 


H < 27.5 


Total 


53.3 


73 


59 


The magnitude limit is 


the 


~5(T detection limit for ob- 



jccts in a 6.35"-diameter aperture. 



dropouts in those samples is provided in Tables [TOlfTTl 
Postage stamps of these dropout candidates are given in 
Figures dOlip 

A convenient summary of the properties of our z ^ 7 
zg5o and z 8 Yios/Yogs dropout samples and our search 
fields is provided in Table [2] The redshift distributions 
derived for the samples are given in Figure ID The mean 
redshift for our zgso-dropout selections is 6.8 and 6.7 for 
our HUDF09 and ERS selections, respectively. Mean- 
while, the mean redshift for our 1^105 /logg-dropout selec- 
tions is 8.0 and 7.8 for our HUDF09 and ERS selections, 
respectively. 

4. ASSESSMENT OF CURRENT LBG SELECTIONS 

4.1. Comparison with earlier LBG selections over our 
search fields 

Before using the present LBG selections to make in- 
ferences about the rest-frame UV LF at z > 7 and its 
evolution across cosmic time, it is instructive to compare 
these samples with previous samples over the same fields. 

4.1.1. Previous z > 6.5 Galaxy Selections within the 
HUDF09 observations over the HUDF 

Nine catalogs of z > 6.5 galaxies have already 
been published using the first-year ultra-deep HUDF09 
WFC3/IR observations over the HUDF (Oesch et al. 
2010a; Bouwens et al. 2010b; McLure et al. 2010; Bunker 
et al. 2010; Yan et al. 2010; Finkelstein et al. 2010; 
Wilkins et al. 2011; Lorenzoni et al. 2011; McLure et 
al. 2011). The z > 6.5 catalogs from this paper in- 
clude all 5 Yio5-dropout candidates in the Bouwens et al. 
(2010b) HUDF09 catalog and 14 of the 16 zgso-dropout 
candidates in the Oesch et al. (2010a) HUDF09 catalog. 
Two of the Oesch et al. (2010a) zg5o-dropout candidates 
are blended with nearby neighbors in the present catalog 
and are therefore not included here. In addition, several 
of the Z85o-dropouts in our previous catalog are now in 
our lio5-dropout catalogs, because we have modified our 
dropout criteria to better match up with the ERS filter 
set (see the discussion in §3.2 and §3.4). 

In comparison to our previous selection of z /F-dropout 
candidates over the HUDF09 (Oesch et al. 2010a; 
Bouwens et al. 2010a; Bouwens et al. 2011a), the present 
sample contains 33 additional z > 7 candidates. The ex- 
panded size of our sample is the result of our use of the 
full two-year HUDF09 observations and our decision (1) 
to extend our selection to lower S/N to include fainter 
sources in the HUDF09 observations (as our improved 
tests for contamination now allow, as discussed in §3.5 
and Appendix D) and (2) to carefully match up our z ^ 7 
z-dropout and z ~ 8 F-dropout criteria so that all cred- 
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ible 2 > 7 sources fall into one of the two samples. Com- 
pared to previous z > 6 catalogs, our dropout sample 
includes 21 additional candidates (see Tables [TOl and [T4)) . 
Essentially all of the candidates have Ji25-band magni- 
tudes of >28.5 AB mag. The lack of completely new 
candidates at brighter magnitudes is not too surprising, 
given the large number of independent z > 6.5 selections 
already performed on this HUDF09 data set. 

4.1.2. Previous z > 6.5 Galaxy Selections within the 

HUDF09-1/ HUDF09-2 Observations 

Several catalogs of z 7 zgso-dropout candidates 
are available over the HUDF09-1 and HUDF09-2 fields 
(Bouwens & Illingworth 2006; Bouwens et al. 2008; 
Wilkins et al. 2011). The first such catalogs were gen- 
erated based upon the ultra-deep NICMOS parallels to 
the HUDF (NICP12 and NICP34) and included two 
2^7 Z85o-dropout candidates NICPARl-3303-4111 and 
NICPAR2-3308-5229. Of the two candidates, only the 
first (NICPARl-3303-4111) is covered by the uUra-deep 
HUDF09 WFC3/IR observations. The source (originally 
reported at 03:33:03.81, —27:41:12.1) again makes it into 
our Z85o-dropout selection, as UDF091z-03791123. Oesch 
et al. (2009) also utihzed the uhra-deep NICMOS data 
in the NICP12 field for a z ~ 7 search, but did not report 
any candidates. 

More recently, ultra deep WFC3/IR observations of 
the HUDF09-1 and HUDF09-2 fields allow for the selec- 
tion of larger samples of z ~ 7 zsso-dropouts. Wilkins 
et al. (2011) performed exactly such a selection. Using a 
•^850 — ^105 > 1 criteria, Wilkins et al. (2011) reported 22 
candidates over the two fields to Hiqq ab ^ 28.2. While 
8 of the candidates from Wilkins et al. (2011) make it into 
our selection, strikingly most of the Wilkins et al. (2011) 
z ~ 7 candidates do not make it into our sample. While 
4 of their candidates miss our selection as a result of 
blending with nearby sources, the most common reason 
their candidates miss our selection (9 cases) is as a result 
of their candidates showing apparent detections in the 
optical i3435V606*775 data. For five of the Wilkins et al. 
(2011) candidates (P34.Z.4288, P34.Z.3053, P34.Z.3990, 
P34. z. 5016, P34.Z. 2397), the optical detections are quite 
prominent, i.e., >1.5tT in >1 band. McLure et al. (2011) 
also fiag the above sources (from Wilkins et al. 2011) as 
showing quite significant detections in the optical data. 
From the optical integration times quoted by Wilkins et 
al. (2011) for the HUDF09-2 field, it is not clear they 
take full advantage of the available ACS data. 

4.1.3. Previous z > 6.5 Galaxy Selections tn the ERS 

Observations 

The WFC3/IR ERS observations provide another 
high-quality, though much shallower dataset, so it is not 
surprising that catalogs of z ~ 7 zgso-dropouts already 
exist (Wilkins et al. 2010, 2011). Wilkins et al. (2010, 
2011) consider z ~ 7 searches over the ERS observa- 
tions using a similar Z85o-dropout criterion to what we 
use here (but see §3.4). Wilkins et al. (2010, 2011) identi- 
fied 13 zgso-dropout candidates over that area. Strikingly 
enough, only two of their thirteen candidates are present 
in our zgso-dropout catalogs (see Table [T^ . 

Most of their sample (11 sources) do not make it 
into our zgso-dropout selection. What is the reason 
for this? Three of their candidates (zD3, ERS.z.26813, 



zDl) miss our color selection due to their having 
2850 — ^098 colors that are too blue (e.g., ERS.z.26813 
is given in Table [20|) . but seem likely to correspond 
to z > 6 galaxies. One other candidate (zD2) was 
blended with a nearby source and therefore not se- 
lected. Another (ERS. z. 87326) was too faint in the 
-ffieo band to be selected. The remaining six z ^ 7 
candidates from Wilkins et al. (2010, 2011) were not 
included in our selection since they seem much more 
likely correspond to low-redshift galaxies or stars. Five 
(zD5,ERS.z.45846,ERS.z.20851,ERS.z.47667,ERS.z.80252) 
show significant detections in the optical (> 1.5cr in 
~2-3 bands) and the sixth (ERS.z. 70546: 03:32:27.91, 
—27:41:04.2) is almost certainly a T-dwarf, given its 
smaU size (~0.1" half-light radius in the WFC3/IR 
data), SExtractor stellarity parameter (^0.95: strongly 
suggesting the source is unresolved) , and extremely blue 
^125 - -H^ieo color (- -0.6). 

4.2. Integrity of Our Lyman-Break Galaxy Selections 

One of the most significant challenges in constructing 
large samples of z > 7 galaxies to very faint magnitudes 
is the limited depth of the optical data. Without deep 
optical data, it is very difficult to discriminate between 
bona-fide high-redshift galaxies and z^l-2 galaxies that 
are simply faint in the optical. How confident are we that 
we have identified all the contaminants in our selections 
(or conversely that we have not discarded a substantial 
number of bona-fide z > 7 galaxies)? 

There are two arguments that suggest the situation 
is largely under control. First, carefully modeling the 
selection efficiency for each of our dropout samples and 
fixing the shape of the LF, we can estimate (p* for the 
z ~ 7 and z ^ 8 LFs on a field-by-field basis (§6.2). If 
we were misidentifying too many low-redshift galaxies as 
high-redshift galaxies - or vice- versa - we would expect 
to observe huge fluctuations in the value of 0* in precisely 
those fields with the shallowest optical data. However, 
we find essentially the same value of (f>* (within 20%) for 
each of the three ultra-deep fields we consider (HUDF09, 
HUDF09-1, HUDF09-2). This suggests that we have a 
reasonable handle on contamination in our samples. 

Second, we can look at the extent to which our zgso- 
dropout and Yios/logs-dropout samples show detections 
in the optical data (quantified here using a xlpt statistic). 
If our selections were largely composed of bona-fide z^7- 
8 galaxies, then we would expect there to be no more 
sources with positive detections in the optical data than 
there are negative detections. Indeed, our samples show 
a similar number of sources with negative and positive 
detections. See Appendix D.4 for more details. 

It is worthwhile noting that much of this discussion 
will likely be highly relevant to future z > 7 selections 
with the James Webb Space Telescope ~ where the lack 
of visible imaging of comparable depths may make the 
selection of z > 7 galaxies challenging. 

5. CONSTRAINTS ON THE REST-FRAME UV LFS AT Z~7 
AND 

The large samples of z~7 and Z'^8 galaxies we iden- 
tified over the wide-area ERS and ultra-deep HUDF09 
fields permit us to place powerful constraints on the 
shape of the UV LF at z~7 and z~8. The available 
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TABLE 3 

Stepwise Constraints on the 
rest-frame uv lf at z ~ 7 and 
2 ~ 8 (§5.1 and §5.3).^ 



Muv.Ab" 4>k (Mpc-3 mag-i) 





2- 


-dropouts (2 ~ 7) 


-21, 


.36 


<0.00002'" 


-20, 


.80 


0.00016±0.00009 


-20, 


.30 


0.00015±0.00009 


-19, 


.80 


0.00050±0.00021 


-19, 


.30 


0.00104±0.00035 


-18, 


.80 


0.00234±0.00068 


-18, 


.30 


0.00340±0.00096 


-17, 


.80 


0.00676±0.00207 




y-dropouts (z ~ 8) 


-21, 


.55 


<0.0000l'" 


-20, 


.74 


0.00011±0.00007 


-20, 


.14 


0.00025±0.00012 


-19, 


.54 


0.00039±0.00017 


-18, 


.94 


0.00103±0.00035 


-18, 


.34 


0.00156±0.00072 


-17, 


.74 


0.00452±0.00207 



^ These stepwise LFs are also 

shown in Figures [8] and [9] 

^ Upper limits here are Icr (68% 

confidence). 

The effective rest-frame wave- 
length is ^ I6OOA for our 2 ~ 7 
selection, ~ 1760A for our 2 ~ 8 
selection. 



constraints are now much stronger than was possible us- 
ing just the early WFC3/IR data over the HUDF (e.g., 
Oesch et al. 2010a; McLure et al. 2010). 

We begin this section by establishing the UV LF at z ^ 
7 considering both a stepwise (§5.1) and Schechter (§5.2) 
representation. In §5.3, we move onto a determination of 
the LF at 2 8. 

5.1. z ^ 7 LF: Stepwise Determinations 

We first consider a stepwise determination of the UV 
LF at z 7. Stepwise determinations of the UV LF 
are valuable since they allow for a relatively model- 
independent determination of the shape of the LF. This 
ensures that our determinations are not biased to adhere 
to specific functional forms such as a Schechter function 
or power law. 

Here the stepwise luminosity function (j)^ is derived 
using a very similar procedure to the Efstathiou et al. 
(1988) stepwise maximum likelihood (SWML) method. 
With this method, the goal is to find the shape of the LF 
which is most likely given the observed distribution of 
magnitudes in our search fields. Since only the shape of 
the distribution is considered in this approach (and there 
is no dependence on the volume density of sources), we 
would expect the LF to be largely insensitive to large- 
scale structure effects and hence more robust. Previously, 
we employed this procedure to derive the LFs at 
(Bouwens et al. 2007). 

In detail, one can write this likelihood as 

C = \lfieid^iP{mi) (2) 

where 

p(m,)=U^^=^ (3) 

\ ^j'^expectedj / 



TABLE 4 

Other Wide-Area Searches for 2 > 7 Galaxies Considered 

Here. 





Area 


Redshift 




Volume 


Instrument 


(arcmin^) 


Range 


< 2 > 


(Mpc^) 


NICMOS + 


~80 


6.3-8.5 


7.1 


1x10^ 


MOIRCS/ISAAC 


248 


6.5-7.5 


7.0 


4x10^ 


Suprime-Cam 


1568 


6.5-7.1 


6.8 


7x10^ 


HAWK-I 


~151 


6.5-7.5 


6.8 


1x10^ 



References: [1] Bouwens et al. 2010c, [2] Ouchi et al. 2009, and [3] 
Castellano et al. 2010a, b 

where nobserved,i is the observed number of sources in 
the magnitude interval i, Uexpectedj is the expected num- 
ber of sources in the magnitude interval j, and where 11 
is the product symbol. The above expression gives the 
likelihood that a given magnitude distribution of sources 
is observed given a model LF. The expected number of 
sources riexpected.i is computed from a model LF as 

^k(t>kVj^k — nexpectedj (4) 

where riexpectedj is the surface density of galaxies in some 
search with magnitude j, 0^ is the volume density of 
galaxies with absolute magnitude fc, and where Vij^k is 
the effective volume over which galaxies in the absolute 
magnitude interval k are selected with a measured magni- 
tude in the interval j. With the Vij^k factors, we implic- 
itly account for the effects of photometric scatter (noise) 
in the measured magnitudes on the derived LFs - allow- 
ing us to account for effects like the Malmquist bias on 
our LF determinations. The Malmquist bias is known 
to boost the number of sources near the selection limit, 
simply as a result of the lower significance sources (^5(t) 
scattering into the selection. See Appendix B of Bouwens 
et al. 2007 (and Bouwens et al. 2008) for a more detailed 
description of the above maximum likelihood procedure. 
We elected to use magnitude intervals of width 0.5 mag 
for our LF determination as a compromise between S/N 
(i.e., the number of sources in each interval) and resolu- 
tion in luminosity. 

We estimate the selection volumes Vj^k for our LBG 
selections using the same set of simulations described in 
Appendix C of this paper. This procedure is the same as 
used in our many previous LF studies (e.g., Bouwens et 
al. 2007, 2008, 2010b) and involves taking the pixel-by- 
pixel profiles of actual HUDF z ~ 4 i3-dropouts, insert- 
ing them at random positions in the actual observations, 
and then attempting to select them using the same pro- 
cedures as we use on the real data (including the Xopt 
non-detection criterion we describe in §3.3 and Appendix 
D). The J/y-continuum slopes /3 adopted for the z 7-8 
model galaxies are chosen to match the observed color 
distribution (Bouwens et al. 2010a). The selection vol- 
umes Vij^k we derived from these simulations have a very 
similar form to those shown in Figure 8 of Bouwens et al. 
(2006) or Figure A2 of Bouwens et al. (2007). We would 
expect these volumes to be reasonably accurate - since 
we reproduce both the observed sizes of z ~ 7-8 galaxies 
(Appendix A) and the C/T^-continuum slopes. 

Due to its relative insensitivity to large scale structure 
effects, the stepwise maximum likelihood procedure de- 
scribed above only provides us with constraints on the 
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TABLE 5 

Determinations of the best-fit Schechter Parameters for 

THE rest-frame UV LFS AT 2 ~ 7 AND 2 ~ 8 AND AT 2~4, 5, 6 

from Bouwens ft al. (2007). 



Dropout 
Sample 


Rcdshift 




0* (10~3 
Mpc-3) 


a 


z 
yb 


6.8 
8.0 


-20.14 ±0.26 
-20.10 ±0.52 


0.861°:^° 
0.591^:°^, 


-2.01 ±0.21 
-1.91 ±0.32 


B 
V 

i 


3.8 
5.0 
5.9 


-20.98 ±0.10 
-20.64 ±0.13 
-20.24 ±0.19 


1.3 ±0.2 
1.0 ±0.3 

1 4+0-6 


-1.73 ±0.05 
-1.66 ±0.09 
-1.74 ±0.16 



Values of M*y are at 1600 A for the Bouwens et al. (2007) z ~ 4, 
2 ~ 5, and z ~ 7 LPs, at ~ 1350 A for the Bouwens et al. (2007) 
2 ~ 6 LP, and at ~ 1750 A for our constraints on the 2 ~ 8 
LP. Since 2~6-8 galaxies are blue {/3 ~ —2: Stanway et al. 2005; 
Bouwens et al. 2006, 2009, 2010a), we expect the value of M*y to 
be very similar (< 0.1 mag) at 1600 A to its value quoted here. 

The derived Schechter parameters depend significantly on 
whether we include or exclude the two brightest (~26 AB mag) 
z ~ 8 galaxies in the HUDP09-2 field. If we exclude the two bright- 
est galaxies (assuming they represent a rare overdensity ) , we derive 
a significantly fainter characteristic luminosity M* and larger nor- 
malization <f>*: M* = -19.54 ±0.56, 4>* = l-btfl ^ Mpc'^ 
o = -1.67 ± 0.40. 



shape of the LF. To provide a normalization for the LF, 
we require that the total number of z ~ 7-8 galaxies 
we predict from Eq. 2] in our search fields matches those 
we actually find (after correcting for contamination: see 
Table [9] and Appendix B for details on how the precise 
corrections depend on search field and the luminosity of 
the sources). 

By applying the above maximum likelihood procedure 
to our ERS, HUDF09, HUDF09-1, and HUDF09-2 z - 7 
samples, we derive the stepwise LF at z ~ 7. The result- 
ing LF is presented in Table [3] and Figure [8] (red circles). 
The z ~ 7 LFs inferred from our z ~ 7 ERS samples 
alone are shown separately {black squares) to demon- 
strate the consistency of our derived LFs in both the 
wide-area and ultra-deep data. Also included in Figure |S] 
are the LF determinations at z ~ 7 from two wide-area 
searches (Ouchi et al. 2009; Bouwens et al. 2010c: see 
Table 2]). These searches are important in establishing 
the shape of the LF at very high luminosities (< —21 AB 
mag). 

5.2. z ~ 7 LF: Representation with a Schechter 
Parameterization 

In the present section, we consider representations 
of the z ~ 7 LF with a Schechter parameterization 
(r(ln(10)/2.5)10-0-4(^^-A^*)"e-io-''-*<"^-'-^''). This pa- 
rameterization features an exponential cut-off at high 
luminosities and a power-law shape at fainter luminosi- 
ties. Such a parametrization has been almost universal in 
characterizing the shape of the galaxy LF since it works 
so well. From the results of the previous section (e.g.. 
Figure [S]), it also appears to be relevant in describing the 
LF of z ^ 7 galaxies. 

As with our stepwise LF determinations, we use a max- 
imum likelihood approach that considers only the shape 
of the LF in deriving the best-fit Schechter parameters. 
The approach is analogous to that developed by Sandage 
et al. (1979), but is formulated in terms of the apparent 
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Fig. 8. — (upper) Rest-frame UV LP determination at z ~ 7 
determined from the present WPC3/IR z ~ 7 samples (§5.1). The 
red circles, red downward arrow, and red line show our stepwise LP 
determination, the Icr upper limit we have on the bright end of the 
LP from our WFC3/IR z ~ 7 searches, and best-fit Schechter LP, 
respectively. Also included in the upper panel are the LPs derived 
from a wide-area Subaru Suprimc-Cam search (Ouchi et al. 2009), 
a wide-area NICMOS+ISAAC+MOIRCS search (Bouwens et al. 
2010c: see also Mannucci et al. 2007; Stanway et al. 2008; Henry 
et al. 2009), and a wide-area HAWK-I search (Castellano et al. 
2010b). The black squares show the 2 ~ 7 LP constraints from 
the wide-area ERS observations. These match up well with the LP 
constraints obtained from the deeper HUDF09 observations. The 
dotted red lines show the approximate asymptotic behavior (linear 
fit) of the 2 ~ 7 LF at the bright and faint ends. Together they 
illustrate the existence of a likely break in the LF at —20 AB mag. 
(lower) 68% and 95% confidence intervals on the characteristic 
luminosity A/* and faint-end slope a from our WPC3/IR samples 
(dotted black lines: §5.2). Stronger constraints on M* and a can 
be obtained by incorporating wide-area z ~ 7 search results (solid 
red lines). 

magnitude distribution (see §3.1 of Bouwens et al. 2007). 
The approach has the important advantage that it is al- 
most entirely insensitive to large-scale structure effects 
(see Appendix C of Bouwens et al. 2007). 

To perform this likelihood analysis, we start with var- 
ious Schechter parameter combinations, calculate the 
equivalent stepwise LF i/i^'s (adopting 0.1 mag bins), and 
then make use of Eq.[2]and|4l This technique allows us to 
set constraints on the characteristic luminosity M* and 
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Fig. 9.— {upper) Rest-frame UV LF at z ~ 8 (see §5.3). The 
red circles show the stepwise LF for the WFC3/IR HUDF09+ERS 
2 ~ 8 sample. The red downward arrow shows the Icr upper limit 
on the bright end of the LF from our WFC3/IR z ~ 8 search. The 
black points show the stepwise LF at 2 ~ 8 derived from the z ~ 8 
ERS search. The dashed black line shows the best-fit power-law 
representation while the red line is for the Schechter fit. The faint- 
end slope a for the 2 ~ 7 LF (a = —2.01) is indicated on this 
figure with the dotted red line. The dotted blue line shows the 
approximate cut-off observed at the bright-end of z < 7 LFs. The 
green upper limit (at —22.25 mag) is based upon a recent wide- 
area z ^ 8 Y-dropout HAWK-I search (Castellano et al. 2010b). 
See also the caption to Figure [S] Note that the brightest points in 
our 2 ~ 8 LF are significantly affected by the three bright sources 
in the HUDF09-2 field (compare the brightest red point with the 
black square: based upon the ERS search), {lower) 68% and 95% 
maximum likelihood contours on the characteristic luminosity M* 
and the faint-end slope a inferred for the LF at 2 ~ 8. 

faint-end slope a. The lower panel of Figure [8] [dotted 
black lines) shows the 68% and 95% confidence intervals 
we are able to obtain on the z ~ 7 LF. 

Inspecting this lower panel, we see that there is a large 
degree of freedom in the Schechter parameters M* and 
a allowed by our search results. The principal reason 
for this latitude (which occurs despite the large number 
of sources and large luminosity range) lies in the almost 
featureless power-law shape exhibited by the LF at z ~ 7 
(top panel: Figure [5]). This restricts us to characteristic 
luminosities M* brighter than ~ —19.7 and approximate 
power-law slopes a of —2.0. 



To obtain much tighter constraints on the Schechter 
parameters a and M*, we need to incorporate observa- 
tions which allow us to resolve more distinct features in 
the LF (e.g., the expected break at brighter magnitudes 
that is not apparent in our small area HUDF09-I-ERS 
data sets). Such constraints can be obtained by combin- 
ing our WFC3/IR results with wide-area searches (e.g., 
Ouchi et al. 2009; Castellano et al. 2010a, b; Bouwens et 
al. 2010c) which more effectively probe the rarer sources 
found at the bright end of the LF. Such searches have 
consistently found a distinct fall-off in the volume den- 
sity of z 7 sources at bright magnitudes (e.g., upper 
panel in Figure H]). Several of the most notable wide-area 
z ^ 7 searches (Table |4|) include a ~85 -t- 248 arcmin^ 
search with NICMOS and ISAAC/MOIRCS by Bouwens 
et al. (2010c: see also Mannucci et al. 2007; Stanway et 
al. 2008; Henry et al. 2009), a '^150 arcmin^ search with 
VLT HAWK-i by Castellano et al. (2010b: but see also 
Castellano et al. 2010a and Hickey et al. 2010), and a 
1568 arcmin^ search with Subaru Suprime-Cam by Ouchi 
et al. (2009). The mean redshift expected for these sam- 
ples is typically ~6. 8-7.0, which is comparable to that for 
the present WFC3/IR zsso-dropout selections. 

Combining the wide-area z ^ 7 search results with 
those available from the ultra-deep WFC3/IR obser- 
vations, we arrive at much tighter constraints on M* 
and a. The best results are M* = -20.14 ± 0.26 and 
a = —2.01 ± 0.21 and are presented in the lower panel to 
Figure [8] with the solid red lines. In combining the wide- 
area search results with the WFC3 /IR results (where the 
overlap in identified z ^ 7 candidates is minimal), we find 
it convenient to modify our LF estimation methodology 
so that the likelihood incorporates information on the 
volume density of the LF (and not simply its shape) - 
comparing the volume density of the sources found with 
what is expected from a model LF. As in our Bouwens et 
al. (2008) LF analysis, likelihoods are computed assum- 
ing Poissonian statistics. We verified that the 68% and 
95% likelihood contours we compute for the WFC3/IR 
search results reproduced those calculated using the STY 
method. 

It is trivial for us to extend the present analysis to 
constrain the normalization (p* of the Schechter function. 
Utilizing the wide-area and ultra-deep WFC3 search re- 
sults, we find a (/)* of 0.86^°;^° x 10"^ Mpc"^. Together 
with the best-fit values of a and M*, we conveniently 
include these parameters in Tabled The effects of large- 
scale structure add somewhat to the uncertainties on 0* . 
Given that most of the sources which contribute to the 
faint counts in the LF are from the CDF South, we es- 
timate that the added uncertainty is 25% on (j)* - which 
corresponds to the sample variance expected assuming a 
'^lOO arcmin^ search area (the approximate area within 
the CDF South where our fields are found), a redshift 
selection volume of width Az ^ 0.8 (from Figured]), and 
a bias factor of ~5 (corresponding to halo volume densi- 
ties of ~10~'^ Mpc~'^). These estimates were made using 
the Trenti & Stiavelli (2008) cosmic variance calculator 
(see also §6.2). 

Of course, the volume density (f>* in the Schechter pa- 
rameterization shows a significant covariance with both 
the characteristic luminosity and faint-end slope a, and 
so there is considerable freedom in the Schechter param- 
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Fig. 10. — Rest-frame UV LF determinations for galaxies at 2: ~ 7 
(§6.1.1). Shown are the UV LF from this paper {red circles and 
red downward arrow) and from Bouwens et al. (2010c: open red 
circles), Oesch et al. (2010a: purple circles), Ouchi et al. (2009: 
black squares), Castellano et al. (2010b: green squares), McLure et 
al. (2010: blue squares) and Bouwens et al. (2008: black crosses). 
Constraints from Wilkins et al. (2011) are similar, but are not 
shown to reduce confusion. See Table [6] for a comparison of the 
Schechter parameters derived here for the z ~ 7 UV LF with that 
found in other analyses. 

eters that can fit current search results. We present the 
likeUhood contours on </>* and the other two Schechter 
parameters in §7.1. 

5.3. z r^8 LF 

We now derive the LF at z ~ 8 based upon the z ~ 8 
selections presented in §3.2 and §3.4. As expected, we 
will find that the z ^ 8 LF is less well-determined than 
that at z 7, but it is still striking to find very in- 
teresting constraints from the combined HUDF09-t-ERS 
datasets. 

We use the same procedure for deriving the stepwise 
LF and Schechter parameters at z ^ 8 as we had earlier 
used at z ~ 7. We base our LF determinations on our 
combined HUDF09 -I- ERS z ~ 8 selections, but also 
fold in the null results from a shallow, ~150 arcmin^ Y- 
dropout search over the CDF-South, BDF, and NTTDF 
(Castellano et al. 2010b). Finally, to account for the 
fact that the mean rcdshift of our wide-area ERS z ~ 8 
selection is Az ~ 0.2 lower than for our HUDF09 z 
8 selection, we adjust the volume densities implied by 
the ERS search downward by a factor of 1.15 (which 
we derive from a fit to the stepwise UV LF evolution 
from z ~ 6 to z ~ 4). We use 0.6 mag intervals for the 
stepwise LFs at z ^ 8 because of the fewer sources in 
~ 8 samples. 

Our results are presented in Table |3] and Figure |9] (red 
points and upper limits: both the HUDF09 and the ERS 
observations are incorporated into the red circles). Also 
included on this figure are our z ~ 8 LF using only the 
ERS observations. Not surprisingly, given the general 
form of the LFs at z ~ 4-7, the z ~ 8 LF maintains an 
approximate power-law form over a ~3 magnitude range 
(red points), with a steep faint-end slope. Brightward of 
—21 AB mag, the LF appears to cut off somewhat. This 
cut-off is clear not only in the upper limits we derive for 
the volume density of luminous sources (red and green 
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Fig. 11. — Rest-frame UV LF determinations for galaxies at 
2 ~ 8 (§6.1.2). Shown are the UV LFs from this paper {solid 
red circles and red upper limits) and from Bouwens et al. (2010b: 
black open circles and black upper limits) and McLure et al. (2010: 
blue squares). LF estimates from Lorenzoni et al. (2011) are sim- 
ilar. Also presented is the la upper limit from Castellano et al. 
(2010b: green downward arrow). See Table[7]for a comparison of 
the Schechter parameters derived here for the 2 ~ 7 UV LF with 
that found in other analyses. 

downward arrows from WFC3/IR and HAWK-I searches, 
respectively), but also in the z 8 ERS results (black 
squares). 

Nonetheless, sufficiently small area has been probed at 
present that the bright end of the z 8 LF remains 
fairly uncertain, and indeed there is significant field-to- 
field variation across our fields. The most notable case 
is over the HUDF09-2 field where 4 bright i/ieo.AS < 27 
z ^ 8 galaxy candidates are found, 2 lying within 3" 
of each other (UDF092y-03781204, UDF092y-03751196) 
and a third Hiqo ^ 27 candidate within 30" (UDF092y- 
04640529: see Table [TB)) . each of which have similar 
^105 ~ ^125 colors and presumably redshifts. As a re- 
sult of such variations, the bright end of the z 8 LF 
is understandably still quite uncertain, with some differ- 
ences between that derived from the ERS observations 
alone (black squares) and from the combined HUDF09 
+ ERS results (red circles). 

Despite the still somewhat sizeable uncertainties, we 
can model the z ~ 8 LF using a Schechter function and 
derive constraints on the relevant parameters. Formal 
fits yield best-fit Schechter parameters of a = — 1.91± 
0.32, M* = -20.10 ± 0.52, and 4>* = O.bdtlfr x lO'^ 
Mpc~'^ for the z ~ 8 LFQ Though the significance is 
only modest, Schechter fits to the z 8 results (red line 
on Figure [5]) are preferred at 2a over fits to a power law 
[dashed line). The implications of such steep faint-end 
slopes are discussed below, but we caution against giv- 
ing too much weight to these slopes at this time, given 
the current large uncertainties. Table [5] gives a compila- 
tion of these parameters, along with a comparison with 

The derived Schechter parameters depend significantly on 
whether we include or exclude the two brightest (~26 AB mag) 
z ~ 8 galaxies in the HUDF09-2 field. If we exclude the two 
brightest galaxies (assuming they represent a rare overdensity) , 
we derive a significantly fainter characteristic luminosity M* and 
larger normalization;/)* : M* = -19.54±0.56, </>* = 1.5tj Q X 10"^ 
Mpc-3, a = -1.67 ±0.40. 
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TABLE 6 

Determinations of the best-fit Schechter Parameters for the 
rest-frame uv lfs at 2 ~ 7 (§6.1.1: see also figure [toll. 

Reference ^^irv 't'* (10~^ Mpc~^) a 

This work -20.14 ±0.26 Omtolg -2.01 ± 0.21 

Castellano et al. (2010a) -20.24 ±0.45 O-SSIqii -1.71 (fixed) 

McLure et al. (2010) -20.04 0.7' -1.71 (fixed) 

Oesch et al. (2010a) -19.91 ±0.09 1.4 (fixed) -1.77 ±0.20 

Ouchi et al. (2009) -20.1 ±0.76 O-^^tl tl -1.72 ±0.65 

Oesch et al. (2009) -19.77 ±0.30 1.4 (fixed) -1.74 (fixed) 

Bouwens et al. (2008) -19.8 ±0.4 l-lto7 -1.74 (fixed) 



Schechter determinations at other redshifts. The trends 
of the Schechter parameters with redshift are also dis- 
cussed more broadly below. 

6. ROBUSTNESS OF LF RESULTS 

6.1. Comparison with previous LF determinations 

6.1.1. LF 

A large number of different determinations of the UV 
LFs at z ~ 7 now exist (Bouwens et al. 2008; Oesch et 
al. 2009; Ouchi et al. 2009; Oesch et al. 2010a; McLure 
et al. 2010; Castellano et al. 2010a,b; Capak et al. 2011; 
Wilkins et al. 2011). Stepwise determinations of the LF 
from these studies are presented in Figure [TUl The cor- 
responding Schechter parameters are given in Table |6l 

It is instructive to compare the present LF determina- 
tion at z ~ 7 with those previously available. We begin 
with those studies where the primary focus was on the 
search for luminous, rarer z ~ 7 galaxies (e.g., Ouchi et 
al. 2009; Castellano et al. 2010a,b; Wilkins et al. 2010; 
Bouwens et al. 2010c). Perhaps the simplest point of 
comparison is the volume density of luminous (~ —21 
AB mag) galaxies at z ^ 7 relative to that at z ~ 6. 
Ouchi et al. (2009), Bouwens et al. (2010c), and Wilkins 
et al. (2010) find a factor of 2 decrease in the volume 
density of luminous (—21 AB mag) galaxies over this in- 
terval while Castellano et al. (2010a, b) find a factor of 
~3.5 decreased In th^ present search, we find a similar 
factor of ^1.8 decrease at the bright end of the LF (i.e., 
-20.5 AB mag). 

At fainter magnitudes, our HUDF09 search fields are 
essentially unique in the constraints they provide on the 
UV LF at Z'~7. Nonetheless, it is useful to compare 
against previous work, particularly from the early ultra- 
deep WFC3/IR pointing. Our results at —19 AB mag 
are in good agreement with the UV LF of Oesch et al. 
(2010a), but are somewhat higher than those found by 
McLure et al. (2010). 

Most studies find best-fit values of M* and at z 
7 that are slightly fainter and of lower density than is 
found at z-4-6 (Bouwens et al. 2008; Oesch et al. 2009; 
Ouchi et al. 2009; Oesch et al. 2010a; McLure et al. 2010; 
Castellano et al. 2010a, b). In the present study, we find 
something similar. The M* we find at z ~ 7 is fainter 
than the z '-.^ 6 value (Bouwens et al. 2007) by 0.1 
mag and the (jf we find at z ~ 7 is lower by a factor of 

Note that while Wilkins et al. (2010) frame the comparison 
between the z ~ 6 and 2 ~ 7 LFs in terms of the 2 ~ 6 Bouwens 
et al. (2006) LF, here we use the more up-to-date 2 ~ 6 LF results 
of Bouwens et al. (2007). 




Fig. 12. — Rest-frame UV LFs derived for galaxies at z ~ 7 {ma- 
genta circles with la errors) at 2 ~ 8 {open red circles), compared 
against similar LF determinations at 2 ~ 4 {blue), 2 ~ 5 {green), 
and z f> {cyan) from Bouwens et al. (2007). The 2 ~ 7 LF re- 
sults incorporate the Bouwens et al. (2010c) NICMOS + ISAAC + 
MOIRCS search results (see Figure [lOt. The upper limits arc Icr. 
The magenta and red lines show the best-fit Schechter functions 
at 2 ~ 7 and 2 ~ 8. The uniformly steep faint-end slopes a of the 
UV LF are quite apparent. Most of the evolution in the UV LF 
from 2 ~ 8 to 2 ~ 4 appears to be in the characteristic luminosity 
(by ~1 mag). 

^1.8. Of course, (/>* and M* are sufficiently degenerate 
at z 7 that there is a modest trade-off between the 
evolution found in either parameter. 

The faint-end slope to our UV LF is steep with a — 
—2.01 ± 0.21. Previously, it was not possible to set tight 
constraints on the faint-end slope. Ouchi et al. (2009) 
estimated a slope a of —1.72 ± 0.65 while Oesch et al. 
(2010a) - fixing (f)* to 0.0014 Mpc'^ - found -1.77±0.20. 
The present values of the faint-end slope are consistent 
with those measured at z^4-6 where this slope was also 
found to be steep at a —1.7 (Bouwens et al. 2007; 
Oesch et al. 2007; Yoshida et al. 2006; Yan & Windhorst 
2004). This new result may suggest a possible steepening 
towards higher redshifts z > 7, but the significance is 
weak, i.e., ~1.5(j. Such a steepening is expected in many 
theoretical models (e.g., Salvaterra et al. 2010; Trenti 
et al. 2010; Jaacks et al. 2011). Clearly we will require 
deeper data to confirm this trend. 

Capak et al. (2011) have suggested that the z ^ 7 
LF may extend to substantially brighter magnitudes 
than have been considered by other groups. Perform- 
ing follow-up spectroscopy on three very bright {H ~ 23 
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TABLE 7 

Determinations of the best-fit Schechter Parameters for. the 
rest-frame uv lfs at 2 ~ 8 (§6.1.2: see also figure [ttil. 



Reference 




(f>* (10-^ Mpc-3) 


a 


This work^ 


-20.10 ±0.52 


n 59+101 

'^■^^-0.37 


-1.91 ±0.32 


McLure et al. (2010) 


-20.04 


0.35 


-1.71 (fixed) 


Bouwens et al. {2010b) 


-19.5 ±0.3 


1.1 (fixed) 


-1.74 (fixed) 


Lorenzoni et al. (2011) 


-19.5 


0.93 


-1.7 (fixed) 



^ The derived Schechter parameters depend significantly on whether we include 
or exclude the two brightest (^^26 AB mag) 2 ~ 8 galaxies in the HUDF09-2 
field. If we exclude the two brightest galaxies (assuming they represent a rare 
overdensity ) , we derive a significantly fainter characteristic luminosity M* and 



larger normalization 



M* 



-19.54 ± 0.56, 



1.5^ 



X 10"'^ Mpc~ 



a = —1.67 ± 0.40 (similar to the values found by Bouwens et al. 2010b and 
Lorenzoni et al. 2011). 
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Fig. 13. — 68% and 95% likelihood contours on the model Schechter parameters for our UV (rest-frame ~1700A) LF determination at 
2 ~ 7 {magenta lines) and 2 ~ 8 {dotted red lines). For comparison, we also include the LF determinations at 2 ~ 4 (blue lines), z 5 
{green lines), and z ^ Q [cyan lines) from Bouwens et al. (2007). We do not include the 2 ~ 8 contours in the center and right panels 
given the large uncertainties on the 2 ~ 8 Schechter parameters. It is striking how uniform the rate of evolution in the UV LF (left and 
middle panels) is as a function of redshift - though there is a significant degree of degeneracy between <p* and M* . Most of the evolution 
in the LF appears to be in M* (particularly from 2 ~ 7 to 2 ~ 4). Within the current uncertainties, there is no evidence for evolution in 
(f)* or a {rightmost panel). 



mag) z ^ 7 candidates over the COSMOS field, Capak 
et al. (2011) claim to have uncovered statistically signif- 
icant detections of line flux. These lines could represent 
possible Lya emission given the existence of strong con- 
tinuum breaks in the sources. If this is the explanation, 
these sources would imply that the UV LF at z ~ 7 
extends to very bright magnitudes indeed, i.e., —24 AB 
mag, and in fact the LF that such sources would suggest 
is essentially power-law in shape. 

What should we make of this intriguing possibility? 
First, such a LF would seem to violate the constraints on 
the bright end of the LF set by deep, wide-area searches 
for 2 ~ 7 galaxies (Ouchi et al. 2009; Bouwens et al. 
2010c; Castellano et al. 2010a, b). Secondly, it is difficult 
to understand how one might create such bright —24 AB 
mag sources in the UV. The challenge is that any source 
forming stars rapidly enough to be so bright in the UV 
(i.e., ^200 Mo/yr) would most likely have built up suffi- 
cient quantities of dust and metals to be quite optically 
thick in the UV. Under such circumstances, it is unlikely 
that the source would be brighter than ~— 22 AB mag. 
Such is the situation at significantly later times (z~2-3) 
where there are exceedingly few system with UV mag- 
nitudes brighter than ~ —23 (e.g., Reddy et al. 2008), 



but many ULIRGs (e.g.. Chapman et al. 2005). It is 
hard to imagine therefore that aX z ^ 7 , just 500-600 
Myr after the first stars, that there is a population of 
galaxies that are substantially brighter in the UV than 
at z 3. Additional work utilizing deeper spectroscopic 
and photometric data is clearly required. 



.1.2. z 



LF 



There are also several published determinations of the 
UV LF at z ~ 8. These LFs are presented in stepwise 
form in Figure 1111 Schechter parameters for these LFs 
are given in Table [T] In general, the present z ~ 8 LF 
determination is in good agreement with previous deter- 
minations, within the quoted errors - though the present 
~ 8 LF is a little higher than the McLure et al. (2010) 
result at the lowest luminosities. 

For the derived Schechter parameters (see Table [7]), 
we find an M* of -20.10 ± 0.52, a 0* of 0.591^;!^^ x lO'^ 
Mpc~3, and faint-end slope a = -1.91 ± 0.32. These 
M* and 0* values are somewhat brighter and lower, re- 
spectively, than what we find at z^Q and z~7. This 
might indicative that some of the evolution in the UV 
LF from 2 ~ 8 to z < 6 can be represented by a change 
in the volume density 0*. However, the uncertainties on 



18 



(j)* and M* are too large to be sure of these conclusions. 
In particular, the existence of a possible overdensity at 
z - 8 within the HUDF09-2 field (see §5.3) illustrates the 
challenges of establishing reliable parameters from these 
initial datasets. 

6.2. Large-Scale Structure Uncertainties 

One of the most important challenges for LF determi- 
nations is large-scale structure. Variations in the volume 
density of galaxies as a function of position can have a 
significant impact on both the shape and overall normal- 
ization of the derived LPs. As a result of our use of the 
stepwise maximum likelihood and STY procedures (Ef- 
stathiou et al. 1988; Sandage et al. 1979) to derive the 
LF, the impact of large scale structure on the shape of 
our LF should be minimal. On the other hand, we would 
expect this structure to have an effect on the overall nor- 
malization of our LFs. One can estimate the size of these 
uncertainties using the Trenti & Stiavelli (2008) cosmic 
variance calculator. Using the approximate volume den- 
sities ~5xl0~^ Mpc"'^ we derived for z~7-8 sources, the 
width of our redshift selection windows Az^l, and the 
lO'xlO' search geometry we use within the CDF South, 
we find that the uncertainty in the normalization of the 
LF at z ^ 1 and z ~ 8 is a very modest ^^25% and ^^20%, 
respectively. 

Large-scale structure variations can also have an effect 
on the faint-end slope a determinations as a result of 
mass-dependent (and hence likely luminosity-dependent) 
bias parameters (e.g., Robertson 2010). We estimate the 
luminosity-dependent bias parameters based upon the 
approximate halo masses that correspond to the volume 
density of sources we are probing at specific points in the 
LF (again using the Trenti & Stiavelh 2008 calculator). 
We then perturb the LFs using these bias parameters 
and the power spectrum at z > 7 and determine the ef- 
fect on the faint-end slope a. We find that large-scale 
structure introduces a la uncertainty of ^0.05 in our 
faint-end slope determinations a at z ^ 7 and z ^ 8. 

We can also look at the issue of field-to-field variations 
from a more empirical standpoint. Fixing the shape of 
the UV LFs at z ~ 7 and z ~ 8 to the best-fit a and 
M*, we can look at the variations in the normalization 
(f)* for our three HUDF09 fields. We find a Icr variation 
in (j)* of just ~20%. This is similar to what we would 
expect from the Poissonian uncertainties. This suggests 
that large-scale structure uncertainties do not pose an 
especially huge problem for LF determinations, in the 
current luminosity range (—21 AB to —18 AB mag). 

Thus, the impact of large scale structure on our z '--^ 7 
and z 8 LFs appears to be quite modest. 

7. DISCUSSION 

7.1. Evolution of the UV LF from z ^ 8 to z ^ A 

In §5, we derived detailed constraints on the UV LFs 
at z ^ 7 and z ~ 8. By comparing these LFs with those 
derived at Z'^4-6, we can examine the evolution with 
cosmic time. In Figure [T^ we compare the stepwise LFs 
derived at z ^ 7 and z 8 with those found at z^A- 
6. The z ~ 7 LF presented in Figure [T^] incorporates 
the wide-area NICMOS + ISAAC + MOIRCS z - 7 
search results from Bouwens et al. (2010c: see the LFs in 
Figure llOp . It is remarkable how uniform the evolution 
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Fig. 14. — The derived faint-end slope a to the UV LF versus 
redshift. The present constraints on a are shown as the large solid 
red circles, with the Icr uncertainties plotted as the error bars. The 
values of a derived at z~4-6 (Bouwens et al. 2007: but see also 
Yoshida et al. 2006; McLure et al. 2009) are plotted as the open 
red circles. Other determinations are from Oesch et al. (2010a: blue 
circles; see also Hathi et al. 2010), Reddy & Steidel (2009: green 
squares), Arnouts et al. (2005: blue crosses), Wyder et al. (2005: 
blue squares), Treyer et al. (1998: black open circle), and Oesch et 
al. (2007, 2010a: green crosses). While the current constraints on 
the faint-end slope « at 2 ~ 7 and 2 ~ 8 are of limited statistical 
weight (see also Ouchi et al. 2009; Oesch et al. 2010a), they suggest 
the faint-end slope a is at least as steep as —1.7 - and possibly as 
steep as ~ —2.0 (§5.2; §5.3). We remark that the derived slopes 
simply describe the shape of the LF to several magnitudes below 
L* . Thus slopes that are even steeper than —2 do not imply a 
divergent luminosity density (since the LF will cut off at some 
faint value). 

of the UV LF is with cosmic time. At higher luminosities 
(-20 AB mag), the UV LF evolves by a factor of --6-7 
from z 8 to z ~ 4 while at lower luminosities (—18 AB 
mag), the decrease is a factor of ^^5. 

Figure [T3l presents the 68% and 95% confidence inter- 
vals derived here for the UV LF at z ~ 7 and z ^ 8. For 
comparison, the LFs at z ~ 4, z ~ 5, and z ~ 6 are also 
shown so that the evolution can readily be inferred. A 
clear trend is evident in the M* /4>* LF fit results {left 
panel) from z ^ 8 to z 4. The characteristic lumi- 
nosity shows evidence for a consistent brightening with 
cosmic time from z ~ 7 to z ~ 4. Ascertaining how 
the 0* and M* evolves out to z 8 is not possible at 
present, given the large uncertainties and degeneracy be- 
tween these two parameters (though there is a hint that 
some of the evolution may be in (jf : see also discussion in 
McLure et al. 2010). Tighter constraints on the volume 
density of bright z 8 galaxies are required for progress 
on this question. 

As far as constraints on the shape of the LF [middle 
panel), again we see a clear trend towards fainter val- 
ues of M* from z '-^ 7 to z ~ 4. The faint-end slope 
a at z 7 and z ~ 8 also is very steep, similar to the 
values a ^ —1.7 found at z^2-6 (Figure [T31) and pos- 
sibly even steeper. This possible steepening of the LF 
with redshift is in qualitative agreement with that pre- 
dicted from various theoretical models for the LF (e.g., 
Trenti et al. 2010; Salvaterra et al. 2010; Bouwens et al. 
2008). Such steep faint-end slopes are potentially very 
important for reionization since it would suggest that low 
luminosity galaxies at z > 7 are major contributors to 
the needed UV ionizing photon flux. We discuss this in 
detail in Bouwens et al. (2011b). 
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Fig. 15. — 68% and 95% likelihood contours on the model 
Schechter parameters M* and (p* for our LF determination at z ~ 7 
{magenta lines) and 2 ~ 8 {dotted red lines). For comparison, we 
also show the constraints on the UV LF at 2 ~ 4 {blue lines), 
2 ~ 5 {green lines), and 2 ~ 6 {cyan lines) from Bouwe ns e t al. 
(2007). This figure is the same as the left panel to Figure [131 but 
now fixing the faint-end slope a to —1.7. A faint-end slope a of 
— 1.7 is similar to the values found at 2~l-6 (Yoshida et al. 2006; 
Bouwens et al. 2007; Oesch et al. 2007; McLure et al. 2009; Reddy 
& Steidel 2009; Oesch et al. 2010c; Arnouts et al. 2005: see Fig- 
ure ll4|l . No constraint on the faint-end slope at 2 ~ 4-6 is imposed, 
and therefore the error contours there are somewhat larger than at 
2 ~ 7. 

No trend is evident as a function of redshift in (jj* and 
a considered together {right panel: Figure I13p . How- 
ever, the uncertainties remain large at ^ ~ 7 and z 8. 
Given these uncertainties, we can potentially gain more 
insight into the evolution of the LF at high redshift by 
fixing the faint-end slope a to —1.7. The faint-end slope 
a is found to be consistent with this value over a very 
wide redshift range z~l-8 (Bouwens et al. 2007; Reddy 
& Steidel 2009; Beckwith et al. 2006; McLure et al. 2009; 
Ouchi et al. 2009; Oesch et al. 2007, 2010a, 2010c), so it 
seems reasonable to fix a to —1.7 at z 7 and z ^ 8. 
The constraints we derive on M* and (f>* fixing a are 
shown in Figure [T51 This figure provides further support 
to the suggestion that the primary mode of evolution in 
the LF from z ~ 7 to z ~ 4 is in the characteristic lumi- 
nosity (see also Bouwens et al. 2006; Yoshida et al. 2006; 
Bouwens et al. 2007; McLure et al. 2009). For z ~ 8 LF, 
the normalization is somewhat lower - suggesting that 
some of the evolution in the LF at z > 4 may be in vol- 
ume density (e.g., McLure et al. 2009; Castellano et al. 
2010b). 

7.2. Shape of the UV LF at z^7 and z^8 

One potentially important open question about the 
UV LFs at z > 7 concerns their overall shape (e.g., see 
discussion in §5.5 of Bouwens et al. 2008). While we 
have generally attempted to represent the LFs here with 
Schechter functions (i.e., a power-law with an exponen- 
tial cut-off at the bright end: see §5.2 and §5.3), it is 
possible that the LFs may have a shape closer to that 



TABLE 8 

UV Luminosity Densities and Star Formation Rate Densities 

TO -17.7 AB MAG.'' 



Dropout 
Sample 


< 2 > 


logio^ 
(ergs s^^ 
Hz-i Mpc-3) 


logjQ SFR density 
{Mq Mpc-3 yr-l) 
Uncorrected Corrected'' 


2 

Y 

Jd 
Jd 


6.8 
8.0 
10.3 
10.3 


25.88±0.10 
25.65±0.11 

24 OQ"^"'^-'" 
<24.42o 


-2.02 ±0.10 
-2.25 ±0.11 

< -3.48<: 


-2.02 ± 0.10 
-2.25 ± 0.11 
-3 61+°-5' 
< -3.48^= 


B 
V 

i 


3.8 
5.0 
5.9 


26.38±0.05 
26.08±0.06 
26.02±0.08 


-1.52±0.05 
-1.82±0.06 
-1.88±0.08 


-0.90 ±0.05 
-1.57 ± 0.06 
-1.73 ±0.08 



Integrated down to 0.05 i*— 3- Based upon LF parameters in 
Table [5] (see §7.3). The SFR density estimates assume > 100 Myr 
constant SFR and a Salpeter IMF (e.g., Madau et al. 1998). Conver- 
sion to a Kroupa (2001) IMF would result in a factor of ~1.7 (0.23 
dex) decrease in the SFR density estimates given here. 
^ The dust correction is taken to be for galaxies at 2'~7-8, given 
the very blue /3's (see also Bouwens et al. 2009; Oesch et al. 2010a; 
Bouwens et al. 2010a). 

Upper limits here are Ict (68% confidence). 

2 ~ 10 determinations and limits are from Bouwens et al. (2011a) 
and assume 0.8 2 ~ 10 candidates in the first case and no 2 ~ 10 
candidates (i.e., an upper limit) in the second case. 

of a power-law, e.g., resembling the halo mass function 
(e.g., Douglas et al. 2009). This is simply because at 
the high redshifts we are probing the typical halo masses 
would be much lower (<10^^ Afo) and therefore many of 
the physical mechanisms likely to be important for trun- 
cating the LFs at z (e.g., transition to hot flows: 
Birnboim & Dekel 2003; Rees & Ostriker 1977; Binney 
1977; Silk 1977 or AGN feedback: Croton et al. 2006) 
may not be as relevant at z > 8. Of course, at suffi- 
ciently high redshifts (i.e., z > 7), we would eventually 
expect some truncation in the LF at high masses (and 
luminosities) due to a similar exponential cut-off in the 
halo mass function (e.g., Trenti et al. 2010). 

Overall, we do not find any evidence that the UV LF at 
z ~ 7 or z ~ 8 has a fundamentally different shape than 
the LF has at lower redshift. Both LFs seem to have 
roughly a power-law-like shape at the faint end and a 
modest cut-off at the bright end - consistent with that of 
a Schechter function (though the evidence for such a cut- 
off in the z ~ 8 LF is weak). Ascertaining whether the LF 
really has a Schechter-like form or has a tail extending to 
bright magnitudes (e.g., Capak et al. 2011) will require 
significantly more statistics, particularly at bright magni- 
tudes, to better establish the shape of the LF there. For- 
tunately, we should soon have such statistics as a result 
of observations being obtained over the GOODS fields as 
part of the CANDELS multi-cycle treasury program, the 
wide-area SEDS/CANDELS fields, and two pure parallel 
programs (Trenti et al. 2011; Yan et al. 2011). 

7.3. Luminosity Densities and SFR Densities 

With the completion of the wide-area ERS and ultra- 
deep HUDF09 observations, we have been able to signif- 
icantly improve our constraints on the UV LFs at z ^ 7 
and z ^ 8. These new LF determinations allow us to up- 
date our luminosity density and SFR density estimates at 
Z'^7-8 from that already given in Bouwens et al. (2010b) 
based on the initial HUDF09 data on the HUDF. These 
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Fig. 16. — Luminosity density and SFR density determinations, as a function of redshift (§7.3). These determinations are integrated 
to —17.7 AB mag (0.05 L*_^). The lower set of points (blue for our measurements) and shaded blue region are the observed luminosity 
densities and SFR density estimates before dust correction. The upper set of points (red for our measurements) and shaded orange region 
indicates the SFR density estimates after dust correction. The SFR density estimates assume > 100 Myr constant SFR and a Salpeter 
IMF. Conversion to a Kroupa (2001) IMF would result in a factor of ~1.7 (0.23 dex) decrease in the SFR density estimates given here. 
At z < 2 dust corrections are from Schminovich et al. (2005), at z~3-6 from Bouwens et al. (2009b), and at z > 7 they are assumed to be 
zero. The solid circles at z ~ 7 and z ~ 8 are the present determinations. For context, we have also included the SFR density estimates at 
z~4-6 from Bouwens et al. (2007), at z~2-3 from Reddy & Steidel (2009: green crosses), and at z<2 from Schiminovich et al. (2005: black 
squares). The estimates at z ~ 10 (Bouwens et al. 2011a) are given with the blue circles and upper limit. See Table[8]for a tabulation of 
the z > 4 luminosity densities and SFR densities presented here. 



densities continue to be of great interest for discussions 
about the role of galaxies in reionization as well as the 
build-up of stellar mass. 

For our luminosity density (and SFR density) esti- 
mates, we simply integrate up our stepwise LF determi- 
nations to the limiting depth of our LF determinations, 
i.e., to -17.7 AB mag (0.05 i*^3). The resuhant lu- 
minosity densities are tabulated in Table [S] They are 
also shown in Figure [16] along with luminosity density 
determinations at lower redshift. 

The luminosity density determinations are converted 
into SFR densities using the Madau et al. (1998) con- 
version factor - assuming a Salpeter IMF (with a stellar 
mass function ranging from 0.1 Mq and 125 Mq). As 
often noted (e.g., Verma et al. 2007; Stanway et al. 2005; 
Bouwens et al. 2009), this conversion factor assumes that 
the SFRs in galaxies are relatively constant for > 100 
Myr prior to observation. Consequently, use of this con- 
version would result in an underestimate of the SFR (by 
factors of ~2) for galaxies with particularly young ages. 

Given the very blue [/V^- continuum slopes /3 of z > 7 
galaxies (e.g., Bouwens et al. 2010a, b; Oeschet al. 2010a; 
Finkelstein et al. 2010; Bunker et al. 2010; see also 
Bouwens et al. 2009), the dust extinction in z > 7 galax- 
ies is likely small, and so no correction for dust is made in 
computing the SFR densities at z > 7. At z<6, we adopt 



the dust corrections given in Schiminovich et al. (2005), 
Reddy & Steidel (2009), and Bouwens et al. (2009). Our 
SFR density estimates are included in Table [8] and on 
Figure [161 

The UV luminosity densities we find to —17.7 AB mag 
at z ~ 7 and z ~ 8 are just 32% and 18%, respectively, 
of that found to the same limiting luminosity at z ^ 4. 
The SFR densities we find to these limits are just 6% and 
4%, respectively, of that at z ^ 4. This provides some 
measure of how substantial galaxy build-up has been in 
the '^l Gyr from z~7-8 to z ~ 4 when the universe was 
nearing its peak SFR density at z 2-3. 

7.4. SFR Densities from the Stellar Mass Density 

It is interesting to compare the present SFR density 
determinations with that inferred from recent estimates 
of stellar mass density at z > 4 (Labbe et al. 2010a, b; 
Stark et al. 2009; Gonzalez et al. 2010). The comparison 
provides us with a very approximate test of the shape 
of the IMF at high redshift as well as allowing us to test 
the basic consistency of these very fundamental measure- 
ments. 

In estimating the SFR densities from the stellar mass 
densities (Stark et al. 2009; Labbe et al. 2010a, 2010b; 
Gonzalez et al. 2010), there are a few issues we need 
to consider. One issue is the stellar mass lost due to 
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Fig. 17. — Comparison of the SFR density determinations at 
z > A with that implied by recent determinations of the stellar 
mass density (§7.4). The red hatched region gives the SFR densi- 
ties implied by current compilations of the stellar mass density at 
2 > 4 (Labbe et al. 2010a,b; Stark et al. 2009; Gonzalez et al. 2010) 
where the width of this region indicates a roughly ztlcr uncertainty. 
The change in density of the cross-hatching at 2 > 7 is intended 
to indicate that the mass-density-derived SFR may be more sub- 
ject to systematic uncertainty from limited measurements and the 
extrapolation required from the measurements at z < 7. At 2 > 8, 
the stellar mass density in galaxies brightward of —18 AB mag is 
assumed to scale with redshift as (1 + z)"**-^ (as found in Labbe et 
al. 2010b). Interestingly, at 2 < 7 the mass-density-derived SFR 
is only about 40% larger than the SFR density estimated from the 
LFs; such a small difference is quite encouraging. This suggests 
that (1) the IMF at high redshift {z ~ 4-7) may not be that differ- 
ent from what it is at later times and (2) that stellar mass density 
estimates at z > 4 are reasonably accurate. 



SNe recycling (Bruzual & Chariot 2003). To correct for 
this, we have multiplied the star formation rate densities 
we infer from the stellar the stellar mass densities by a 
factor of 1/0.72 1.39 (appropriate for the Salpeter IMF 
assumed in the Labbe et al. 2010b mass estimates and a 
~300 Myr typical age). 

A second issue we need to consider is that new galax- 
ies will be entering magnitude- limited samples at all red- 
shifts simply due to the growth of galaxies with cosmic 
time (e.g., galaxies not present in a z ~ 7 sample would 
grow enough from 2 7 to z ^ 6 to enter magnitude- 
limited selections at z ~ 6). Therefore, it is not quite 
possible to determine the SFR density by subtracting 
the stellar mass density determination in one redshift in- 
terval from that immediately below it. To account for 
those galaxies entering magnitude-limited samples due 
to luminosity growth, galaxies are assumed to brighten 
by 0.33 mag per unit redshift (the approximate M* evo- 
lution: see §7.5). This results in a typical 20% correction 
to the SFR density. Based on these corrections and con- 
siderations and the stellar mass density at z ~ 8 (Labbe 
et al. 2010b), we derive the SFR density down to z ~ 4. 
The results are shown in Figure [IT] ( /latc/ied red region). 
Fl Interestingly, the SFR densities inferred from the stel- 
lar mass density are only ~40% higher than that inferred 
from the UV LFs at z < 7 (where the stellar masses are 

The stellar mass density estimated at z ~ 8 (Labbe et al. 
2010b) - although uncertain - suggests a certain SFR density at 
2 > 8. Of course, since the stellar mass is an integral quantity 
(and there are a large range of plausible SF histories that yield the 
observed stellar mass density), a redshift dependence needs to be 
assumed. We adopt the same scaling of the stellar mass density 
with redshift at 2 > 8, i.e.,(l -|- 2)-*-4, as Labbe et al. (2010b) 
found to fit the observational results at 2~4-8. 



more well determined) E3 This suggests that (1) the IMF 
for SF galaxies at high redshift (z ~ 4-7) may not be that 
different from what it is at later times (z < 3: but see dis- 
cussion in Chary 2008) and (2) the stellar mass densities 
at z > 4 inferred from the Spitzer IRAC photometry are 
reasonably accurate. This latter conclusion is important 
since nebular emission line fluxes have been suggested 
to dominate the measured light in z > 7 galaxies with 
IRAC (Schaerer & de Barros 2010: but see Finlator et al. 
2011). These comparisons suggest that nebular emission 
is not as dominant as has been suggested (Schaerer & de 
Barros 2010). 

7.5. Extrapolations to z > 8: A LF-Fitting Formula 

The present LF determinations provide us with reason- 
ably accurate measures of luminosity density and SFR 
density at z ^ 7 and z 8 - which are of considerable 
importance for assessing the impact of galaxies in the 
build-up of mass and metals in the universe, as well as 
understanding their role in reionization. However, know- 
ing the LF results at z^7-8 only provide us with part of 
the information we want. To address some of the most 
outstanding of today's questions (i.e., whether galaxies 
reionize the universe, how stellar mass and metals are 
built up versus cosmic time), we really need to know 
what the UV LFs look like at even higher redshift. 

For these questions, perhaps the best approach is sim- 
ply to extrapolate the present LF results to even higher 
redshifts. Certainly, given the very uniform rate of evo- 
lution of the UV LF from z ~ 8 to z '--^ 4, it seems quite 
plausible that such an extrapolation would be reason- 
ably accurate. As in our previous work (§5.3 of Bouwens 
et al. 2008), we frame the evolution of the LF in terms 
of a Schechter function and assume that each Schechter 
parameter M* , a, and (j)* can be expressed as a linear 
function of redshift. We then identify linear coefficients 
that maximize the likelihood of reproducing the observed 
contours in Figure [T51 at z '-^ 4, z ~ 5, z ~ 6, z ~ 7, and 
z ^ 8. The result is 

M^y = (-21.02 ± 0.09) + (0.33 ± 0.06)(z - 3.8) 

(j)* = (1.14 ± 0.20)10(" ™3±0.055)(^-3.8)io-3Mpc-3 

a = (-1.73 ± 0.05) + (-0.01 ± 0.04)(z - 3.8) 

We see that the evolution in M^y is significant at the 
5(7 level. However, the evolution in (jf and a is still 
not significant. The fitting formula above are similar to 
those found in Bouwens et al. (2008) based upon a large 
z - 4-6 ACS sample and z - 7 NICMOS sample. 

7.6. Implications for z ^ 10 Ji2^-dropout Searches 

In the previous section §7.5, we combined the present 
constraints on the UV LFs at z 7 and z 8 with that 
available at z^4-6 to derive fitting formula for the evolu- 
tion of the UV LF with redshift. These formula arguably 
provide us with our best means for extrapolating the UV 
LF to even higher redshifts and therefore estimating the 
form of the UV LF at z 10. This enterprise has impor- 
tant implications for z 10 galaxy searches in present 
and future WFC3/IR observations. 

The difference increases slightly at earlier times where the 
galaxy stellar mass function is more poorly determined and ex- 
trapolations are used. 
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What surface density of z ^ 10 Ji25-dropout can- 
didates do we expect to find in current ultra-deep 
WFC3/IR observations? We can estimate this by taking 
advantage of the fitting formula we derive in the previ- 
ous section and the detailed simulations carried out by 
Bouwens et al. (2011a) to estimate Ji25-dropout selec- 
tion efficiencies. We predict 1.0 and 2.3 Ji25-dropouts 
arcmin"^ to an Hiqq magnitude of 28.5 and 29.5 AB mag, 
respectively. These numbers are similar albeit somewhat 
larger than that found by Bouwens et al. (2011a) over 
the HUDF to Hieo ^ 29.5. Certainly, the fact that 
a few Ji25-dropout galaxies are expected in the HUDF 
search suggests that the z 10 galaxy candidate iden- 
tified by Bouwens et al. (2011a) is plausible. Of course, 
deeper observations are clearly needed for more defini- 
tive statements. We do not address the Yan et al. (2010) 
'^125-dropout numbers here, given the challenges in in- 
terpreting those numbers (both as a result of very low 
3(7 detection thresholds used in that study and the likely 
high contamination rates expected from the proximity of 
the many Yan et al. 2010 candidates to bright foreground 
galaxies: see discussion in Bouwens et al. 2011a). 

8. SUMMARY 

The ultra-deep (14 arcmin^) HUDF09 and deep wide- 
area (40 arcmin^) ERS WFC3/IR observations are a sig- 
nificant resource for identifying large samples of star- 
forming galaxies at z > 7. The depth of these ob- 
servations in three distinct wavelength channels com- 
bined with their availability over regions of the sky with 
ultra-deep optical data make such identifications possible 
with well-established selection procedures like the Lyman 
Break Galaxy technique. 

We utilize LBG selection criteria over these fields to 
identify 73 candidate z ~ 7 galaxies and 59 z ^ 8 can- 
didate galaxies (Table [2]). These candidates extend over 
a ~3 mag magnitude range in luminosity from Hiqq^ab 
~26 to ^29.4 mag - thus reaching to —17.7 AB mag 
(0.05 ^2=3). The use of LBG criteria allows us to iden- 
tify likely high-redshift sources from the extremely deep 
WFC3/IR data and then to exploit fully the deep, high 
spatial resolution ACS optical data to remove contami- 
nants in a very robust way. Our sample of 132 galaxies in 
the reionization epoch is the largest currently available. 

To determine the contamination rate we thoroughly 
assess the impact of lower redshift sources, photomet- 
ric scatter, AGN, spurious sources, low mass stars, and 
transients (e.g., SNe) on our selection (§3.5). The con- 
tamination rate that we find is just ^7% for the HUDF 
sample and ^^30% for the brighter, but less crucial, ERS 
sample (where the shallower optical data provide weaker 
constraints). 

Through simulations, we have found that very strict 
optical non-detection criteria are essential to obtaining 
contamination-free z > 7 samples. We set a requirement 
that an accepted source must not show a 2cr-detection 
in any optical filter (nor a 1.5CT-detection in >1 optical 
filter). However, through extensive simulations we have 
found that this alone is not rigorous enough; it is likely 
that -30-40% of our LBG samples over the HUDF09-1, 
HUDF09-2, and ERS fields would have consisted of con- 
taminants if that had been the only optical constraint. 
What is much more effective is to combine the informa- 
tion in the entire optical dataset to derive a xtpt value 



for each source. In addition, we recommend validating 
the selection by looking at the xtpt distribution for the 
entire sample (§4.2: Figure [TS]) . The use of this tech- 
nique results in much smaller contamination and is the 
reason why we can reach as faint as 29 AB mag and still 
reach an overall contamination rate of —11% (—7% in the 
HUDF09 selection). This is discussed extensively in §3.5 
and Appendix D, while our concerns about contamina- 
tion in some other recent studies are discussed in §4.1.2 
and §4.1.3. 

The high surface density of galaxies in our samples has 
direct implications for the LFs of galaxies in the first 800 
Myr of the universe and the build-up of galaxies with 
cosmic time. We derive selection volume for our LBG- 
selected galaxy samples through careful simulations, and 
use these selection volumes to derive luminosity functions 
with well-understood errors. The techniques we use to 
derive the LFs are insensitive to large-scale structure bi- 
ases and also explicitly account for other observational 
biases. We combine our LF determinations with other 
wide-area search results (see Table S]) to derive more ro- 
bust constraints on the overall shape of the LF. We then 
compare these LFs with those found at z— 4-6 to quan- 
tify the evolution of the UV LFs in the first 1.6 Gyr of 
the universe. 

Here are our primary findings: 

• The faint-end slope a for the UV LF at z — 7 
is very steep, with values of —2.01 ± 0.21. The 
current faint-end slope determination is consistent 
with previous estimates at z — 7 (Ouchi et al. 2009; 
Oesch et al. 2010a) and is a logical extension of that 
found at lower redshift where a~— 1.7 from z — 6 
to z — 1 (Oesch et al. 2010c; Bouwens et al. 2007; 
Reddy & Steidel 2009; Oesch et al. 2007; Yoshida 
et al. 2006; Beckwith et al. 2006; Yan & Windhorst 
2004). The faint-end slope a at z ~ 8 is more un- 
certain, but also consistent with being very steep: 
a = — 1.91±0.32. We have verified that the present 
faint-end slope determinations are not significantly 
affected by large-scale structure (see §6.2). Such 
steep slopes do not necessarily lead to a divergent 
integrated luminosity density since it is expected 
that the galaxy LF will almost certainly truncate at 
some physically-established low-mass or luminosity 
limit. These results provide increasingly strong ev- 
idence that very low luminosity galaxies provide 
the dominant contribution to the luminosity den- 
sity at z > 6. Strikingly, >75% of the luminosity 
density is from galaxies fainter than —18 AB mag 
(if a < —1.7). Interestingly, such large contribu- 
tions are expected from theory (e.g., Salvaterra et 
al. 2010; Trenti et al. 2010; Barkana & Loeb 2000). 

• The characteristic luminosity M* and the normal- 
ization (j)* we find at z ~ 7 are —20.14 ± 0.26 and 
0.86t^:^^ X 10~3 Mpc-3, respectively. These values 
are consistent with other recent determinations of 
the UV LF at z ^ 7 (Ouchi et al. 2009; Oesch 
et al. 2010a; Castellano et al. 2010a) and are also 
consistent with lower redshift trends (i.e., Bouwens 
et al. 2007, 2008). If we represent the z - 8 LF 
with a Schechter function, the characteristic lumi- 
nosity M* and normalization 0* that we determine 
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is -20.10 ± 0.52 and OMtlil x lO"^ Mpc'^, re- 
spectively. This suggests that the evolution in the 
UV LF at z > 4 occurs not only through changes 
in M* (as preferred from the z ^ 7 to z ^ A results: 
§7.1), but also through changes in <j>* . The result 
however is quite sensitive to the surface density of 
very luminous galaxies [in this case the two bright 
2: ~ 8 galaxies in the HUDF09-2 field which would 
appear to represent an overdensity: §5.3]. The un- 
certainties are still large at present. 

The 2: ~ 7 LF appears to be much better repre- 
sented by a Schechter function than by a power 
law and shows evidence for an exponential cut-off 
at the bright end (§7.2). The z ^ 8 LF also seems 
to be better represented by a Schechter function - 
though the evidence is weaker. It is unclear what 
physical explanation would impart such a cut-off 
in the LF (particularly if M* becomes fainter at 
higher redshift: §7.5), but it is possible this may 
be due (at least partially) to a similar exponential 
cut-off in the halo mass function (e.g., Trenti et al. 
2010). 

By integrating the UV LFs derived here at Z'^7 and 
8, we estimate the UV luminosity density and SFR 
density at 2; ~ 7 and 2 ~ 8 to a limiting luminosity 
of -17.7 AB mag (0.05 L*=3). Consistent with 
the very blue UV-continuum slopes /3 observed for 
2~7-8 galaxies (e.g., Bouwens et al. 2010a, b; Oesch 
et al. 2010a; Finkelstein et al. 2010; Bunker et al. 
2010), we assume that the dust extinction is zero. 
The UV luminosity densities at Z'-~^7 and z^8 are 
~32% and ~18% of the value at z ~ 4, while the 
SFR densities estimated at z~7 and Z'-^8 to these 
same luminosity limits are ~6% and ^4% of the 
value at z - 4 (see Figure [16] and Table [D §7.3). 

The SFR densities we estimate from published stel- 
lar mass density determinations at z ~ 4-7 (e.g., 
Labbe et al. 2010b; Stark et al. 2009; Gonzalez et 
al. 2010) are only ~40% higher than that derived 
directly from the LFs. See §7.4 and Figure[l7l This 
suggests that (1) the initial mass function for star 
formation at very high redshift may not be that 
different from what it is at later times and (2) the 
stellar mass densities at z > 4 inferred from the 
IRAC photometry are reasonably accurate. The 
latter conclusion is important since it has been sug- 
gested that nebular emission line fluxes may dom- 
inate the measured light in 2; > 7 galaxies with 
Spitzer IRAC (Schaerer & de Barros 2010: but see 
also Finlator et al. 2011). The present comparison 
suggests otherwise. 

We use the current contraints on the LF at z~7- 
8 to update the Bouwens et al. (2008) fitting for- 
mula describing the evolution of M*, 0*, and a as 
a function of redshift (§7.5). The revised formula 
areMljy = (-21.02±0.09)-f (0.33±0.06)(z-3.8), 

4>* - (1.14 ± 0.20)10(0 ™3±0 055)(-3.8)io-3Mpc-3, 

and a ^ (-1.73 ± 0.05) + (-0.01 ± 0.04)(2 - 3.8) 
(§7.5). 



• Given the interest in the existence of massive black 
holes at z > 6 we have considered whether our z ^ 
7 and z ~ 8 samples include any obvious evidence 
in the WFC3/IR images for AGNs. It appears that 
the frequency of obvious AGN (unresolved "stellar" 
sources that are comparable to, or brighter than, 
the underlying galaxy) is quite small (i.e., ^2%) in 
these z ~ 7-8 samples (§3.5.3). 

• These new LF results allow us to estimate the ex- 
pected number of 2 ^ 10 galaxies from a search of 
the current IIUDF09 data. We predict 1.0 and 2.3 
>/i25-dropouts arcmin^^ to an iJigo-band magni- 
tude of 28.5 and 29.5 AB mag, respectively. These 
numbers are consistent with the single z ^ 10 can- 
didate found by Bouwens et al. (2011a) over the 
HUDF to ~ 29.5 (§7.6). See Oesch et al. (2011) 
for an updated discussion using a larger set of ob- 
servations. 

The present WFC3/IR observations allowed us to select 
132 galaxies at z^7 and 2^8, over a 3 mag range in 
luminosity, reaching as faint as 29 AB mag in the near- 
IR. These observations have provided us with our best 
insight yet into the build-up and evolution of galaxies at 
early times - providing us with a clear measurement of 
rate of evolution in the UV LF from 2 '--^ 8 to z ~ 4 in 
the first 1.6 Gyr of the universe. Surveying the results, 
it is remarkable how uniform the evolution of the LF is 
with cosmic time. These new observations also provide 
us with further insight to the contributions of galaxies to 
reionization (see Bouwens et al. 2011b for new estimates 
based on these deep LFs). 

Further progress should be made over the next few 
years in the field of high-redshift LFs from the wide- 
area HST Multi-Cycle Treasury (MCT) CANDELS pro- 
gram. The CANDELS program will substantially im- 
prove our constraints on the prevalence of bright galax- 
ies at 2 > 7 and the overall shape of the LF. However, 
the crucial observations needed for ultimately making 
progress on characterizing the role of galaxies in reioniz- 
ing the universe in the period from z^ll to z^6 and on 
the nature of the dominant population of low-luminosity 
galaxies at Z'~7-10 are further ultra-deep observations 
with WFC3/IR. Such observations would significantly 
improve our determinations of the faint-end slope a, the 
J/y-continuum slopes /? of very low luminosity galaxies 
(e.g., Bouwens et al. 2010a), and the prevalence of 2 ~ 10 
galaxies (e.g., Bouwens et al. 2011a). 
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Fig. 18. — Comparison of the half-light radii of 2 ^ 7-8 galaxies observed over the HUDF09 field {histogram) with those from our fiducial 
galaxy model (thick solid line: which arc used to estimate the selection volumes, see Appendix C). The left panel is for our 2 ~ 7 selection 
and the right panel is our for 2 ~ 8 selection. Galaxies are selected from our simulations in exactly the same way as from the observations. 
The magnitude intervals plotted are chosen to be sufficiently brightward of the faint limit to allow for the selection of galaxies over a range 
of sizes, while still including a significant fraction of our HUDF09 z ~ 7-8 sample. Also shown is the half-light radius distribution for a 
galaxy population with sizes that are 40% smaller and 50% larger than our fiducial model (dashed lines and dotted lines, respectively). 
The observed half-light radius distribution is in excellent agreement with the sizes from our fiducial model. The fit to the other size models 
is less good (the disagreement ranges from ~1-2(t), but the uncertainties are large. Deeper observations will be required to establish the 
size distribution more accurately to very low luminosities. We find no evidence that our model size distribution is too large, and therefore 
that our selection volumes are too low, as suggested by Grazian et al. (2011: see Appendix A). 
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APPENDIX 

A. SIZES OF OUR MODEL GALAXY POPULATION AND A COMPARISON WITH THE OBSERVATIONS 

An essential input to determinations of the luminosity functions are estimates of the selection volume. For deter- 
mining these selection volumes, one of the most important inputs is the size (or surface brightness) distribution of the 
model galaxy population (e.g., Bouwens et al. 2003a, b; Bouwens et al. 2004; Bouwens et al. 2006; Oesch et al. 2007). 
Use of a model galaxy population with sizes that are too large results in selection volumes that are too low, while use 
of a model galaxy population with sizes that are too small results in selection volumes that are too high. 

Fortunately, we do not expect this to be an issue because of our approach. We expect that the sizes (or surface 
brightnesses) of our model galaxy population will match that seen in the observations since we base the sizes of our 
model galaxy population on the sizes already measured for z ^ 7-8 galaxies in the HUDF (Oesch et al. 2010b) and on 
the size evolution observed from z ~ 8 to z ~ 4. Nonetheless, it is useful to verify that there is good agreement between 
our model galaxy population and observations using our current z ^ 7-8 samples and the full two-year HUDF09 data. 

We provide such a check in Figure [TS] - comparing the sizes of z 7 galaxies seen in the observations and from our 
simulations. The simulations are based on z ~ 4 galaxies from the HUDF and the size-redshift evolution seen from 
z ~ 7 to z 4. Sources from the simulations are selected in an identical manner. For this comparison, we only utilize 
galaxies between Ji25,ab ~ 27.8 and Ji25,ab ^ 28.5 for our z ^ 7 selection and i?i6o,AB 27.8 and Hiqo ab ^ 28.6 
for our z ~ 8 selection - since these magnitude intervals allow for the selections of galaxies over a range of sizes and 
thus usefully probe the size distribution. Faintward of this, only the smallest galaxies are selected, and so comparisons 
between the observed and model size distributions do not allow us to effectively test our model (see discussion in §3.7 
and Figure 5 of Bouwens et al. 2006; Appendix A of Bouwens et al. 2008; Figure 3 of Grazian et al. 2011). Also shown 
in this figure is the half-light radius distribution for a galaxy population with sizes that are 40% smaller and 50% 



26 



TABLE 9 

Estimated Number of Contaminants in our ^ ~ 7 and z --^ 8 samples. 



Mag Range'' HUDF'^ HUDF09-1 HUDF09-2d ERS 



z 1^ 1 Samples 



25.0-25.5 


0.00 


0.00 


0.00 


0.01 


25.5-26.0 


0.00 


0.00 


0.00 


0.24 


26.0-26.5 


0.00 


0.02 


0.03 


0.66 


26.5-27.0 


0.00 


0.09 


0.06 


1.05 


27.0-27.5 


0.02 


0.15 


0.15 


0.85 


27.5-28.0 


0.13 


0.41 


0.33 


0.12 


28.0-28.5 


0.20 


0.86 


0.63 




28.5-29.0 


0.36 


0.58 


0.71 




29.0-29.5 


0.27 


2 ~ 8 Samples 






25.0-25.5 


0.00 


0.00 


0.00 


0.01 


25.5-26.0 


0.00 


0.00 


0.00 


0.19 


26.0-26.5 


0.00 


0.00 


0.00 


0.47 


26.5-27.0 


0.00 


0.00 


0.00 


0.70 


27.0-27.5 


0.00 


0.06 


0.02 


0.86 


27.5-28.0 


0.02 


0.25 


0.15 




28.0-28.5 


0.17 


0.49 


0.39 




28.5-29.0 


0.50 


0.45 


0.38 




29.0-29.5 


0.41 


0.01 


0.06 





For our LBG selections over our three HUDF09 fields, our estimated contamination rates 
are made assuming that the faint photometric samples have the same color distributions as 
found in the high S/N HUDF09 26.5 < -ffieo.AS < 28.0 subsample. For our ERS selections, 
our contamination estimates are based upon the high S/N HUDF09 26 < i^iso AB < 27.0 
subsample. Noise is then added to the photometry, and our 2: ~ 7 or z ~ 8 selection criteria 
are then applied. Any source that is selected by our LBG criteria (after adding noise), but 
detected in the optical in the original high S/N observations is counted towards the total 
contamination level. See §3.5 for details. 

^ J125 and Hiao magnitudes used for our 2^7 and 2 ~ 8 selections, respectively. 

It is striking how much lower the estimated contamination levels are in our HUDF selections 
than in our HUDF09-1/ HUDF09-2 selections. This difference is a direct result of the much 
greater depth of the optical data over the HUDF (by ~1 mag) than the HUDF09-1/HUDF09- 
2 selections. 

As a result of the ACS parallel observations acquired as part of the HUDF09 WFC3/IR 
program, the effective contamination rate over the HUDF09-2 field decreases by a factor of 
~3 for those regions with deeper ACS data (the cyan "PI" region shown on Figure [T]|. The 
effect of these additional data is rcfiected in the numbers we quote in this table. 



larger than our fiducial model. 

As Figure [TS] illustrates, the observed size distribution is in excellent agreement with that predicted by our fiducial 
model. The fit to the other size models is less good, particularly the model with the smaller size distribution (the 
disagreement ranges from '^l-2tT), but the uncertainties are large. The poor fit of the observations to the models with 
smaller sizes suggest that this issue is unlikely to bias us towards shallower slopes. Deeper observations will be required 
to establish the size distribution more accurately to very low luminosities. We find no evidence that our model size 
distribution is too large, and therefore that our selection volumes are too low, as suggested by Grazian et al. (2011). 

Incidentally, it would appear that Grazian et al. (2011)'s claim is based on a mistaken understanding of our method- 
ology (shown in their Figures 2 and 4) . They suggest that we model all galaxies using spectroscopically confirmed (and 
therefore predominantly bright) z 4 galaxies from the HUDF. Regretably, this is a misunderstanding of our approach 
and leads them to the wrong conclusion about our results. We explicitly consider luminosity when modeling source 
sizes and structures, i.e., the structure of luminous z 7 galaxies are modeled with luminous z ^ 4 galaxies (from 
the HUDF) and the structure of lower luminosity z 7 galaxies are modeled with lower luminosity z ^ 4 galaxies 
(from the HUDF). Since lower luminosity 2: ^ 4 galaxies are typically much smaller than more luminous galaxies, we 
model the sizes of faint galaxies a.t z ^ 7-8 with much smaller galaxies than has been suggested. Moreover, as shown 
in Figure [T8l these sizes appear to be in excellent agreement with the observations. 

B. EXPECTED CONTAMINATION RATES 

In §3.5, we consider many different sources of contamination for our z ^ 7 and z ^ 8 selections, including (1) 
low- mass stars (T dwarfs), (2) transient sources (SNe), (3) AGN, (4) low redshift galaxies, (5) photometric scatter, 
and (6) spurious sources. After careful consideration of the source characteristics and simulations, we find that the 
only meaningful source of contamination for our z > 7 samples are lower-redshift galaxies photometrically scattering 
into our selection. In Table [9l we give the expected contamination rates for the z ^ 7 and z 8 selections as a function 
of magnitude for all four fields considered here. 
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Fig. 19. — Apparent detection significance xtpt (Appendix D: blue lines) expected for z ~ 7 galaxy candidates from simulations for three 
different data sets considered here: ERS {left), HUDF09 (center), and HUDF09-1 (right). Since z ~ 7 galaxies should not be detected 
in the optical, these distributions (blue lines) are, as expected, approximately centered on 0, with ~50% above and ~50% below (due 
to noise). Also presented here (red lines) are the distributions expected for lower-redshift contaminants in our fields (after excluding the 
obvious contaminants detected at >2a in >1 bands). The latter distributions (red lines: for the low-redshift contaminants) are estimated 
through the photometric scatter simulations described in §3.5.5. A dramatic demonstration of the added value from our full Xopt approach 
can be seen in these figures. Note the large number of low-redshift galaxies we expect to satisfy our z ~ 7 LBG criterion that are not 
excluded by a 2a optical non-detection criterion. This is effectively illustrated by the area under the red line. For fields with optical data 
of limited depths (e.g., the ERS or HUDF09-1 fields), this number can be comparable to the number of bona- fide z ~ 7 galaxies. To 
effectively exclude these low redshift contaminants, we derive Xopt values for sources in our selections and then only select sources with 
small values of Xopt> '-^-i Xopt < xfim' ^^'^ HUDF09 and HUDF09-1 fields, the limiting values of Xopt shown in this diagram are for 
sources selected ~0.7 mag brighter than the selection limit (~28.5 AB mag but this depends on the field). The dashed vertical lines show 
these Xopt limits. The sources which exceed the Xopt limits (but otherwise meet our selection criteria) are given in Tables [TSl I19l and l20l 



C. EFFECTIVE SELECTION VOLUMES 

To derive the rest-frame UV LFs from the observed surface densities of z 7-8 candidates, we require a good 
estimate of the effective selection volume for our selections, as a function of magnitude. To obtain these estimates, we 
construct mock catalogs for our search fields. Pixel-by-pixel images are generated for each of the sources in our catalogs, 
these sources are added to our observational data, and finally source catalogs are constructed from the simulated data 
in the same way as on the observations. Dropouts are then selected using the criteria described earlier in this section. 
In generating two-dimensional images for each source in our catalogs, we randomly select a similar luminosity z 4 
i?-dropout from the Bouwens et al. (2007) HUDF sample and then scale its size (physical) as (1 -I- z)~^ to match the 
observed size evolution to z ^ 8 (Oesch et al. 2010b; see also Ferguson et al. 2004; Bouwens et al. 2004, 2006). In 
setting the UV colors of the sources, we assume that the C/F-continuum slopes /? are distributed as found by Bouwens 
et al. (2010a): i.e., with a mean value of -2.0 at L*^3 (-21 AB mag), a mean value of —2.7 at 0.1L*^3 (-18.5 AB 
mag), and with a mean /3 between these two extremes for luminosities in the range 0.1-1LJ^3 (e.g., Bouwens et al. 
2009; Bouwens et al. 2010a; Finkelstein et al. 2010; Bunker et al. 2010; R.J. Bouwens et al. 2011, in prep). These 
simulations are performed for our z ^ 7-8 selection over the ERS observations and for all three of the IIUDF09 search 
fields. 

D. USE OF THE Xopt DISTRIBUTION TO CONTROL FOR CONTAMINATION 

The best way to control for contamination in the z > 7 samples is to take advantage of the optical imaging 
observations over these samples. Typically, this is done by imposing an optical non-detection (< 2a) criterion in one 
or more bands or in a stack of the optical data - since high-redshift sources should not show a detection at optical 
wavelengths. Unfortunately, for the present z > 7 selections, the optical observations available over our WFC3/IR 
fields are sufficiently shallow (in some cases) that it is not clear that the above procedure will be effective in excluding 
all low redshift contaminants from our selections (or even the majority). This concern is particularly germane for the 
ERS, IIUDF09-1, and HUDF09-2 fields - where the optical data are not appreciably deeper than the near-IR data. 

D.l. Definition of xtpt 

To make maximal use of the optical imaging data that exist over the HUDF09 and GOODS fields, we clearly need 
to use a more sophisticated procedure. To this end, we define the quantity Xopt that effectively combines information 
in all the optical imaging data. We define xlpt 



(Dl) 
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where fi is the flux in band i in our smaller scalable apertures, ai is the uncertainty in this flux, and SGN(/i) is equal 
to 1 if /j > and —1 if /,; < 0. Included in the computed Xopt values are the i?435, Vqqq, 1775 band observations for the 
z ^ 7 zgso-dropout selection and the -B435, Vgoei *775i a-nd zgso band observations for the z 8 Yios-dropout selection. 

D.2. Xopt Distribution for z>7 Galaxies and Low-Redshift Contaminants 

How shall we use the measured x^^j values for sources in our selection to exclude low-redshift contaminants? To 
answer this question, we must know the Xopt distribution expected for bona-fide z >7 galaxies and also that distribution 
expected to come from possible low-redshift contaminants. The Xopi distribution for z ^ 7 candidates is calculated 
assuming no optical flux in z ~ 7 candidates and ideal noise properties. Meanwhile, the x^pj distribution for the 
low-redshift contaminants is calculated in the same way as we calculate contamination from photometric scatter in 
§3.5, i.e., starting with the color distribution for an intermediate magnitude subsample of real sources in our data and 
then adding noise to match the fluxes of faint sources in our search flelds. While discussing the use of x^^i to control 
for contamination in our z ^ 6 selections, we will assume that we have already removed sources that are detected at 
> 2(7 in individual images. 

We show the model Xopt distributions in Figure [19] for the z ^1 selections we perform over the ERS, HUDF09, and 
HUDF09-1 fields (which represent a good sampling of the selections considered here). For the HUDF09 and HUDF09-1 
fields, these distributions are shown ^^0.7 mag from the approximate selection limit for each field. As Figure fT9l shows, 
the Xopt values for z > 7 galaxies are centered on and have lower values than for low-redshift galaxies. However, due 
to the effects of noise, both z > 7 galaxies and low-redshift contaminants are spread over a wide range in Xopt values. 

Indeed, it is striking how low the Xopt values are expected to become for some low-redshift galaxy contaminants that 
would appear in our selections. At relatively low Xopt values of 3, for example, the total number of contaminants in 
our ERS or HUDF09-1 selections (red lines) is expected to exceed the number of bona-fide z > 7 galaxies {blue lines: 
Figure fT9)) . A few contaminants are even expected to have observed Xopt values that are negative, particularly in those 
fields with the shallowest optical observations (i.e., the ERS or HUDF09-1 data). The prevalence of contaminants to 
even very low values of Xopt underlines the importance of adopting very conservative optical non-detection criteria in 
our selection. 

D.3. Use of the measured Xopt selecting z > 7 galaxies 

We only select galaxies to a value of Xopt where we expect the number of z 7 galaxies to be greater than the 
number of contaminants. This critical value Xum depends on both the depth of the optical data and the magnitude 
of the source - since we would expect our fainter selections to admit more contamination than our brighter selections. 
We have found that xlm values of 5, 3, and 5, respectively, worked well for our HUDF09, HUDF09-1, and HUDF09-2 
data sets, respectively, at brighter ( J125, -ffieo ^ 28.5) magnitudes; however, at fainter magnitudes ( J125, -ffieo ^ 29.2), 
xfim needed to be 2.5, 1.5, and 2.5, respectively, for the same fields. For sources with magnitudes between 28.5 and 
29.2, the upper limits on Xopt interpolation between these extremes. For the small fraction (^-^30%) of sources 

in the HUDF09-2 field without additional optical/ACS coverage from our HUDF09 program (those not within the PI 
field shown in Figure 1), it was appropriate to adopt the same Xopt hmits as for the HUDF09-1 field. 

We derive Xopt values for our z > 7 candidates in three different apertures: 0.35"-diameter circular apertures, small 
0.18"-diameter circular apertures, and Kron-style scalable apertures (Kron factor of 1.2). The use of three different 
apertures allows us to be much more comprehensive in our identification of possible optical light in sources (which may 
be spatially very compact or extended) and therefore rejection of these sources as low-z contaminants. No smoothing 
of the ACS optical images is performed, prior to making these formal flux measurements. This ensures that we obtain 
the highest possible S/N for these flux measurements - effectively taking advantage of the sharper PSF in the ACS 
images. 

D.4-- Validation of the technique 

The Xopt distributions observed for our samples also provide a convenient way of ascertaining the actual number of 
bona-fide z > 7 galaxies in our search fields and therefore the contamination rate (by inference). Looking at Figure [TOl 
we see that an exceedingly small number of contaminants are expected to have xtpt ^ from the simulations, and 
so essentially all of sources satisfying our LEG criteria with Xopt < bona-fide z~7-8 galaxies. However, we 
also know the xlpt < requirement introduces an incompleteness of ~50% (simply because ^50% of the sources will 
have positive values of Xopt ^^d ^50% negative values - given our definition of Xopt)- Therefore the actual number 
of bona-fide z > 7 candidates is approximately double the number of candidates with Xopt < 0. Moreover, since we 
quantify the Xopt distribution for each search field in three different apertures (0.35"-diameter circular, 0.18"-diameter 
circular, scalable Kron), this provides us with three different estimates of the actual number of z > 7 candidates. 

Reassuringly, the number of high-redshift candidates in our samples is in good agreement with the numbers derived 
from the above argumentation. Multiplying the number of z ~ 7 candidates with xlpt < by 2, we expect 28.7, 15.3, 
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15.3, and 8.7 z ~ 7 candidates in the HUDF09, HUDF09-1, HUDF09-2, and ERS fields, respectively, vs. the 29, 17, 14, 
and 13 candidates in our observed samples. Similarly, we expect 28, 15.3, 10.7, and 6 z 8 candidates, respectively, 
in the above fields vs. the 24, 14, 15, and 6 candidates in our observed samples. 

D.5. Limitations of high-redshift selections not making use of quantities like Xopt 

It is instructive to contrast the present criteria for rejecting low-redshift contaminants with other frequently-used 
methods. Perhaps the most common of these is one relying just upon the 2a detections of the sources in individual 
optical bands for removing low-redshift contaminants from high-redshift selections (e.g., Bouwens et al. 2007; Wilkins 
et al. 2011). The difficulty with this approach is that it relies upon the optical data being very deep in each individual 
band (which is rarely the case). 

As an example of the limitations of this approach we found that had we only used a 2tT-optical rejection criterion in 
selecting z ~ 7 and z ^ 8 galaxies, a large fraction (^-^30-40%) of our z ^ 7-8 samples over the HUDF09-1, HUDF09-2, 
and ERS fields would have consisted of contaminants. Such large contamination rates significantly reduce the value 
of the samples for quantitative work. 

The red lines in Figure \W\ show the Xopt distribution expected for low-redshift contaminants not detected in any 
optical band at >2a. The area under the red curves in Figure [TOl is comparable to that under the blue curves for data 
sets like HUDF09-1 or the ERS fields. Clearly, contamination from low-redshift sources can be quite substantial and 
have a significant impact in fields with shallower optical data (e.g., recall the discussion of contamination in the Wilkins 
et al. 2010, 2011 sample in §4.1.2 and §4.1.3). This is perhaps an even more important concern than uncertainties 
about the size distribution at faint magnitudes (see, e.g., Bouwens et al. 2004; Grazian et al. 2011). Careful use of 
measures like Xopt required to limit contamination (or correct for it statistically in a robust way). 

E. Z~7-8 SAMPLES 

In this section, we provide detailed catalogs and postage stamps for our z ^ 7 and z 8 samples from the three 
ultra-deep WFC3/IR HUDF09 fields and ERS fields. 

A description of the HUDF09 z ^ 7 and z ^ 8 selections is given in §3.2. A similar description of the ERS z ~ 7 
and z ^ 8 selections is given in §3.4. The coordinates, total magnitudes, colors, and sizes measured for our z 7 
candidates over our ultra-deep IIUDF09 and ERS fields are provided in Tables [TOlfTSl Small postage stamp cutouts of 
the z ^ 7 candidates are given in Figure [20l The properties of our z 8 candidates in the three ultra-deep HUDF09 
and ERS fields are given in Tables [HllTTl Postage stamps of the z '-^ 8 candidates are given in Figure [2TJ The selection 
of these candidates is described in §3.2 and §3.4. 

Sources that just missed our HUDF09 or ERS z > 7 selections as a result of flux in the optical bands (Xopt) or colors 
that just missed our selection criteria are given in Tables [T8ll20l (see Appendix D). 




Fig. 20. Coaddcd yi05/>098 + J125 + Hum images (2.34" x 2.34") of the 73 2 ~ 7 candidates identified over our HUDF09 and ERS 
searcli fields. Eacli candidate is annotated with its source ID and Ji25-band magnitude. The candidates are ordered in terms of their 
Ji25-band magnitude, with those from the HUDF09 program given first and the ERS observations second. 
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Fig, 21,— Coaddcd J125 + H160 images (2,34" X 2,34") of the 59 2 ~ 8 identified over our HUDF09 and ERS search fields. Each candidate 
is annotated with its source ID and Hieo-band magnitude. The candidates are ordered in terms of their ifigo-band magnitude, with those 
from the HUDF09 program given first and the ERS observations second. 
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TABLE 10 

Z ~ 7 CANDIDATES IDENTIFIED IN THE ULTRA-DEEP HUDF09 WFC3/IR OBSERVATIONS OVER THE HUDF. 



Object ID 




R.A. 




Dec 




-,'^125 




2850 - Yios'^ 


Yl05 - ^125" 


Ji25 — Hiao 






UDFz- 


-42566566 


03: 


;32 


:42, 


,56 


-27:46:56, 


.6 


25.9±0, 


.0 


1.6±0.1 


0.2±0.0 


-O.liO.O 


0. 


,26" 


1,2,3,4,6,8,9,10,11,12 


UDFz- 


-42567S14 


03: 


:32 


:42, 


,56 


-27:47:31, 


,4 


27.0±0, 


.1 


1.7±0.5 


0.3±0.1 


-O.liO.l 


0. 


.19" 


1,3,4,6,8,9,10,11 


UDFz- 


-42196278 


03: 


:32 


:42, 


19 


-27:46:27, 


,8 


27.6±0, 


.1 


0.9±0.2 


O.OiO.l 


-O.litO.l 


0. 


17" 


8,11,12 


UDFz- 


-'^67775'^6 


03: 


:32 


:36. 


,77 


-27:47:53, 


,6 


27.8±0, 


.1 


1.0±0.2 


-O.liO.l 


-O.litO.l 


0. 


14" 


6,8,12 


UDFz- 


-44 '^06452 


03: 


:32 


:44. 


,30 


-27:46:45, 


,2 


27.9±0, 


.1 


0.8±0.3 


0.1±0.2 


-0.3±0.1 


0. 


18" 


8,11,12 


UDFz- 


-■^744651 2 


03: 


:32 


37 


,44 


-27:46:51, 


,2 


27.9±0, 


.1 


1.1±0.3 


-O.liO.l 


-0.2±0.1 


0. 


14" 


6,8,10,11,12 


UDFz- 


-■^9576564 


03: 


:32 


39 


,57 


-27:46:56, 


,4 


27.9±0, 


.1 


1.2±0.3 


O.OiO.l 


-O.litO.l 


0. 


16" 


6,8,10,11,12 


UDFz- 


-4n5fi64'^5 


03: 


:32 


40, 


,56 


-27:46:43, 


,5 


27.9±0, 


.1 


1.3±0.5 


-O.OiO.l 


-0.2±0.1 


0. 


19" 


6,8,9,10,11,12 


UDFz- 


-44746449 


03: 


:32 


:44, 


,74 


-27:46:44, 


,9 


28.0±0, 


.1 


>1.6 


0.1±0.2 


-0.0±0.2 


0. 


,27" 


8,10,11 


UDFz- 


-'^6'^7716S 


03: 


:32 


:36. 


,37 


-27:47:16, 


,3 


28.0±0, 


.1 


1.1±0.3 


-O.liO.l 


-0.4±0.2 


0. 


15" 


6,8,9,10,11,12 


UDFz- 


-41 8261 1 2 


03: 


:32 


:41. 


,82 


-27:46:11, 


,2 


28.0±0, 


.1 


1.1±0.6 


0.0±0.2 


0.0±0.2 


0. 


,23" 


10,11,12 


UDFz- 


-42406242 


03: 


:32 


:42. 


,40 


-27:46:24, 


,2 


28.1±0, 


.1 


>2.0 


-0.1±0.2 


-0.2±0.2 


0. 


24" 


10,11 


UDFz- 


-41 597(144 


03: 


:32 


:41. 


,59 


-27:47:04, 


,4 


28.3±0, 


.1 


1.1±0.4 


-0.1±0.2 


-0.3±0.2 


0. 


15" 


12 


UDFz- 


-■^972621 S 


03: 


:32 


39 


,72 


-27:46:21, 


,3 


28.5±0, 


.1 


>1.9 


0.2±0.2 


-0.3±0.2 


0. 


12" 


6,8,9,10,11,12 


UDFz- 


-S7'^56245 


03: 


:32 


:37. 


,35 


-27:46:24, 


.5 


28.5±0, 


,2 


1.1±0.7 


-0.2±0.2 


-0.4±0.3 


0. 


22" 




UDFz- 


-41866322 


03: 


:32 


:41, 


,86 


-27:46:32, 


,2 


28.6±0, 


.2 


>1.4 


0.1±0.3 


-0.4±0.3 


0. 


19" 




UDFz- 


-34026504 


03: 


:32: 


:34. 


,02 


-27:46:50, 


.4 


28.7±0, 


,1 


>1.6 


0.3±0.2 


-0.6±0.3 


0. 


,12" 




UDFz- 


-37588002 


03: 


:32: 


:37, 


,58 


-27:48:00, 


,2 


28.7±0, 


.2 


>0.9 


0.2±0.4 


0.3±0.3 


0, 


,24" 


10 


UDFz- 


-36678068 


03: 


:32: 


:36, 


,67 


-27:48:06, 


.8 


28.9±0, 


.2 


>1.4 


0.2±0.3 


-0.3±0.3 


0, 


,14" 




UDFz- 


-37015533 


03: 


:32: 


:37, 


,01 


-27:45:53, 


.3 


28.9±0, 


,2 


0.7±0.6 


0.0±0.3 


-1.0±0.8 


0, 


,17" 




UDFz- 


-34566493 


03: 


:32: 


:34. 


,56 


-27:46:49, 


.3 


29.0±0, 


,2 


>1.0 


0.2±0.4 


0.1±0.3 


0. 


,19" 




UDFz- 


-39747451 


03: 


:32: 


:39, 


,74 


-27:47:45, 


.1 


29.0±0, 


,2 


>1.2 


0.1±0.3 


-0.1±0.3 


0, 


,14" 




UDFz- 


-40707346 


03: 


:32: 


:40, 


,70 


-27:47:34, 


.6 


29.0±0, 


.2 


>1.4 


-0.2±0.3 


-0.3±0.4 


0, 


,16" 




UDFz- 


-37546018 


03: 


:32: 


:37, 


,54 


-27:46:01, 


.8 


29.1±0, 


,2 


>0.9 


0.2±0.4 


-0.1±0.4 


0, 


,17" 




UDFz- 


-33247173 


03: 


:32: 


:33, 


,24 


-27:47:17, 


.3 


29.1±0, 


.2 


>1.0 


0.3±0.4 


-0.2±0.4 


0, 


,17" 




UDFz- 


-38586520 


03: 


:32: 


:38, 


,58 


-27:46:52, 


.0 


29.1±0, 


.2 


1.4±0.8 


0.1±0.3 


0.0±0.3 


0, 


,14" 




UDFz- 


-40376559 


03: 


:32: 


:40, 


,37 


-27:46:55, 


.9 


29.2±0, 


,2 


>1.3 


-0.2±0.4 


-0.3±0.4 


0, 


,15" 




UDFz- 


-38896589 


03: 


:32: 


:38, 


,89 


-27:46:58, 


.9 


29.2±0, 


,3 


>0.8 


0.2±0.5 


-0.1±0.4 


0, 


,16" 




UDFz- 


-38537518 


03: 


:32: 


:38, 


,53 


-27:47:51, 


.8 


29.4±0, 


,2 


>1.5 


0.1±0.3 


-0.6±0.3 


0. 


,09" 


6 



^ Lower limits on the measured colors are the la limits. Magnitudes are AB. 

Half-light radius calculated using SExtractor from the square root of image. The median difference between the half-light 
radii estimates in the first-year observations and that estimated from the full two-year observations are ~0.02". 
" References: [1] Bouwcns et al. (2004), [2] Bouwens & lUingworth (2006), [3] Labbe et al. (2006), [4] Bouwens et al. (2008), [5] 
Oesch et al. (2009), [6] Oesch et al. (2010a), [7] Bouwens et al. (2010b), [8] McLure et al. (2010), [9] Bunker et al. (2010), [10] Yan 
et al. (2010), [11] Finkelstein et al. (2010), [12] McLure et al. (2011) 



TABLE 11 

Z 7 CANDIDATES IDENTIFIED IN THE ULTRA-DEEP WFC3/IR OBSERVATIONS OVER HUDF09-1. 



Object ID 




R.A. 




Dec 




Jl25 




2850 - ^105" 


yi05 - J125'' 


J125 — ^^160 




ReF 


UDF091Z- 


-58980501 


03: 


:32: 


:58, 


,98 


-27:40: 


:50, 


.1 


26.3±0, 


.1 


1.2±0.2 


O.litO.l 


-O.litO.l 


0. 


,36" 




UDF091Z- 


-59731189 


03: 


:32: 


:59, 


,73 


-27:41: 


:18, 


.9 


27.0±0, 


.1 


>2.6 


0.2it0.1 


-O.OitO.l 


0. 


,24" 




UDF091Z- 


-08071156 


03: 


:33: 


:08, 


,07 


-27:41: 


:15, 


.6 


27.2±0, 


,1 


1.1±0.3 


-O.litO.2 


0.2it0.2 


0. 


,19" 




UDF091Z- 


-02421312 


03: 


:33: 


:02, 


,42 


-27:41: 


:31, 


.2 


27.3±0, 


,1 


2.5±0.7 


O.litO.l 


-O.OitO.l 


0. 


,16" 


2 


UDF091Z- 


-58530234 


03: 


:32: 


:58, 


,53 


-27:40: 


:23, 


,4 


27.4±0, 


.1 


>2.5 


0.2±0.2 


-O.litO.l 


0, 


,16" 


2 


UDF091Z- 


-56930504 


03: 


:32: 


:56, 


,93 


-27:40: 


:50, 


,4 


27.8±0, 


.1 


1.4±0.4 


0.1±0.2 


-0.5it0.3 


0, 


,16" 




UDF091Z- 


-59171517 


03: 


:32: 


:59, 


17 


-27:41: 


:51, 


,7 


27.9±0, 


.2 


1.8±0.8 


0.0±0.3 


0.1±0.2 


0, 


,22" 




UDF091Z- 


-57811065 


03: 


:32: 


:57, 


,81 


-27:41: 


:06, 


,5 


28.1±0, 


.2 


>1.6 


0.2±0.3 


0.0±0.3 


0, 


,19" 




UDF091Z- 


-02122003 


03: 


:33: 


:02, 


12 


-27:42: 


:00, 


,3 


28.2±0, 


.2 


>2.0 


-0.0±0.3 


0.2±0.2 


0, 


,16" 




UDF091Z- 


-57860514 


03: 


:32: 


:57, 


,86 


-27:40: 


:51, 


,4 


28.2±0, 


.2 


1.0±0.6 


0.4±0.4 


-0.9it0.5 


0, 


,16" 




UDF091Z- 


-07300599 


03: 


:33: 


:07, 


,30 


-27:40: 


:59, 


,9 


28.4±0, 


.2 


>1.6 


0.3±0.4 


-0.4it0.4 


0, 


,15" 




UDF091Z- 


-58740215 


03: 


:32: 


:58, 


,74 


-27:40: 


:21, 


,5 


28.4±0, 


.2 


1.1±0.5 


-0.0±0.3 


-O.litO.3 


0, 


,15" 




UDF091Z- 


-02822023 


03: 


:33: 


:02, 


,82 


-27:42: 


:02, 


,3 


28.6±0, 


.2 


>1.3 


0.4±0.4 


-0.4it0.4 


0, 


,13" 




UDF091Z- 


-57511299 


03: 


:32: 


:57, 


,51 


-27:41: 


:29, 


,9 


28.6±0, 


.2 


>1.5 


0.3±0.4 


-0.2it0.4 


0, 


,13" 




UDF091Z- 


-04582038 


03: 


:33: 


:04, 


,58 


-27:42: 


:03, 


,8 


28.6±0, 


.2 


1.1±0.6 


0.1±0.3 


-0.3it0.4 


0, 


,13" 




UDF091Z- 


-00471180 


03: 


:33: 


:00, 


,47 


-27:41: 


:18, 


,0 


28.6±0, 


.2 


1.3±0.7 


0.1±0.3 


0.0±0.3 


0, 


,12" 




UDF091Z- 


-03791123 


03: 


:33: 


:03, 


,79 


-27:41: 


:12, 


,3 


28.7±0, 


.2 


>1.8 


-0.1±0.3 


-0.4it0.4 


0. 


,13" 


1 



Lower limits on the measured colors are the la limits. Magnitudes are AB. 
^ Half-light radius calculated using SExtractor. See footnote b from Table [TOl 
" References: [1] Bouwens et al. (2008), [2] Wilkins et al. (2011) 



33 



TABLE 12 

Z ~ 7 CANDIDATES IDENTIFIED IN THE ULTRA-DEEP WFC3/IR OBSERVATIONS OVER HUDF09-2. 



Object ID 




R.A. 




Dec 




^125 




2850 - y'los'' 




J125 - H^&Q^ 




ReF 


UDF092Z-01320531 


03: 


:33;01, 


,32 


-27:50:53, 


.1 


26.5±0, 


1 


1.1±0.3 


0.2±0.2 


-O.litO.l 


0.25" 




UDF092Z-09770485 


03: 


:33:09, 


,77 


-27:50:48, 


.5 


26.9±0. 


,1 


1.4±0.5 


-O.OitO.l 


0.3±0.1 


0.26" 


1 


UDF092Z-09151554 


03: 


:33:09, 


15 


-27:51:55, 


.4 


27.0±0. 


,1 


>2.6 


-O.liO.l 


-O.OiO.l 


0.17" 


1,2 


UDF092Z-00801320 


03: 


:33:00, 


,80 


-27:51:32, 


.0 


27.8±0. 


,1 


0.8±0.3 


0.1±0.2 


-0.0±0.2 


0.16" 




UDF092Z-05401189 


03: 


:33:05, 


,40 


-27:51:18, 


.9 


27.9±0. 


1 


1.7±0.9 


0.3±0.2 


-0.4±0.2 


0.12" 


1,2 


UDF092Z-09710486 


03: 


:33:09, 


,71 


-27:50:48, 


.6 


28.1±0. 


1 


>1.5 


0.1±0.3 


-0.1±0.2 


0.18" 




UDF092Z-06571598 


03: 


:33:06, 


,57 


-27:51:59, 


.8 


28.4±0. 


1 


>2.1 


-0.1±0.2 


-0.3±0.3 


0.13" 




UDF092Z-00921119 


03: 


:33:00, 


,92 


-27:51:11, 


.9 


28.7±0. 


,2 


>1.4 


0.0±0.4 


-0.2±0.4 


0.15" 




UDF092Z-07420235 


03: 


:33:07, 


,42 


-27:50:23, 


.5 


28.7±0. 


,2 


1.1±0.5 


-0.3±0.3 


0.0±0.3 


0.16" 




UDF092Z-09901262 


03: 


:33:09, 


,90 


-27:51:26, 


.2 


28.8±0. 


,2 


>1.3 


-O.liO.3 


0.2±0.3 


0.12" 




UDF092Z-04412250 


03: 


:33:04, 


,41 


-27:52:25, 


.0 


28.8±0. 


,2 


1.2±0.9 


0.1±0.4 


-0.6±0.6 


0.15" 




UDF092Z-04760294 


03: 


:33:04, 


,76 


-27:50:29, 


.4 


28.8±0. 


2 


1.7±0.7 


-0.1±0.3 


-0.3±0.4 


0.11" 




UDF092Z-04651563 


03: 


:33:04, 


,65 


-27:51:56, 


.3 


28.9±0. 


2 


>1.4 


0.3±0.3 


-0.2±0.3 


0.11" 




UDF092Z-04110148 


03: 


:33:04, 


11 


-27:50:14, 


.8 


28.9±0. 


2 


1.2±0.9 


0.1±0.4 


-0.0±0.4 


0.13" 





Lower limits on the measured colors are the Icr Umits. Magnitudes are AB. 
^ Half-hght radius calculated using SExtractor. See footnote b from Table [TO] 
" References: [1] Wilkins ct al. (2011), [2] McLure et al. (2011) 



TABLE 13 

2; ~ 7 CANDIDATES IDENTIFIED IN THE DEEP WIDE- AREA ERS OBSERVATIONS. 



Object ID 




R.A. 




Dec 




J125 




2850 — •/l25^ 


■^098 - ^125" 


^125 — ^^160 




Ref= 


ERSz- 


■2240942138 


03: 


:32:24, 


.09 


-27:42:13, 


,8 


25.8±0, 


,1 


1.6±0.9 


0.4±0.3 


-0.1±0.2 


0.30" 


4,5,6 


ERSz- 


■2432842478 


03: 


:32:43, 


.28 


-27:42:47, 


,8 


25.9±0, 


,1 


1.2±0.4 


0.7±0.3 


O.litO.l 


0.24" 


1 


ERSz- 


■2354442550 


03: 


:32:35, 


.44 


-27:42:55, 


,0 


26.1±0, 


,1 


2.2±0.8 


1.1±0.3 


-O.litO.l 


0.16" 




ERSz- 


■2160041591 


03: 


:32:16, 


.00 


-27:41:59, 


,1 


26.3±0, 


.1 


1.6±0.7 


0.3±0.3 


-0.3it0.2 


0.26" 


6 


ERSz- 


■2150943417 


03: 


:32:15, 


.09 


-27:43:41, 


,7 


26.4±0, 


.1 


>2.2 


1.0±0.5 


-O.litO.2 


0.22" 




ERSz- 


■2068244221 


03: 


:32:06, 


.82 


-27:44:22, 


,1 


26.5±0, 


.1 


1.1±0.4 


-0.0±0.2 


-0.4it0.2 


0.15" 


6 


ERSz- 


■2295342044 


03: 


:32:29, 


.53 


-27:42:04, 


,4 


26.6±0, 


.1 


1.4±0.6 


0.3±0.3 


-O.litO.2 


0.17" 


3,4,5,6 


ERSz- 


■2111644168 


03: 


:32:11, 


.16 


-27:44:16, 


,8 


26.9±0, 


.1 


>2.0 


1.1±0.5 


-0.3it0.2 


0.14" 




ERSz- 


■2150242362 


03: 


:32:15, 


.02 


-27:42:36, 


,2 


27.0±0, 


.2 


>1.2 


0.4±0.5 


0.4it0.3 


0.20" 




ERSz- 


■2056344288 


03: 


:32:05, 


.63 


-27:44:28, 


,8 


27.0±0, 


,2 


1.4±0.7 


0.2±0.4 


0.3it0.2 


0.16" 




ERSz- 


■2161941497 


03: 


:32:16, 


.19 


-27:41:49, 


,7 


27.1±0, 


,2 


>1.4 


0.6±0.5 


-O.litO.3 


0.17" 


6 


ERSz- 


■2225141173 


03: 


:32:22, 


.51 


-27:41:17, 


,3 


27.2±0, 


,2 


1.1±0.5 


0.1±0.3 


O.litO.3 


0.13" 




ERSz- 


■2352941047 


03: 


:32:35, 


.29 


-27:41:04, 


,7 


27.6±0, 


,2 


1.1±0.9 


0.3±0.5 


0.2it0.3 


0.14" 





" Lower limits on the measured colors are the la limits. Magnitudes are AB. 
^ Half-light radius calculated using SExtractor. See footnote b from Table [TO] 

" References: [1] Bouwens et al. (2010c), [2] Hickey et al. (2010), [3] Castellano et al. (2010a), [4] Wilkins et al. 
(2010), [5] Wilkins et al. (2011), [6] McLure et al. (2011) 
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TABLE 14 

Z ~ 8 CANDIDATES IDENTIFIED IN THE ULTRA-DEEP WFC3/IR OBSERVATIONS IN THE HUDF. 



Object ID 




R.A. 




Dec 








no5 - 


-^125 — ^^^160 


r 


b 

hi 


Ref<: 


UDFy-39537174 


03: 


;32:39 


.53 


-27:47:17, 


.4 


26.5±0, 


.0 


0.7±0.1 


O.litO.O 


0. 


.32" 


1,3,6,8,9,10,11,13 


UDFy-38807071 


03: 


:32:38, 


.80 


-27:47:07, 


.1 


26.8±0, 


.0 


0.5±0.1 


O.litO.O 


0. 


15" 


1,2,3,4,5,6,8,9,10,11,13 


UDFy-44706443 


03: 


:32:44, 


.70 


-27:46:44, 


.3 


27.1±0, 


.1 


0.9±0.2 


-O.litO.l 


0. 


19" 


6,8,9,10,11,13 


UDFy-43136284 


03: 


:32:43, 


.13 


-27:46:28, 


.4 


27.6±0, 


.1 


0.5±0.2 


O.litO.l 


0. 


18" 


6,8,9,10,11,13 


UDFy-37218061 


03: 


:32:37, 


.21 


-27:48:06, 


.1 


27.6±0, 


.1 


0.6±0.2 


-O.OitO.l 


0. 


17" 


6,8,9,10,11,13 


UDFy-38125539 


03: 


:32:38, 


.12 


-27:45:53, 


.9 


27.9±0, 


.1 


1.9±0.7 


O.litO.l 


0. 


17" 


7,8,9,10,11,12,13 


UDFy-42876344 


03: 


:32:42, 


.87 


-27:46:34, 


.4 


28.0±0, 


.1 


0.8±0.2 


O.OitO.l 


0. 


16" 


7,8,9,10,11,12,13 


UDFy-37636014 


03: 


:32:37, 


.63 


-27:46:01, 


.4 


28.0±0, 


.1 


1.8±0.6 


-O.litO.l 


0. 


,22" 


7,8,9,10,11 


UDFy-37796001 


03: 


:32:37, 


.79 


-27:46:00, 


.1 


28.2±0, 


.1 


1.4±0.4 


-O.litO.l 


0. 


15" 


7,8,9,10,11,12 


UDFy-43076241 


03: 


:32:43, 


.07 


-27:46:24, 


.1 


28.3±0, 


.2 


1.4±0.8 


-O.litO.3 


0. 


,28" 


8,13 


UDFy-38356118 


03: 


:32:38, 


.35 


-27:46:11, 


.8 


28.5±0, 


.1 


0.5±0.2 


-0.2it0.1 


0. 


12" 


9,10,11,13 


UDFy-33126544 


03: 


:32:33, 


.12 


-27:46:54, 


.4 


28.5±0, 


.2 


0.7±0.4 


-O.litO.2 


0. 


14" 


8,9,12,13 


UDFy-39468075 


03 


:32:39 


.46 


—27:48:07, 


.5 


28.6±0, 


.1 


>1.7 


0.5it0.2 


0. 


12 




UDFy-43086259 


03: 


:32:43, 


.08 


-27:46:25, 


.9 


28.6±0, 


.2 


0.8±0.5 


-O.litO.2 


0. 


,20" 




UDFy-40328026 


03: 


:32:40, 


.32 


-27:48:02, 


.6 


28.6±0, 


.2 


1.2±0.6 


O.OitO.2 


0. 


,17" 


11 


UDFy-35517443 


03: 


:32:35, 


.51 


-27:47:44, 


.3 


28.9±0, 


.2 


>1.4 


0.2it0.3 


0. 


,15" 




UDFy-34626472 


03: 


:32:34, 


.62 


-27:46:47, 


.2 


29.0±0, 


.2 


0.9±0.8 


0.2it0.3 


0. 


,14" 




UDFy-42406550 


03: 


:32:42, 


.40 


-27:46:55, 


.0 


29.0±0, 


.2 


1.1±0.7 


-O.OitO.3 


0. 


,16" 




UDFy-43086276 


03: 


:32:43, 


.08 


-27:46:27, 


.6 


29.1±0, 


.2 


1.2±0.5 


-0.4it0.2 


0. 


,12" 


7,8,10,11,12 


UDFy-33436598 


03: 


:32:33, 


.43 


-27:46:59, 


.8 


29.2±0, 


.2 


>1.6 


0.2it0.3 


0. 


,13" 


14 


UDFy-36816421 


03: 


:32:36, 


.81 


-27:46:42, 


.1 


29.3±0, 


.3 


>1.3 


-0.2it0.4 


0. 


,18" 




UDFy-39347255 


03: 


:32:39, 


.34 


-27:47:25, 


.5 


29.4±0, 


.2 


0.8±0.7 


0.2it0.4 


0. 


,11" 




UDFy-39106493 


03: 


:32:39, 


.10 


-27:46:49, 


.3 


29.4±0, 


.3 


0.9±0.5 


-0.7it0.3 


0. 


,19" 




UDFy-39505451 


03: 


:32:39, 


.50 


-27:45:45, 


.1 


29.4±0, 


.4 


l.Oil.O 


-0.2it0.5 


0. 


,15" 





^ Lower limits on the measured colors are the Icr Umits. Magnitudes are AB. 
^ Half-light radius calculated using SExtractor. See footnote b from Table [Tol 

" References: [1] Bouwens et al. (2004), [2] Bouwens & lUingworth (2006), [3] Labbe et al. (2006), [4] Bouwens 
et al. (2008), [5] Oesch et al. (2009), [6] Oesch et al. (2010a), [7] Bouwens et al. (2010b), [8] McLure et al. 

(2010) , [9] Bunker et al. (2010), [10] Yan et al. (2010), [11] Finkelstein et al. (2010), [12] Lorenzoni et al. 

(2011) , [13] McLure et al. (2011), [14] Bouwens et al. (2011) 



TABLE 15 

Z ~ 8 CANDIDATES IDENTIFIED IN THE ULTRA-DEEP WFC3/IR OBSERVATIONS OVER HUDF09-1. 



Object ID 




R.A. 




Dec 










. — Jl25^ 


J125 - 




UDF091y-02059566 


03: 


:33: 


;02, 


.05 


-27:39: 


:56, 


.6 


27.1it0. 


,1 


1, 


.OitO.3 


O.litO.l 


0, 


.19" 


UDF091y-02812107 


03: 


:33: 


:02, 


.81 


-27:42: 


:10. 


.7 


27.2it0. 


,1 


0, 


.7it0.3 


0.2±0.1 


0, 


.19" 


UDF091y-55761063 


03: 


:32: 


:55, 


.76 


-27:41: 


:06. 


.3 


27.4it0. 


,1 


0, 


.5it0.2 


0.2±0.1 


0. 


.14" 


UDF091y-56451003 


03: 


:32: 


:56, 


.45 


-27:41: 


:00. 


.3 


27.6it0. 


,1 


0, 


.6it0.2 


-O.litO.2 


0, 


.19" 


UDF091y-59372014 


03: 


:32: 


:59, 


.37 


-27:42: 


:01. 


.4 


27.6it0. 


,1 


0, 


.6it0.3 


0.1±0.2 


0, 


.19" 


UDF091y-02020311 


03: 


:33: 


:02, 


.02 


-27:40: 


:31. 


.1 


28.2it0. 


,3 


1, 


.4itl.0 


-O.OitO.4 


0, 


.20" 


UDF091y-02992182 


03: 


:33: 


:02, 


.99 


-27:42: 


:18. 


.2 


28.3it0. 


,2 


0, 


.8it0.5 


-O.OitO.3 


0, 


.15" 


UDF091y-04340518 


03: 


:33: 


:04, 


.34 


-27:40: 


:51. 


.8 


28.4it0. 


,2 


0, 


.8it0.6 


0.2±0.3 


0, 


.14" 


UDF091y-02081326 


03: 


:33: 


:02, 


.08 


-27:41: 


:32. 


.6 


28.4it0. 


,3 




>1.6 


-0.2it0.4 


0, 


.20" 


UDF091y-03721514 


03: 


:33: 


:03, 


.72 


-27:41: 


:51. 


.4 


28.4it0. 


,3 


1, 


.OitO.5 


-0.4it0.3 


0, 


.18" 


UDF091y-00531465 


03: 


:33: 


:00, 


.53 


-27:41: 


:46. 


.5 


28.5it0. 


,2 


0, 


.5it0.4 


0.2±0.3 


0, 


.12" 


UDF091y-03021569 


03: 


:33: 


:03, 


.02 


-27:41: 


:56. 


.9 


28.5it0. 


,2 


1, 


.3it0.6 


-0.4it0.3 


0, 


.14" 


UDF091y-01480311 


03: 


:33: 


:01, 


.48 


-27:40: 


:31. 


.1 


28.8it0. 


,3 


1, 


.2itl.0 


0.2±0.4 


0, 


.13" 


UDF091y-05391244 


03: 


:33: 


:05, 


.39 


-27:41: 


:24. 


.4 


29.1it0. 


,4 




>1.4 


-0.3it0.4 


0. 


.12" 



See footnote b in Table [TT] 
^ Lower limits on the measured colors are the Icr limits. Magnitudes are AB. 
Half-light radius calculated using SExtractor. See footnote b from Table [Tol 
References: [1] Lorenzoni et al. (2011) 
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TABLE 16 

2: ~ 8 CANDIDATES IDENTIFIED IN THE ULTRA-DEEP WFC3/IR OBSERVATIONS OVER HUDF09-2. 



Ohiprt in 




R. 


A. 




Dec 








^ lOZ 


T ^ 
1 "^125 


^125 ^160 






UDF092y-03781204<i 


03: 


:33: 


:03. 


78 


-27: 


:51: 


:20, 


,4 


26, 


.1±0, 


.0 


0, 


.6±0.1 


O.liO.O 


0, 


,20" 


1,3 


UDF092y-03751 ige'^ 


03: 


:33: 


:03. 


75 


-27: 


:51: 


:19, 


,6 


26, 


.3±0, 


.0 


0, 


.7±0.1 


-O.OitO.l 


0, 


,23" 


3 


UDF092y-07580550 


03: 


:33: 


:07. 


58 


-27: 


:50: 


:55, 


,0 


27, 


.1±0, 


,1 




>2.4 


0.3±0.1 


0, 


,22" 


3 


UDF092y-04640529'^ 


03: 


:33: 


:04. 


64 


-27: 


:50: 


:52, 


.9 


27, 


.2±0, 


,1 


0, 


.8±0.2 


— O.litO.l 


0, 


.15" 


3 


UDF092y-09661163 


03: 


:33: 


:09. 


,66 


-27: 


:51: 


:16, 


.3 


27, 


.2±0, 


,1 


1, 


.3±0.6 


0.5±0.2 


0, 


,24" 


2 


UDF092y-06321217 


03: 


:33: 


:06. 


,32 


-27: 


:51: 


:21, 


.7 


27, 


.7±0, 


.1 


0, 


.6±0.2 


O.OitO.l 


0, 


,15" 




UDF092y-09611126 


03: 


:33: 


:09. 


,61 


-27: 


:51: 


:12, 


.6 


27, 


.9±0, 


.2 


0, 


.8±0.5 


0.1±0.2 


0, 


,22" 




UDF092y-06970279 


03: 


:33: 


:06. 


,97 


-27: 


:50: 


:27, 


.9 


28, 


.0±0, 


.2 




>1.7 


-O.litO.3 


0, 


.22" 




UDF092y-09291320 


03: 


:33: 


:09. 


,29 


-27: 


:51: 


:32, 


.0 


28, 


.2±0, 


.2 


0, 


.8±0.5 


-0.2±0.3 


0, 


.20" 




UDF092y-03391003 


03: 


:33: 


:03. 


,39 


-27: 


:51: 


:00, 


.3 


28, 


.2±0, 


.2 


1, 


.2±0.6 


0.1±0.2 


0, 


,14" 




UDF092y-04242094 


03: 


:33: 


:04. 


,24 


-27: 


:52: 


:09, 


.4 


28, 


.3±0, 


.1 


1, 


.0±0.6 


0.5±0.2 


0, 


,11" 




UDF092y-03811034 


03: 


:33: 


:03. 


,81 


-27: 


:51: 


:03, 


.4 


28, 


.3±0, 


.2 


0, 


.5±0.4 


0.1±0.2 


0, 


,16" 




UDF092y-06391247 


03: 


:33: 


:06. 


,39 


-27: 


:51: 


:24, 


.7 


28, 


.3±0, 


.2 


0, 


.7±0.3 


-0.5±0.2 


0, 


.16" 




UDF092y-08891430 


03: 


:33: 


:08. 


,89 


-27: 


:51: 


:43, 


.0 


28, 


.4±0, 


.2 


0, 


.5±0.5 


0.4±0.3 


0, 


,12" 




UDF092y-01971379 


03: 


:33: 


:01. 


,97 


-27: 


:51: 


:37, 


.9 


28, 


.4±0, 


.2 


0, 


.8±0.6 


0.0±0.3 


0, 


,20" 




UDF092y-07830324 


03: 


:33: 


:07. 


,83 


-27: 


:50: 


:32, 


.4 


28, 


.7±0, 


.3 




>1.1 


0.0±0.4 


0, 


,16" 





^ Lower limits on the measured colors are the Icr limits. Magnitudes are AB. 

^ Half-light radius calculated using SExtractor. See footnote b from Table [TOl 

= References: [1] Wilkins et al. (2011), [2] Lorenzoni et al. (2011), [3] McLure et al. (2011) 

All three sources here are within 30" of each other (two being within 3" of each other), have consistent 
^105 ~ J125 colors (and thus redshifts), and are among the brightest z ~ 7.5 galaxies in the field. They 
would appear to be part of a possible overdensity. 



TABLE 17 

Z ~ 8 CANDIDATES IDENTIFIED IN THE DEEP WIDE- AREA ERS OBSERVATIONS. 



Object ID 


R.A. 


Dec 


^^160 


2^850 - ^125'' 


Vo98 - ^125'' 


^125 — -f^ieo 




Ref 


ERSY-2306143041'^ 


03:32:30.61 


-27:43:04.1 


26.0±0.1 


2.0±0.8 


1.4±0.7 


0.3±0.2 


0.23" 




ERSY-2399642019 


03:32:39.96 


-27:42:01.9 


26.4±0.2 


>1.4 


>1.5 


0.3±0.2 


0.22" 




ERSY-2251641574 


03:32:25.16 


-27:41:57.4 


26.8±0.2 


>1.7 


>1.5 


O.OibO.3 


0.21" 




ERSY-2029843519 


03:32:02.98 


-27:43:51.9 


26.8±0.2 


>1.7 


>1.6 


0.1±0.3 


0.14" 


1,2 


ERSY-2377942344'^ 


03:32:37.79 


-27:42:34.4 


26.8±0.2 


>1.9 


>1.7 


0.2±0.2 


0.14" 




ERSY-2354441327 


03:32:35.44 


-27:41:32.7 


27.2±0.2 


>1.6 


>1.7 


-0.1±0.3 


0.14" 


2 



Lower limits on the measured colors are the la limits. Magnitudes are AB. 
^ Half-light radius calculated using SExtractor. See footnote b from Table [TOl 
" References: [1] Lorenzoni et al. (2011), [2] McLure et al. (2011) 

Source is formally "detected" at ~ Icr in the Veoe band and ~ 1.3cr in the zg5o band, so is more likely to be at 
low-redshift than the other candidates 

Source has Hieo — 3.6/i and Hieo — 4.5/i colors of ^^2 mag, which is quite red for a 2 ~ 8 galaxy and may indicate 
it has a redshift z ~ 2 — 3. 
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TABLE 18 

Possible candidates over the HUDF09 fields that did not meet our selection criteria.'' 







rt. 


A 




Dec 




Jl25 




^85C 


1 — J 105 


J^IOS — '.'125 


Jl25 — -"160 






UDF092Z-01101160 


03; 


;33; 


:01, 


.10 


-27; 


;51: 


:16, 


.0 


25, 


.9±0, 


.0 





.8±0.1 


0.2±0.1 


-O.OitO.O 


0. 


.26" 


4 


UDF091Z-56701076 


03: 


:32: 


;56, 


.70 


-27; 


:41: 


:07, 


,6 


26, 


,7±0, 


,1 


1 


.4±0.3 


0.4±0.1 


-O.litO.l 


0, 


,20" 


3 


UDF091Z-59481461 


03: 


:32: 


:59, 


,48 


-27; 


:41: 


;46, 


1 


27, 


.6±0, 


1 


1 


.1±0.5 


0.1±0.2 


-0.6±0.3 


0, 


,25" 




UDF091Z-59270343 


03: 


:32: 


:59, 


,27 


-27; 


:40: 


:34, 


,3 


27, 


,8±0, 


1 




>2.0 


0.2±0.2 


-0.2±0.2 


0, 


,19" 




UDFz-37296175 


03: 


:32: 


:37, 


,29 


-27; 


:46: 


;17, 


,5 


28, 


,0±0, 


1 




>1.7 


0.1±0.2 


0.2±0.2 


0, 


,28" 


2 


UDF091Z-58830381 


03: 


:32: 


:58, 


,83 


-27; 


:40: 


:38, 


,1 


28, 


,1±0, 


1 


1 


.0±0.4 


0.1±0.2 


0.1±0.2 


0, 


,15" 




UDF092Z-07961350 


03: 


:33: 


:07, 


,96 


-27; 


:51: 


:35, 


,0 


28, 


,1±0, 


,1 


1 


.7±0.7 


0.1±0.2 


O.litO.l 


0, 


,12" 




UDFz-36916516 


03: 


:32: 


:36, 


,91 


-27; 


:46: 


:51, 


,6 


28, 


,2±0, 


,2 




>1.1 


0.2±0.4 


0.0±0.3 


0, 


,44" 




UDFz-36467324 


03: 


:32: 


:36, 


,46 


-27; 


;47: 


;32, 


,4 


28, 


,3±0, 


1 


1 


.6±0.6 


-0.1±0.2 


0.4±0.1 


0, 


,16" 


1,2 


UDF092Z-08551377 


03: 


:33: 


:08, 


,55 


-27; 


:51: 


:37, 


,7 


28, 


,3±0, 


,2 




>1.2 


-0.1±0.3 


0.2±0.3 


0, 


,22" 




UDF091Z-04781529 


03: 


:33: 


:04, 


,78 


-27; 


:41: 


:52, 


,9 


28, 


,3±0, 


,2 




>1.5 


0.2±0.4 


-0.1±0.4 


0, 


,19" 




UDF091Z-59561338 


03: 


:32: 


:59, 


,56 


-27; 


:41: 


;33, 


,8 


28, 


,3±0, 


,2 


1 


.8±0.6 


-0.3±0.3 


-0.4±0.4 


0, 


,19" 




UDF091Z-01020006 


03: 


:33: 


:01, 


,02 


-27; 


:40: 


:00, 


,6 


28, 


,4±0, 


,2 


1 


.2±0.7 


-0.1±0.3 


-0.5±0.5 


0, 


,18" 




UDFz-42247087 


03: 


:32: 


:42, 


,24 


-27; 


:47: 


;08, 


,7 


28, 


,5±0, 


,2 


1 


.5±0.8 


-0.2±0.2 


0.4±0.2 


0, 


,20" 


2 


UDF091Z-59451483 


03: 


:32: 


:59, 


,45 


-27; 


:41: 


:48, 


,3 


28, 


,6±0, 


,2 




>1.4 


-0.1±0.4 


-0.1±0.4 


0, 


,16" 




UDF091Z-02071130 


03: 


:33: 


;02, 


,07 


-27; 


;41: 


;13, 


,0 


28, 


,7±0, 


,2 





.7±0.5 


-0.0±0.4 


0.0±0.4 


0, 


,15" 




UDF092Z-10561177 


03; 


;33; 


;10. 


,56 


-27; 


;51: 


:17. 


,7 


28, 


,7±0. 


,2 





.9±0.6 


-0.1±0.3 


0.1±0.3 


0. 


,15" 




UDF092Z-09690571 


03; 


;33; 


;09. 


,69 


-27: 


;50: 


:57. 


,1 


28, 


.7±0. 


,3 




>0.8 


0.1±0.5 


-0.0±0.4 


0. 


,18" 




UDFz-40506512 


03; 


;32; 


;40. 


,50 


-27; 


;46: 


:51. 


,2 


28, 


.8±0. 


,2 




>1.3 


0.3±0.3 


-0.6±0.3 


0. 


,14" 




UDF091Z-59842044 


03; 


;32; 


;59, 


,84 


-27; 


;42: 


:04. 


,4 


28, 


.8±0. 


,2 


1 


.4±0.9 


0.2±0.4 


0.1±0.3 


0. 


,14" 




UDF091Z-05751079 


03; 


;33; 


;05. 


,75 


-27; 


;41: 


:07. 


,9 


28, 


.9±0. 


,2 




>1.2 


0.3±0.5 


-0.0±0.4 


0. 


,13" 




UDF092Z-07192310 


03; 


;33; 


;07. 


,19 


-27; 


;52: 


:31. 


,0 


28, 


.9±0. 


,3 




>0.9 


0.2±0.5 


-0.0±0.4 


0. 


,15" 




UDF092Z-07582097 


03; 


;33; 


;07, 


,58 


-27; 


;52: 


:09. 


,7 


28, 


.9±0. 


,3 





.9±1.1 


0.3±0.6 


0.3±0.4 


0. 


,16" 




UDF092Z-06900317 


03; 


;33; 


;06, 


,90 


-27: 


;50: 


:31. 


,7 


29, 


.0±0. 


,2 




>1.4 


-0.2±0.4 


0.4±0.3 


0. 


,14" 




UDF092Z-06960354 


03; 


;33; 


;06, 


,96 


-27; 


;50: 


:35. 


,4 


29, 


.0±0. 


,3 




>0.8 


0.2±0.6 


0.2±0.4 


0. 


,14" 




UDF092Z-04992391 


03; 


;33; 


;04, 


,99 


-27; 


;52: 


:39. 


,1 


29, 


.0±0. 


,3 




>0.9 


0.1±0.5 


0.3±0.4 


0. 


,15" 




UDFz-41357078 


03; 


;32; 


;41, 


,35 


-27; 


;47: 


:07. 


,8 


29, 


.2±0. 


,2 





.9±0.9 


-0.1±0.4 


0.1±0.3 


0. 


,16" 




UDFz-37316420 


03; 


;32; 


;37, 


,31 


-27; 


;46: 


:42. 


,0 


29, 


.2±0. 


,2 


1 


.5±0.6 


-0.5±0.3 


-0.9±0.7 


0. 


,15" 





* EELch of the sources here satisfy the color criteria for our dropout selections, but show sUght detections in the Xopt images 
generated from the B435, Veoe, ^775 data. These sources therefore do not make it into our Lyman-Break galaxy samples. See §3.3 
and Appendix D for a discussion of Xopt '''^'^ selection criteria. While a few of these sources will correspond to 2 > 7 galaxies, 

the majority will not according to our simulations (Appendix D). For a few sources, the Xopt detection comes almost entirely from 
the i775-band, suggesting that some might nevertheless have redshifts slightly above 2~6. 
^ Lower limits on the measured colors are the Icr limits. Magnitudes are AB. 

Half-light radius calculated using SExtractor 
d References: [1] McLure et al. (2010), [2] Finkelstein et al. (2010), [3] Wilkins et al. (2011), [4] McLure et al. (2011). 



37 



TABLE 19 

Possible 2:~8 candidates over the HUDF09 fields that did not meet our selection criteria.'' 



Object ID 


R.A. 


Dec 


J125 





J125 — H\QQ 






ReP 


VLfr Uyzy-lUO ( liDD 


no . QQ . 1 n c;7 
Uo. oo. iU.O / 


97- f^T ■ 1 C< 
— Z / .01. iu.D 


97 Q-Un 1 
z i .oitu. i 


n K-Uo 9 

u.omu.z 


U.4mU.l 


1 c;" 
u.io 




u ur uyzy-u / ^yzuou 


nQ-QQ-07 zlQ 


97. c;9.0(^ 
— Z ( .OZ.UD.U 


97 ^-\-C\ 1 

z 1 . omu. i 


^1 Q 

^i.y 


f^-l-O 9 

u.omu.z 


90" 
u.zu 






UO.OO.UU-OO 


97-Af\-AA Q 


97 7-1-0 
Z 1 . i ^u.o 


u.omu.o 


9-1-0 A 

u. zmu.'i 


u.zo 




TTDTTflQlA/ '\Q77^A'if, 
uur^uyiy-oy* / ±^ou 


0'i-'^9-f^Q 77 


97. /LI ■A'i 

Z i . 'I: J_ . 'iO . U 


9S 0-1-0 9 
Zo.U^U.Z 


u.omu.o 


u.omu.o 


1 7" 
u. 1 1 




u ur uyzy-u ( uy ± idu 


no. Qo.f|7 fiQ 

Uo.OO.U ( .uy 


97. rri .1 n 
— Z i .01. iO.U 


9S 0-1-0 9 

Zo.udiu.z 


1 1 -1-0 A 


4-1-0 9 

— u.^imu.z 


1 Q" 

u. ly 


9 
z 


U Ur y-oOl ( 1 1 (u 


no. Qo. q c: 1 7 


97- /I 7- 1 7 
— Z / .4 ( . i i .U 


9Q 1 -Un 1 
zo.iitU.l 


U.omU.4 


9-l-A 9 

U.ZmU.Z 


9q" 

U.Zo 


1 
1 


u -Ur uy zy-uo / zui ly 


nq.qo.nc 79 
Uo.oO.UO. / Z 


97. rrn. i 1 Q 

— z / .ou. ± i.y 


90 1 _i_n 9 

ZO. l^U.Z 


u.omu.o 


Q-l-0 

u.omu.o 


90" 
u.zu 




U JJr y-4o4iOoDU 


no. 00. /I Q /1 1 
Uo.oz.4o.41 


— z ; .40.00. U 


zo.zitU. i 


1 n-un c; 
i.UmU.o 


n O-LA 

U.ZmU.Z 


n 1 
U.ly 


1 
i 


U JJr Uyzy-Uo ( zlUyO 


Uo. 00. Uo. ( z 


07. "1 .no 
— z ; .ol.Uy.o 


Zo.oitU.z 


U.ymU.O 


U.zmU.z 


A 1 c;" 
U.io 




TTFM?nno,r nQfiviri/ii 
UJJr Vi)Zy-voD ( iU41 


no.QQ.nQ P.1 
Uo.oo.Uo.D / 


— z ( .oi.U4.1 


zo.oitU.z 


1 AJ-A c; 
i.UmU.o 


n 1 -LA 
U.lmU.z 


A 1 i^" 

U.io 




u ur uyzy-±uoziuoo 


UO. 00. J-U.OZ 


97. ci .fiSi 
— Z 1 .Ol.Uo.D 


98 Q-l-O Q 

zo.omu.o 


u.ymu.y 


u. imu.4 


n on" 
u.zu 




u ur uy ly-oy luiulo 


uo.oz.oy.iu 


97-41 -m 8 

— Z 1 .'il.Ul.O 


Zo.omu.o 


1 4-i-n Q 
i.'rimu.y 


n 0-1-0 4 
— u.zmu.'i 


n on" 
u.zu 




U Ur y-4tl 1 o 1 o4t4l: 


(IQ.QO./II 7Q 
UO.OZ.4i. (£5 


— Z 1.41 .o4.4 


9Q yl-Un 1 


U.omU.o 


n A-i-A 
U.UmU.z 


u.iy 




U Ur UyZy-UZ i olOD4: 


Uo.oo.Uz. / 


— Z / .Ol.OD.4 


OQ /i-un 1 


U.ymU.O 


yi-un 
U.4mU.Z 


n 1 n" 
u.iu 




uur uy iy-Dyozzuio 


Uo.oz.oy.oz 


97-/19-01 
— Z i .4Z.Ui.O 


9fi c;-i-o q 
zo.O±U.o 


u.Dmu.o 


A O-UO A 

— u.umu.4 


A 1 Q" 
U. io 




u ur uyzy-uoDoUo / u 


no. 00. no CO 
UO. 00. Uo. Do 


97. c^0-'^7 
— Z i .OU.O / .U 


9Q (^4-0 9 

zo. omu. z 


1 94-0 Q 

i.zmu.y 


1 4-0 

u. imu.o 


"1 f^" 
U. 10 




TTFtPOQI-^^ 009/10/170 

u ur uy ly-uuz'^iU'^i / u 


no.qo.nn O/i 
UO.OO.UU.Z^ 


97-A0-A7 

Z ( .-yrU.'^: ( .U 


Zo.omu.o 


.>1.0 


9-1-0 A 

— u.zmu.-y: 


1 Q" 

u. ly 




T T'nTrOQ9-(^ 070009 1^ 

u ur uy zy-u i uyuzio 


nq.qq.n7 no 

UO.OO.U ( .uy 


97.crn-91 S 
— Z i .OU.Zl.o 


9S (^-1-0 

Zo.omu.o 


f^-l-0 A 

u.omu.'4: 


Pi-I-O Q 

— u.omu.o 


1 s" 
u.io 




TTTIPOQIa/ 009^^900/1 
uur uy -Ly-uuzozuu -i 


OQ-QQ-00 9^^ 

uo.oo.uu.zo 


97-49-00 A 
Z i .'4lZ.UU.'4: 


9S fi-l-0 

Zo.omu.o 


fi-kO 4 
U.DmU.4: 


— u.omu.o 


1 V 
u. 10 




kjur uy z y-uuo'^± -Luoo 


O'^-'^'^-OO RA 
uo.oo.uu.oi 


— 97-'^l -OK 

Z / . -L . UO. 


90 a-\-f\ q 
zo . umu. 


0.6ib0.4 


— 0.6di0.4 


0.18" 






O'^-Q'^-OiS 7Q 

L/O.OO.L/U. 1 U 


Z> ( .'±W.'JtJ.O 


9fi (s-hn 4 


>1.4 


— 0.4di0.5 


0.19" 




u ur y-oouyoio-^i 


nq.qo.qo nn 
UO. oZ.Oo. UU 


97. /I 0.1 r: A 


9S 7-1-0 9 

Zo. ( mu.z 


u.omu.o 


— u.omu.z 


1 q" 
u. ly 




U Ur y-oD4y ; 40o 


oq-q9-qfi /10 
Uo.oZ.v5D.4y 


97- A'7-A^ q 
— Z i .4 i .40.0 


9fi 7-1- 9 

Zo. ( mu.z 


1 9-1-1 

i.zmi.u 


Q-LA q 

u.omu.o 


1 c\" 
u. iy 




uur uy ly-uDz ± lu^D 


nq.qq.ni^ 91 
UO.OO.UO.Zl 


97-/1 "1-0/1 
— Z / .■^il.U'^i.D 


9S 7-1-0 

Zo. / mu.o 


u.omu.o 


9-1-0 A 
— U. ZzCU.-^i 


1 fi" 

U . 1 D 




u ur uyzy-uyooi^ ( ^ 


nq. qq. nn f;q 

UO. 00. uy. DO 


97. ./I 7 /i 


9S 7-1-0 Q 

Zo. ( mu.o 


1 1 -LI n 
1. imi .u 


94-0 A 

u.zmu.4 


1 4" 
U. 1^ 




u Ur uyzy-u t \)^i.^^o 


Uo.OO.U 1 .y* 


97. CI .AA c; 


9R S-l-O 9 

zo.omu.z 


u.omu.o 


Q-l-fl Q 

u.omu.o 


nil" 
u. 11 




Uur y-o ( oycsi^y 


Uo.oz.o 1 .oy 


97. /I a. 1 9 Q 
— z 1 .4o.lz.y 


OQ Q-UO 9 

zo.omu.z 


i.OmU.y 


n Q-UA Q 

— U.omU.o 


n 90" 

U.zU 






uo . oo . U X . O 1 


— 97'fi1 '(14 4 


00 e-Lfi Q 
zo.omu.o 


0.7ib0.8 


0.2ib0.4 


0.14" 




u ur uy ry uoy'iAiJOO 


uo . OO . uo . y^ 


—97-41 '^R '\ 


zo.omu.o 


0.9ib0.7 


O.OibO.4 


0.13" 




T TDFnQ2v-n '^S42nQ4 


V70. OO. L/O. Ot: 


_97-f^9-OQ 4 


9« 8+0 4 

iiO .O— 1— VJ.'± 


l.OibO.9 


— 0.2di0.4 


0.18" 




UDFnq2v-nsn6i 458 




— 27-51 -45 8 


28 9+0 3 

. C/— 1— W . vJ 


>1.4 


— 0.3di0.4 


0.15" 




UDFy-37517553 


03:32:37.51 


-27:47:55.3 


29.0±0.2 


>1.5 


0.0±0.3 


0.16" 




UDF092y-07220444 


03:33:07.22 


-27:50:44.4 


29.0±0.3 


>1.1 


0.0±0.4 


0.12" 




UDF092y-02271378 


03:33:02.27 


-27:51:37.8 


29.0±0.3 


0.6±0.6 


-0.2+0.4 


0.15" 




UDF092y-10231298 


03:33:10.23 


-27:51:29.8 


29.0±0.3 


0.8±0.6 


-0.3+0.4 


0.12" 




UDF092y-10121280 


03:33:10.12 


-27:51:28.0 


29.0±0.3 


0.9±0.6 


-0.4+0.4 


0.14" 




UDFy-39405459 


03:32:39.40 


-27:45:45.9 


29.0±0.3 


1.3±0.9 


-0.4+0.4 


0.29" 




UDF092y-09351539 


03:33:09.35 


-27:51:53.9 


29.0±0.4 


>1.2 


-0.1+0.5 


0.14" 




UDFy-43347121 


03:32:43.34 


-27:47:12.1 


29.1±0.3 


1.1±0.8 


-0.1+0.3 


0.17" 




UDFy-38396038 


03:32:38.39 


-27:46:03.8 


29.2±0.2 


1.0±0.7 


0.1+0.3 


0.14" 




UDFy-39807470 


03:32:39.80 


-27:47:47.0 


29.3±0.2 


>1.0 


0.3+0.4 


0.15" 




UDFy-34437343 


03:32:34.43 


-27:47:34.3 


29.3±0.4 


0.8±0.8 


-0.4+0.5 


0.20" 




UDFy-37706328 


03:32:37.70 


-27:46:32.8 


29.4±0.3 


0.9±0.7 


-0.2+0.4 


0.15" 





* Each of the sources here satisfy the color criteria for our Lyman-Break galaxy selections, but show slight detections in the Xopt 
images generated from the B435, Veoe, 1775, 2^850 data. These sources therefore do not make it into our Lyman-Break galaxy 
samples. See §3.3 and Appendix D for a discussion of Xopt '^^'^ our selection criteria. While a few of these sources will correspond 
to 2: > 7 galaxies, the majority will not according to our simulations (Appendix D). 

Lower limits on the measured colors are the Icr limits. Magnitudes are AB. 

Half-light radius calculated using SExtractor 
d References: [1] Yan et al. (2010), [2] Wilkins et al. (2011) 
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TABLE 20 

Possible 2:~7-8 candidates that did not meet our selection criteria.'' 



Object ID 


R.A. 


Dec 


Jl25 


Z850 — Jvl^ 




J125 — ifieo 




Ref<^ 






Possible z r^l candidates 










ERSz-2320541505f 


03:32:32.05 


-27:41:50.5 


25.9±0.1 


1.9±0.7 


1.2±0.4 


0.3±0.1 


0.28" 




ERSz-2226543006S 


03:32:22.65 


-27:43:00.6 


26.1±0.1 


1.1±0.3 


0.1±0.2 


0.4±0.1 


0.20" 


1,2,3,5 


EKSz-2279041042 


r^o oo 

03:32:27.90 


-27:41:04.2 


26.2±0.1 


2.7±0.8 


0.7±0.2 


— 0.6±0.1 


nil// 
0.11 


5 


ERSz-209754559lf 


03:32:09.75 


-27:45:59.1 


26.2±0.2 


>1.8 


0.7±0.5 


0.0±0.3 


0.24" 




ERSz-2030845122'''S 


03:32:03.08 


-27:45:12.2 


26.3±0.1 


1.1±0.4 


1.1±0.4 


0.2±0.1 


0.13" 




ERfaz-2104045407» 


03:32:10.40 


-27:45:40.7 


26.5±0.1 


0.8±0.4 


0.4±0.3 


0.2±0.2 


0.23" 




ERSz-1599542348'^ 


03:31:59.95 


-27:42:34.8 


26.5±0.2 


0.9±0.6 


0.7±0.6 


0.3±0.2 


0.27" 




ERfaz-2417341555= 


03:32:41.73 


-27:41:55.5 


26.6±0.2 


>1.2 


0.5±0.5 


0.4±0.3 


0.31" 




ERSz-2313142239f 


03:32:31.31 


-27:42:23.9 


26.6±0.2 


1.0±0.5 


0.8±0.5 


0.2±0.2 


0.21" 




ERSz-2338641089f 


03:32:33.86 


-27:41:08.9 


26.6±0.2 


1.4±0.8 


0.6±0.4 


0.1±0.2 


0.27" 




ERSz-2413843168S 


03:32:41.38 


-27:43:16.8 


26.8±0.1 


0.9±0.3 


0.3±0.2 


0.0±0.2 


0.13" 




ERSz-2308840450'' 


03:32:30.88 


-27:40:45.0 


26.8±0.2 


1.1±0.6 


0.8±0.6 


0.1±0.3 


0.22" 




ERSz-2379242208' 


03:32:37.92 


-27:42:20.8 


27.0±0.2 


>1.4 


1.1±0.9 


0.3±0.3 


0.24" 




ERSz-2202443342S 


03:32:20.24 


-27:43:34.2 


27.0±0.2 


0.8±0.5 


0.3±0.3 


-0.1±0.3 


0.19" 




ERSz-2346843052S 


03:32:34.68 


-27:43:05.2 


27.0±0.2 


1.7±0.9 


0.5±0.4 


0.1±0.3 


0.19" 




ERSz-2400642142' 


03:32:40.06 


-27:42:14.2 


27.2±0.2 


>1.5 


1.0±0.7 


0.2±0.3 


0.17" 




ERSz-2153943286S 


03:32:15.39 


-27:43:28.6 


27.2±0.2 


0.9±0.5 


— 0.0±0.3 


— 0.3±0.3 


0.13" 




ERSz-2232742361'^ 


03:32:23.27 


-27:42:36.1 


27.2±0.2 


0.9±0.5 


0.2±0.4 


0.0±0.3 


0.16" 




liRt3Z-242DD43478* 


r^o Ac\ e*f^ 

03:32:42.66 


-27:43:47.8 


27.2±0.2 


1.2±0.8 


>1.2 


0.3±0.3 


n 1 c// 
0.15 




ERSz-22293440998 


03:32:22.93 


-27:44:09.9 


27.3±0.2 


1.0±0.5 


0.1±0.4 


0.3±0.3 


0.14" 


5 


ERSz-2267841314f'' 


03:32:26.78 


-27:41:31.4 


27.4±0.2 


>1.6 


0.8±0.6 


-0.0±0.3 


0.16" 




ERSz-2270041428S 


03:32:27.00 


-27:41:42.8 


27.5±0.2 


0.8±0.5 


0.1±0.3 


-0.2±0.3 


0.15" 




ERSz-2202541197f 


03:32:20.25 


-27:41:19.7 


27.5±0.2 


1.3±0.9 


0.3±0.5 


0.3±0.3 


0.15" 








Possible z ~ 8 candidates 










ERSY-21462423691 


03:32:14.62 


-27:42:36.9 


26.0±0.1 


>1.7 


1.3±0.7 


0.5±0.2 


0.25" 




ERSY-2416741275f'™ 


03:32:41.67 


-27:41:27.5 


26.1±0.1 


>1.4 


>1.4 


0.6±0.2 


0.21" 




ERSY-2165043527'' 


03:32:16.50 


-27:43:52.7 


26.2±0.1 


>2.2 


1.3±0.5 


0.0±0.2 


0.21" 




ERSY-2376440061° 


03:32:37.64 


-27:40:06.1 


27.0±0.3 


1.3±0.9 


>1.3 


0.1±0.3 


0.41" 




ERSY-236754134lj 


03:32:36.75 


-27:41:34.1 


27.1±0.2 


>1.4 


>1.4 


0.3±0.3 


0.17" 




ERSY-2107343270f 


03:32:10.73 


-27:43:27.0 


28.2±0.3 


>1.1 


>1.4 


-0.3±0.4 


0.12" 





" Sources in this table either slightly miss our selection criteria or show slight detections in the Xopt images generated from the 
^4351 ^606) *775) •2^850 data. These sources therefore do not make it into our Lyman-Break galaxy samples. See §3.3 and Appendix 
D for a discussion of Xopt «md our selection criteria. 

^ Lower limits on the measured colors arc the Icr limits. Magnitudes arc AB. 

Half-light radius calculated using SExtractor 
" References; [1] Bouwens ct al. (2010c), [2] Hickey et al. (2010), [3] CastcUano ct al. (2010a), [4] Wilkins et al. (2010), [5] Wilkins 
ct al. (2011) 

° The size/structure of this source and its colors arc consistent with the source being a T-dwarf. The case is particularly compelling 
for ERS-2279041041 

' Show slight (> 1.5(t) detections in one or more optical bands, strongly suggesting they do not correspond to bona-fide z~7-8 
galaxies. For ERSz-2097545591, the detection is only apparent in a small 0.2"-diameter aperture, 
s If a star-forming galaxy, more likely at z~6-6.6 

^ This source shows a red Y — J ~ 1 color - and hence a possible z ~ 7 galaxy if the Lyman-BreaJc falls between the Y and J 
bands - but the source is also detected at >1.5it in the i and z bands. 

' Despite not satisiying our selection criteria, these sources may correspond to z > 6.5 galaxies. 
^ Small 2ct detection evident in a small 0.2"-diameter aperture in the i343B band 

la detection in i band and apparent \a detection in z band inconsistent with source being at z ~ 8 
' la detection in Veoe band in 0.2"-diameter aperture 

Red J — H 0.6 color suggests this source is likely a low-redshift, intrinsic red galaxy 
" Slight detection in 4775, zsso band in 0.2"-dianieter aperture 



